For Reference 


NOT TO BE TAKEN FROM THIS ROOM 


1IBRIS 
NISERSTCACIS 
i RCALNSIS 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Alberta Library 


https://archive.org/details/Martin19/72_ 0 


a 
b pi, ba? 
idl) 


‘ 
eit 2, is “chp alas Ve renege 
i i ot ee. oo oh, Pir F ra) i 
ba ism, ' 7 ii i 8 a : v i a i { , i ; Hi) “ 


a 
a q 


t 


TP 
beter 
4 


oy Vy “aa i Aes 
es aa 


THE UNIVERSITY OF ALBERTA 


A SIMULATION APPROACH TO 


COST ANALYSIS OF FOREST FIRE DETECTION SYSTEMS 


by 


GREGORY ARCHILLE MARTIN 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 


OF MASTER OF BUSINESS ADMINISTRATION 


EDMONTON, ALBERTA 


SPRING, 1972 


Lie 


ABSTRACT 


A simulation model for the evaluation of alternative 
forest fire detection systems is beseri beds An analysis of 
nine years of fire-weather data and nine years of fire data for 
the Footner Forest in Northern Alberta was used in the develop- 
ment of the simulation. A three variable fire growth model is 
developed and described. A Monte Carlo sampling technique is 
used to generate the values of stochastic variables. Chi- 
square and Kolmogorov - Smirnov tests were used to substantiate 
probability distributions for significant variables. As a result 
the simulation was judged to be realistic with respect to tower 
and air route operation. A sequential stopping rule was used 
to reduce the number of simulated fire seasons by testing the 
mean of the total seasonal cost at a 95% confidence level. An 
example of how the output from the model can be used to show the 
relationship between various design configurations and the total 
cost of fire detection and fire suppression is presented. Also 
presented is a run-write-up of the model complete with procedures 


for modifying the design and the input data. 
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CHAPTER 1 


THE MODEL: DESCRIPTION, VERIFICATION AND RESULTS 


Introduction 
This study was undertaken as the first step in a series of 
research projects on the prevention - detection - suppression - refo- 
restation systems managed by the Alberta Department of Lands and 
Forests. The primary objective of the study about which this paper 
reports was to develop, refine and adapt a quantitative decision 
model to the realities of an operating forest management service 
so that current techniques may be assessed in terms of their useful- 
ness in improving forest fire detection. 
The current state of technology is such that the design 
of fire detection systems revolves around the effective use of 
three major system components; these are manned tower, surveillance 
aircraft, and the general public. A number of attempts have been 
made to incorporate the use of Management Science techniques into 
the area of forest fire detection [5] [6] [12]. The work presented 
here differs in several important aspects: 
(1) The model developed is based on actual fire and 
weather data which was accumulated by the Alberta 
Department of Lands and Forests for the past nine 
years. This large volume of data was used to deve- 


lop probability distributions for use in a Monte 


Carlo decision model by utilizing the system 
of computer programs provided by SPSS [10]. 

(2) The fire-growth model which was developed was 
based on published fundamental-empirical research 
into various forest environmental parameters and 
the resulting effect on the burning of forest 
fuels [2]. The model considers the propensity to 
burn of fuels in the forest, the burning charac- 
teristics of four different fuel types and diurnal 
fluctuation of fire-weather factors!, This 
approach differs from the works in which analytical 
fire-growth models are developed by fitting curves 
fOr datas tromeactual fire histories. [12]. 

(3) The model aids the detection-system designer by 
providing information about the effectiveness of 
particular design configurations with respect to 
fire discoveries and about the overall design 
performance as measured by total system cost. The 
model is specifically aimed at the search for a detec- 
tion design that produces the lowest total cost of 


the following criterion function: 


'1l. The idea for the fire-growth model was germinated by some 
excellent unpublished work of R.S. Miyagawa and E.V. Stashko 
of the Forest Protection Branch of the Alberta Forest Service. 
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(4) In addition, the advice of Mitroff was heeded and 
a great deal of thought was given to how the work 
should be validated [7]. Given the previous work 
done in this field together with the ability to 
analyze what was judged to be very good empirical 
data, it was felt that any modeling should pass 
generally accepted statistical tests for goodness 


Of LIt. 


A Description of the Model 


The Footner Forest of Alberta was chosen as the prototype 
for a model which could eventually include forests throughout the 
entire province. Footner Forest consists of approximately 32,000 
square miles of water covered and forested lands containing coniferous, 
mixed wood, lichen-muskeg, and grass fuels. The northern boundary 
of the forest adjoins the Northwest Territories and the western 
boundary adjoins British Columbia. The Footner Forest is currently 
protected by a formal detection system consisting of 19 manned towers 
and 5 air routes, which are scheduled on both a seasonal and a daily 
basis, and an informal public system which includes forestry officers, 


forest personnel, other aircraft, licensees, government officials, 


2. At the present time data adequate to support the cost of damage 
is not available; therefore the cost of damage has been excluded 
in the example of results presented in this thesis. 
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the general public, and others. The 19 manned towers are composed 
of 14 within the forest and 5 outside the forest but close enough 
to discover fires within the forest boundaries. 

The simulation model uses the Monte Carlo technique [8]. 
Experiments have been based on historical data used in a static 
fashion. For example, no attempt has been made to forecast 
changes in such parameters as the average number of fires or the 
average number of fire days per season. However, the design of 
the model is such that it can be readily changed to handle fore- 
casts as well as historical data if the designer so desires. 

The model works in the following manner. A fire season 
begins. The total number of fires expected during the fire season 
is computed by the Monte Carlo procedure, i.e., during computation 
uniform - pseudorandom numbers are generated and compared with an 
empirical probability distribution in order to determine values 
assumed by a variable. These fires are then spread throughout the 
days of the fire season in Monte Carlo fashion. If one or more 
fires occur on a given day then a location, cause, fuel type, and 
time-of-day is assigned to each fire in a similar manner. Daily 
fire-weather factors and visibility are also obtained by the Monte 
Carlo process. Given the visibility a check is made to see if the 
fire can be spotted from a fire-tower or from a surveillance air- 
craft. A fire might be discovered by either a tower or a surveillance 


aircraft or both. When this is the case the time that the fire burns 
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until discovery is taken as a minimum of the elapsed time between 
fire ignition and cor. A fire might not be discovered by 
either fire-towers or aircraft. In this instance, the length of 
time that the fire burns until discovered by a public agency is 
determined by the Monte Carlo Beceedutes After having established 
the length of time that the fire burns prior to discovery it is 
possible, using the fire-growth model, to compute the size of the 
fire at discovery. The fire. suppression aout is a function of the 
size of the fire at the time of discovery and is determined by the 
Monte Carlo method. This general procedure is followed for each 
day during the fire season. The number of fire seasons examined, 
during an experimental run under a given set of design conditions, 
is determined by the sequential sampling - stopping rule described 
by Naylor [8]. Once the conditions of the stopping rule have been 
satisfied, results of the experiment are printed and the computer 
run terminates. 

The fundamental logical framework supporting the model is 
straight forward and is portrayed by the simplified flowchart of 
Figure 1-1. Although the logic is relatively simple the functional 
relationships between variables and the interaction of variables 
in the computer program is complex. A summary of model relation- 
ships is presented in Table 1-l. 

Armed with this model the designer of forest-fire detection 


systems retires to his drawing board after having acquired a map 
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FIGURE 1-1 


GENERALIZED FLOWCHART 


GENERATE DAILY FOREST FIRE-WEATHER FACTORS 


READ DATA 


GENERATE NUMBER OF DAILY FIRE OCCURRENCES 


S THERE A FIRE TODAY? NO 


| 


YES 


GENERATE FIRE LOCATION 


GENERATE FIRE CAUSE 


GENERATE TYPE OF FUEL 


GENERATE TIME OF DAY OF FIRE OCCURRENCE 


GENERATE DAILY VISIBILITY 


COMPUTE TIME TO FIRE DISCOVERY 


COMPUTE FIRE SIZE AT DISCOVERY 


GENERATE COST OF FIRE .SUPPRESSION 


NO TS eS ei eS ENDVOFSTHE FIRE SEANSON2 


YES 


TABULATE SEASON'S DESIGN STATISTICS AND COMPUTE 


TOTAL COST FOR THE FIRE SEASON 


STOPPING RULE CHECK 
S THIS THE END OF THE EXPERIMENT? 


NO YES 


COMPUTE AVERAGE TOTAL 
COST FOR ALL FIRE 
SEASONS IN SAMPLE 


INITIATE ROUTINE 


FOR REMOVAL OF 
AUTOCORRELAT ION 


PRINT TOTAL COST 
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The Most Important 
Probabilistic 
Variables 
Vv 


1. Number of Fires 
Per Fire Season 


2. Number of Fires 
Per Day 

3. Spread Index 

4. Buildup Index 

S$. Burning Index 


6. Fire Location 

7. Fire Cause 

8. Fuel Type 

9. Fire Time of Day 
10. Visibility 

11. Fire Discovery By 


Tower 


12. Fire Discovery By 
Aircraft 


13. Time of Burn If 
Discovered By Tower 
or Aircraft 


14. Time of Burn If 
Discovered By 
Public Agency 


18. Fire Size At 
Discovery 


16. Cost of Fire 
Suppression 


TABLE 1-1 


Number Of 
Possible Values 
Variable Might 

Assume 


Average of About 
50 


Upper Bound of About 
7 


Number of Towers 
+1 


Number of Air Routes 
+1 


Upper: Bound of 
Approximately 
12 Hours 


Upper Bound of 
30 Hours 


Stratified Into 
15 Ranges 


11 


SUMMARY OF MODEL RELATIONSHIPS 


Number Of 
Probability Distributions 
Used To Model Variable 
Behavior 


None 


None 


None 


15 


Independent Variables, x, 
Used To Define The 
Functional Relationship 
VefiGor. 5 hox),) 


Nine Year Average 
Seasonal Occurrence 


Nine Year Averages of 

Fires Occurring on A 
Specific Day 

Specific Day In Fire Season 
Specific Day In Fire Season 


Day, Month, Spread Index, 
Buildup Index 


Month, Forest District 
Month, Forest District 
Month, Forest District 
Cause 

Spread Index 

Tower Location, Blind Area, 
Visibility, Fire Location, 
Time of Day of Fire Start 
Location of Air Route, 
Visibility, Fire Location, 
Time of Day of Fire Start, 
Time of Day of Flying 
Minimum of Time To 
Discovery By Tower or 
Aircraft 


Fuel 


Fuel, Burning Index, Time 
of Day of Fire Start 


Fire Size At Discovery 
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of the forest, a protractor, a ruler, and a piece of onion-skin 

graph paper. While at the drawing board the designer uses his 
judgment, analytical See intuition, imagination, and experience 
in creating a layout of towers and air routes that he thinks might 
provide appropriate detection coverage for the forest. He uses 

his ruler to obtain the map coordinates of each tower location and 

to measure the length of each leg of each air route. He uses his 
protractor to measure the angle that each leg of an air route makes 
with the horizontal, see Figure 1-2. If a "seen-area'’ map is avai- 
lable the designer uses the piece of onion-skin graph paper to 
estimate the amount of blind area within each of the twenty segments 
that make up the maximum observable area of a fire tower, see Figure 1-3. 
If a "seen-area'' map is not available blind area data is obtained from 
field reconnaissance or from a topographic map. 

Physical characteristics with respect to location of towers 
or air routes are not the only design factors which the designer 
might wish to investigate or control. For example, the designer 
might wish to establish policies with respect to actual operation 
of the system. For instance he might wish to fly his air routes 
once, twice, or three times a day or he may simply wish to vary the 
starting time of air surveillance. The model allows the designer 
to choose from a reasonably realistic set of management policies, 
or decision rules, and to investigate the results obtained through 


the implementation of a chosen combination of policies acting in 


FOREST MAP 


@ = Angle of Air Route Leg 
With Horizontal 


Length of Air Route 


Pat es ea Leg 
Air Route dearer ee / 


ON Aes Location 


(x,y) 


x 


(0,0) 


Reference Point For 
Map Coordinates 


FIGURE 1-2 


EXAMPLE OF DESIGN LAYOUT 
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FIGURE 1-3 PARTITIONING OF FIRE TOWER OBSERVABLE AREA 
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conjunction with a given design. Decision rules with which the 
designer might experiment are shown in Table 1+2. After preparing 

a design and selecting a set of decision rules the designer uses 

the model to experimentally investigate the behavior of his creation. 
The computer output produces an average total cost of fire detection 
and fire suppression for any given physical design - decision rule 
combination. In addition, information is provided with respect to 
the effectiveness of tower placement and air route location. This 
provides the designer with guidelines for layout improvement. Also, 
the output is presented in such a manner that statistics with respect 


to important variables are readily discernable and easy to analyze. 


Goodness of Fit 

As mentioned previously the Footner Forest is presently 
protected by a formal detection system consisting of 19 manned 
towers, 5 air routes, and a diverse public system. The approach 
used to test the model was to run it under what was understood to 
be those management operating policies which produced the nine 
years of data used to develop and operationalize the model. The 
simulated output was compared with the historical data to test 
the hypothesis that the set of data generated by the model had 
the same frequency distribution as the set of historical data. 
After this verification process, the model was then used to 


investigate other design configurations. 
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TABLE 1-2 


DECISION RULES AVAILABLE TO DESIGNER! 


Fly each air 
Fly each air 
Fly each air 
prescribed. 


Fly each air 


route once each day during the months prescribed. 
route twice each day during the months prescribed. 


route three times each day during the months 


route only if the burning index is equal to or 


greater than a level prescribed by the designer; for example 


fly only if the burning index is greater than or equal to 2. 


Fly each air 


route only if a fire has been caused by the 


agency identified by the designer; for example fly only if a 


fire has been caused by lightning. 


Fly each air 


route if the prescribed burning index is equaled 


or exceeded or if a fire has been caused by a prescribed agency. 


The model has been dimensioned to handle from 0 to 30 towers 
and from 0 to 15 air routes or any combination thereof. Air 
routes can be coded so that they will be active only during 

prescribed months of the fire season. In addition the daily 
flying schedule is prescribed by the designer. 
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It was judged that the most important factor with respect 
to daily operation concerns the decision made by the Fire-Control 
Officer as to whether or not to operate the air routes. He uses 
his judgment and experience with respect to interpretation of 
fire-weather conditions in the forest in making the decision. If 
forest conditions are extremely hazardous, or if settlers have 
been clearing and burning land, or if lightning storms have moved 
through the forest he will probably operate the air routes. It 
is this decision process against which the model was tested. Given 
the functional relationships described in Table 1-1, it was judged 
that the most important variables to test for goodness of fit were 
cause, fuel, visibility, fire size at discovery, fire suppression 
cost, and overall frequency of discovery by towers, aircraft, 
and public agency. The output was designed for easy analysis of 
these variables. It is worth pointing out that goodness of fit 
tests for the cause, fuel, and visibility variables are primarily 
a test of data, while the goodness of fit tests for fire size, fire 
cost, and overall frequency of discoveries are essentially a test 
of the model. 

The results of both Chi-Square and Kolmogorov-Smirnov 
tests for goodness of fit are presented in Table 1-3. It is noted 
that the overall frequency of discovery by towers, aircraft, and 
public agency did not pass either test at the 95% confidence level. 


The actual and simulated overall discovery frequencies are shown 
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TABLE 1-3 


RESULTS OF STATISTICAL TESTING OF SIMULATED DATA 


AT THE 95% CONFIDENCE LEVELL 


Variable Passed . «Variable Passed. 

Variable Chi-Square Test? Kolmogorov-Smirnov Test? 
Cause No Yes 
Fuel Yes Yes 
Visibility Yes Yes 
Fire Size at Discovery No Yes 
Fire Suppression Cost Yes Yes 
Overall Frequency of 

Discovery by Towers, 

Aircraft, and Public 

Agency No No 


1. The sample sizes were 2649 for Fuel and 2667 for all others. 
The design alternative used was 19 towers and fly if the burning 
index is greater than or equal to 2 or if a fire has been caused 
by lightning. 
TABLE 1-4 
COMPARISON OF HISTORICAL AND SIMULATED FREQUENCY OF DISCOVERY BY 
TOWERS, AIRCRAFT, AND PUBLIC AGENCY 


Simulated Frequency 


Historical Frequency of Discovery 
of Discovery Fire Sample Size - 2667 
Towers 41.8 40.3 
Aircraft ere, TG 7. 
Public Agencies 3959 43.0 


Total 100.0 Total 100.0 
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in Table 1-4. It was concluded that the difference between these 
two frequencies could be explained in the aircraft scheduling 
procedure and the aircraft flight paths. It would appear that 

the aircraft routes actually find more fires than what is predicted 
by the model because of the pilot's flexibility in his course direc- 
tion and his current information as to where fires may exist based 
on the route he flew yesterday or storm movements from the previous 
night or settlement activity. Some of the difference could also 

be due to some behavioral aspects of fire reporting about which no 
empirical evidence exists. Nevertheless, the overall goodness of 
fit of all but one of the most important variables of the model, at 
a large sample size (N greater than 2000) and the visual closeness 
of the actual and simulated frequencies of discovery attests to the 


validity of the model. 


Example of Results 


The model was tested using data from the Footner Forest. 
A total of 18 experiments were performed under various design- 
policy combinations. A typical computer run using an object deck 
requires between 72 and 99 seconds of CPU time on an IBM 360/67 
computer and costs about $40.00 at current commercial rates. The 
result of the computer runs is presented in Figure 1-4. Detection 
cost includes charges of $6,000.00 per year per tower and $48.00 


per flying hour. Fire suppression cost includes the total of fire 
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Results of Footner Forest Analysis 
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identification of potentially suitable detection strategies. The 
model's only purpose is to provide decision makers with a method 
for easily evaluating various detection Est bo ee Prin erabres, 
the results cannot be directly related to the detection strategy 
employed in another forest. However, it is important to realize 


that the model can be easily adapted for use in different forests. 


Accordingly, the implications of a large number of different detection 


strategies can be readily explored in many different forests. 


Conclusions 

The results strongly suggest that a quantitative decision 
model may be successfully refined and adapted to the realities of 
an operating forest management service. It may be further concluded 
that a quantitative decision model may be useful in improving the 
art of forest fire detection. However, this later conclusion must 
be more fully considered. A well designed program of implementation 
must be undertaken to assure success in operating the model. Such 
a program of implementation would be directed toward (1) formally 
evaluating the usefulness of the model in terms of improving the 
art of forestry management and (2) identifying those factors which 
are important for the successful introduction of such a model in 


an operating forest management service. 


4, A paper by P.M. Maher, B.M. Harnden, and G.A. Martin has 


been prepared in order to describe such a program of implementation. 
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Further steps in the quantification of forest fire 
detection and suppression are judged to be era First, the 
importance of continuing the work being done to establish values 
for public lands is clear [1] [11]. These efforts may be furthered 
if serious attempts are made, with the aid of the decision model 
described in this thesis, to examine the implications of various 
evaluation strategies. Second, the unexpected results obtained 
from the analysis of the existing detection system in the Footner 
Forest would suggest that an analysis of the fire-prevention sub- 
system be undertaken in order to examine the relationships between 
this subsystem and the public agency portion of the existing fire 
detection systems. Finally, now that a comprehensive model is 
available which facilitates a better understanding of the overall 
system, it may be possible to circumvent some of the complexities 
of the present model to develop useful analytical models which 


may lead to identifiable optimal detection strategies. 


4 ok? Sess ee tobte 


yea Nhe a oe i Tee ae a eee lly. 
r y : vf j +. nS a. ; eile ad ' 
OA ada 2 capa , \igiid ibaa? a8 a is oo tees oo ve” bits nok: 
oe | i hk hat i 
» one : i ee arp pips ‘Sat : ; ik! Ms ers, 
een . @einy iieBtéeres ‘Or wt ‘stow allt at ton “hic sien oie 
et | Pky Aa i a a oy Mi a Ak ) ao a 
; beset wut od Yer sy7ol is “shea fen ir) “nals eb at = vip 
.  £sbom nekakosh ‘sis i abe ait sabe sabe te a9 
ms A 7 a) ra a 


ee 4 
abot ey to a rot seoEte mi otle ‘uthaxs: os 1 peeite ane! nt baddi-x oeab 
: ~ Ae F ae oe Sat 


4 a ys 4 

cs 7 a a : sae i a a, 

bsetsado eatuegs be j98q sin ‘sda ba0998 - 9eGhgpssiia norton sited 

; : = aa aa) _ j 

| Dims Cpe. mt Vireo 
tsatoc% silt at mageye re): kisoaat a nye odd id eka isa + hanes ve 
; mar iS) Stated Le i Ss 

y var ‘ol : 


& 


fie AG anavarahatt®, and to. ateyies 


meawied eieusust’i siox eet atm Rg ae iginmiacd ‘“ tel 

a en 

: ay ae ‘es 
bréoq any, payee a a 


9tii antiaions 53 to “a0 t3z0q) ORS gs: | 
_) — 
ee 


‘ 


is . “ar 7 t = 
! 2h feboo subtest adsaqiap5' B/ ids wor , eehiaaed eee sak : a 
ane se ll ee 
Veen a 
flevevo sa 43 gatbuste bathe teased s eosastntoes ‘dokdw atte Ee 5 
( ve \ : ma ie » - 
. itixelgqnion sis 26. soe thawantoryt of agtaety ¢ ad yam tm aie te eae 
M ) 34 > . : 
dotitw gist bom Teoks elsns. ‘fuses a oLaveb ‘ca slow 40 


A ay = ost 


7 | Py aiskge3 sexta) dokioasab teatsg9 Vecatesaees » a 


i pa ee , 
7 EN aaa > ince a) Ma Dh 
Sy i: a 
; 
Ad 


—e 


_— 5! 


' 5 7 


| yg 


20 


References - Chapter 1 


1. Department of Regional Economic Expansion, "The Canada Land 
Inventory Objectives Scope and Organization", Report No. 1, 
Second Edition, 1970. 


2. Fishman, G.S., "Estimating Sample Size In Computing Simulation 
Experiments," Management Science, Vol. 18, No. 1, September, 
LOT Leaps. 2s 30% 


3. Hardy, Charles E., and Franks, James W., "Forest Fires In Alaska," 
U.S. Forest Service Research Paper INT-5, Inter-mountain Forest 
and Range Experiment Station, United States Department of 
Agriculture, Ogden, Utah, 1963, pps. 41 - 48. 


Aw elavs.e Wels. and Winkler, R., Statistics, Vol. 11, Holt, 
Rinehart and Winston, Inc., New York, 1970. 


5. Kourtz, P.H., and O'Regan, W.G., "A Cost-Effectiveness Analysis 
of Simulated Forest Fire Detection Systems,'' Hilgardia, Vol. 39, 
No,wal2.vAucust, 1968. 


oF ; -A Model for a Small Forest Fire .... To Simulate 
Burned and Burning Areas for Use in a Detection Model," 
BOReStoC1ence,, VOL. 1/,) NO. 2, 19/1, pps. 163. = (169. 


7. Mitroff, Ian I., "Fundamental Issues in the Simulation of Human 
Behavior: A Case in the Strategy of Behavioral Science," 
Management Science, Vol. 15, No. 12, August, 1969, pps. B-635 - 
B-649, 


8. Naylor, Thomas H., Computer Simulation Experiments With Models 
of Economic Systems, John Wiley & Sons, Inc., New York, 1971, 
pps. 293 - 294. 


9. Naylor, T.H., and Finger, J.M., “Verification of Computer 
Simulation Models," Management Science, Vol. 14, No. 2, 
October 1967, pps. B-92 - B-101. 


10pem ites Noben Bent. -D. Ha wand Bull. (C,H... Statistical Packsce: For 
The Social Sciences, McGraw-Hill, New York, 1970. 


Ll? Nowicki, J.J. and White, S.R., et al., ‘Methodology for the 
Foothills Resource Allocation Study," Canada Land Inventory 
Program, Forest Land Use Branch, Department of Lands and 
Forests, Alberta, March, 1971. 


3 bese 


iS 


Zu 


Parks, George M., "Development and Application of a Model for 


Suppression of Forest. Fires, " Management Science, Vol. 10, 
No. 4, July 1964, pps. 7/60 —" 760. 


Van. Horn, Richard. L., "Validation of Simulation Results," 
Management Science, Vol. 17, No. 5, January, TOV Ie ppS.8 2 4p 220. 


oe i: i" 


a a 


9 bean, a 


. i 
" wows 


tay nS aT 


a ea ea 
i" me a ane 


:) 
m p> we EY WW aie 
wl 4c 


; | mn. 
1 Wein we me ie eal hig 


* art. 
ey sa | 


AE 
eh Uae : 
uy 


fox Aisle : 


a 


pit, 


a 


49 
is i 
. hia) y 


is Om oF aA 4 


bak aie cae all 
oeel ho 


\ we r’ *, ( (oPaes ne) 
marae Ne ‘wee eta 


viet ee 


ae Bind peer 


ov ase 


hia 


Tet 


igi 


| Bik: i ae Sih, oa 


i a. 


» 7 re on 


4 ee i nue anee.) i 


ic pee ira tee 


oe 


Pie 


serine 


as , 
a a. 
: Lal — 

. y 


oe a 
Wy 


ad, 


ee 


nt : 4 ee ©. a. 7 
x | Se. ae 
eee reg wees <4 fila: 3. 
Ree ant Oe 


22 


CHAPTER 2 


BACKGROUND TO THE STUDY 


plvereontat ge Hesrddy paste? 

The Faculty of Business Administration and Commerce at 
the University of Alberta in Edmonton has acquired a number of 
faculty members who have the interest and necessary technical skills 
to play an important role in the development of Canada and in parti- 
cular the Province of Alberta. In an attempt to realize the poten- 
tial of this new faculty, a series of meetings was held with 
interested business and government organizations. 

One series of meetings was held within the Department of 
Lands and Forests of the Province of Alberta. These meetings 
culminated in a research agreement between the Department of Lands 
and Forests and the Faculty of Business Administration and Commerce. 

The research agreement was an open ended agreement in that 
the department agreed to support an unspecified research project. 
As a result, a number of strategy meetings were held among key 
personnel of the Department of Lands and Forests and the Faculty 
of Business. The purpose of these meetings was to identify key 
personnel that would be interested in participating in a study, to 
set out some general guidelines as to what the project should accom- 


plish and to identify possible areas open to study. 
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During the period from April 26 to May 20, 1971, 
G.A. HecEaAS a graduate student in the Faculty of Business Adminis- 
tration and Commerce, gathered information on all aspects of the 
Alberta Forest Service Fire Protection Branch in an effort to define 
the area of research. Sources of information were interviews with 
Mee (C. he @iact, supervisor Orethe Pire Control oection; Mr% CS. 
McDonald, Fire Control Officer; Mr. W.J. McLean, Supervisor of the 
Weather Section; Mr. J. Niederleitner, Fire Control Technician; and 
Mr. G.H. Deans, Fire Control Dispatcher. 

Interviews were also conducted with Mr. S.R. Hughes, Head 
of the Forest Protection Branch and Mr. H.M. Ryhanen, Head of the 
Administration Branch of the Alberta Forest Service; and Mr. J. Griegel, 
Research Officer with the Canadian Forestry Service in Edmonton. Also 
included as information sources were numerous papers and studies 
located in the Department of Lands and Forests Library. Three more 
meetings were held in the Spring of 1971, on April 27, May 5 and 
May 19. Professors B.M. Harnden and P.M. Maher along with Messrs. 
D. Cauvin, Section Head of the Economics Section of the Timber Branch; 
R.S. Miyagawa, Supervisor of Fire Research for the Forest Protection 
Branch; and G.A. Martin met at the University of Alberta. The purpose 
of these meetings was to discuss projects applicable to the Alberta 
Forest Service. 

It soon became apparent from these meetings and discussions 
that there was a high degree of interdependence among the different 


projects being considered. For example, analysis of a suppression 
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System cannot be done without prior knowledge of the detection system 
and the size of the fires at the time of initial attack. A second 
area of concern involved the long term implications of the research 
ncn cets Ideally the project sought was one that would lay the 
ground work upon which future research could be based. 

On the advice of Messrs. Miyagawa and Cauvin, a proposal 
was drafted on May 25, 1971 and recommended by Mr. S.R. Hughes and 
approved by Mr. R.G. Steele, Director of Forestry, shortly thereafter. 
The purpose of the proposal was to outline a sequence of research 
studies that could be carried out by the Alberta Forest Service. 

The primary objective set out in the proposal was to provide the 
management of the Alberta Forest Service with a clearer picture of 
alternatives available for strategic decision making. The proposal 
which was recommended began with a study of fire incidence and was 
to proceed through the sequence shown in Figure 2-1. 

Figure 2-1 represents a schematic of the possible research 
studies recommended in the May 25 proposal and shows the order the 
studies should assume. The schematic begins with a study of forest 
fire incidence. This study could be followed by either a study on 
fire prevention or fire detection, The output from the historical 
fire incidence study would be required as input for both the prevention 
and detection study. Output from the detection study would be required 
for the proposed presuppression study and output from the presuppression 


study could be used in the suppression study. A separate study of 
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suppression forces would provide input to the presuppression study. 
The overall administration system study would be done last. Output 
for the other six studies would provide information required to 
conduct this last ae 

The previously mentioned problem of interdependence among 
the 7 proposed studies was minimized by organizing the planned sequence 
shown in Figure 2-1. The sequence was arranged so that the results 
of one project would supply some of the information required to solve 
the next scheduled assignment. Accordingly, as one study is finished, 
the next is brought into much sharper focus because of the new infor- 
mation supplied from the completed assignment. 

The research began with a study of fire incidence in Alberta. 
The objective of this study was to arrange and program existing fire 
data into a system whereby meaningful analysis could be done on 
forest fire causes, occurrences, growth, etc. Under the guidance of 
Mr. R.S. Miyagawa, systematic analysis of Alberta's fire and weather 
data from 1961 to 1969 was undertaken. The data was analysed using 
an IBM 360-65 computer and the Statistical Package for the Social 
Sciences (SPSS). In addition, two field trips into the Whitecourt 
and Rocky-Clearwater Forests were made to gather information on how 
the fire control and detection systems operate in Alberta. The 
study of fire incidence was essentially completed once the data was 
organized and the computer was providing statistical information on 


Alberta forest fires. The formal fire incidence study was completed 
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by mid August and Miyagawa is currently updating the data bank 

with fire data from the 1970 and 1971 seasons and extracting meaningful 
statistics on Alberta forest Pinca The formal research proceeded 

from the incidence study to a study of a fire detection system. 
Professor BLM. Harnden and G.A. Martin began programming a simulation 
model for the detection system in the Footner Forest in Northern 
Alberta. The objective was to develop a general purpose computer 
simulation which could eventually be used to study the detection 

system for all of Alberta. The results of the detection study can 

be used as a guideline in allocating resources and measuring the 


effectiveness of the present detection methods. 


Pre-Requisites to the Model Design 


The accumulation and organization of the data required to 
operationalize a large simulation model is a very time consuming task. 
Figure 2-2 shows the sequence of main events leading to actual program- 
ming of the simulation. Prior to the beginning of this research, the 
Alberta Forest Service had accumulated fire data for each of the 4,668 
fires occurring in Alberta between 1961 and 1969 inclusive. This data 
had been transfered from fire report forms onto IBM cards and magnetic 
tape for use on an IBM 7070 computer. There were 67 different pieces 
of information about a fire available from the fire report forms. A 
list of this information is shown in Table 2-1. The first step was 


to transfer the data from the IBM 7070 tape to a nine track tape for 


, yi alee © { 
“es ce SL 
; i & i) ‘s ae, " F 
\ 1 “Fi Me : mA, , A 
18) 8 an my 
o Dy te. ve )4 i v 
t} : M : : na oy Ay 
1 ey = a 
o a a 
Sins BIBS 3a gnitebqu [sis tIus 
| A ee 
iuigninsa aiionsaits bes i at ; 
| . Bad eu ha ches vn ‘ 
bobosso Goltenees Lemitot, bs dl Bod si 833 a 
tte shee a mob3se39b out? 1 & a sabes 5 od ‘ehuse: sonobto ant ott 
ae f % a i op Lit rw AP | 
xotssiume & grime rgo tq m8; aed bart As a bas steer ana | roses 
i Nit mr . 
rads 10" at saan? wanyoo! sift: nt ‘maseue, noksosash sit ww? t 
} ‘ 7 ‘ cae : sh ne > : 
tesigmos seoq we Kewanes s 5 alowsb oF BBW , avitastde ont. A sented sey 
eaat ¥ a et Lay 7 i 


ro rt30835b ads yhlise a3 bowe ad = sudo v5 aad Hoke insantin 


ABD vhid a ‘qobsdeteb ert a epiuast “gar 


ye ay es pas ; - 
- sniisbtng & es oe 


a3 id Hino bas | Lasioiuods ants 


pet) ee : 7 ou ©) 
-aboritem aoktoss abbidinsesxd wits 26 aauae 
fa , a . q 
: = ae a? | ici 
- , PL is J A don 
Rr ens agtesd Tebow aid oF ead tetunerics 


fi af 


13 betiupss sisb saz to soktssinagso bas “ois of umudoe sa? 


vised gatimenca ombs vrav 8 sek tsbom notsslombe satel, fe ssninot en a ‘ 


AX 1 Ms 


“m8 DBO Isosoe ot. gatbeo! eansvs se aaa to sanoupee » ‘of evs Sf ica : 


. a 
ott Prrerre atta 2 gainntsed ads ‘od zoxt nobis bebidas sitg Jo-gitke 
BOG, A ory to a 2 103 sonb on eset bert ‘sabiaee abt08 sonecth ‘ 7 : 
stab aldT svteubean caer bre Taet asowsd eansdiA eet ‘ ie — 
r seahecied bars, abut. wax we ste ee. ida 


+a 
soatir sors ‘a oe mt Hqmo> a am it eee fit 1 a 


- 


ae t i 


a 
o 


ci 


» how Ud 


ae io 


ee qoie seth a as 


i oe | i eof x 
oa ‘i es an ny 


oa) 
i 


y a a 


A eaieto® sroaes of} a a i . 


FIGURE 2-2 28 


A TIME SEQUENCE OF KEY EVENTS 


KEY EVENTS PRELIMINARY & INTERMEDIATE 
REQUIREMENTS 
APRIL 
Research offer 2 
Preliminary meeting 13 Broad problem areas outlined 
where research may be applicable 
Offer accepted 1s 
Research started 26 
2nd meeting of U of A staff Start of literature survey 
and forestry officials 27 and indoctrination period 
MAY 
3rd meeting of U of A staff 5 Initial research proposal writ- 
and forestry officals ten and reviewed by Miyagaws and 
Cauvin. Proposal nodified 
4th meeting of U of A staff 19 Modified proposal reviewed and 
aod forestry officials _ remodified. Final draft of 


research proposal written for 
submiesion to Alberta Forest 
. Service, approved by R.G. Steele 
Research Proposal submitted 25 
to Alberta Forest Service 


FIRE INCIDENCE STUDY BEGINS Introduction to IBM system and 
review of existing tape files 
Organization of fire data on 28 and data. Reports written to 
TBM 7070 tape begins explain and request computer usage 
and time required. 
JUNE 
Computer usage report issued 2 Fire data reorganized and examined 
also time requested for IBM- 
360 
Fire data transfered to IBM- 3 SPSS program being written to 
360 tape create an SPSS file "Forstudy” 
Computer programs to create 18 Literature was studied to decide 
SPSS file written which variables to select in a 
simulation model for detection 
study 
JULY 
SPSS file FORSTUDY created 15 Cross tabulations were programmed 
to aid in determination of variables 
needed for detection study. 
Organization of weather data 16 Weather data investigated and 7070 
for transfer to 360 tape begins tape reorganized. Intermediate 
programs to transfer weather data 
by forest used to reduce computer 
requirements. 
AUGUST 
Footner Forest data segregated 
onto an IBM-360 tape. ‘ 
SPSS file WEATHER created for 12 Cross tabulations were programmed 
Footner Forest weather data to accumulate probability distri- 
butions on fire variables and 
weather variables. 
FIRE DETECTION STUDY BEGINS Plow chart drawn for a detection 
model using a simulation technique 
SEPTEMBER 
Pootner Forest maps cxamined 
Start of programming of 3 Miyagawa supplied consulting 
simulation model in FORTRAY IV service 
by C. Martin and B. Harnden 
Coding and debugging of program 
MARCH 
1972 
Footnec Forernt detection 
design runs completed 
APRIL 
1972 


Martin's Masters thesia and 
report for Alberta Forest 
Service completed, 
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usage on the IBM 360-65 computer located in the Provincial Data 

Center in Edmonton. The data was further processed and transfered 

to the IBM 360 tape by the end of June. Using the Statistical Package 
for the Social Sciences (SPSS) available at the Provincial Data 
Processing Center, the data was stored at the Provincial Data 
Processing Center in an SPSS file called FORSTUDY. 

As well as the file FORSTUDY, a file called WEATHER was 
created. The WEATHER file contained the Footner Forest weather data 
from 11,246 weather reports. Each weather report contained information 
on 26 variables. A listing of the 26 variables is shown in Table 2-2. 

SPSS was selected because of its availability. It is under- 
standable to people unfamiliar with computer programs and is designed 
to accomodate ordinal and cardinal values. This is an important 
aspect since many of the variables involved in forest fires are assigned 
ordinal values or categorized by qualitative attributes. SPSS allowed 
the examination of data pertaining to each forest or ranger district 
independently. The subroutine CROSSTABS was used extensively to 
examine the relationships between key variables and to obtain frequency 
distributions from which probability distributions could be calculated. 
An example of CROSSTABS output is shown in Figure 2-3. 

The creation of the SPSS files FORSTUDY and WEATHER has 
provided the Alberta Forest Service with the capability of analysing 
many variables pertaining to forest fires with relatively little effort. 


For example, a simple program consisting of less than twenty-five IBM 


am i eis 


&S ofdel ab awode ei watdntaov or to Ef A seca, 2890 


7 
| 3 

| wieben ab 47 Coke Theva eat Yo aeusaihh Bea 2stod Sam = 
bongiash ek bas aosrgeg : scl aie wok Loe ter Wopen oF camera - 


sasyxoqak m6 ef akay 
bangiees sas esti? sa010%" ak 


gsulev Ishtixes bab tambo oasbonosos ox 

aaldetvav old 30 Baal esate Josqse | 

bowolle 8292 snes A icdscbr iar ud bosivogstea) to asusev Isntbx0 | 
jokuteth asgasx yo sesz0t byead canal oe to : seceaeanntne: six 


i) ab 
od Ylevienbtxo bseu ew (2HATe sdsosbnsqabat a 
woneups fi nisido of ~ a sixe) oe asouied aes mee one r. 
wag ; as i a i on Asa 

basedatctnd’ od Mien ? 


a 


pei Bas} pines eae phish ps5 vat 
eer a oan 


baped eee om . veo 
mye Ale ies Sad A , 


a1 


TABLE 2-2 


A LISTING OF THE VARIABLES IN THE SPSS FILE WEATHER 


VARIABLE DEFINITION 

1. Station Weattier station reporting the data 
2. Month Month in which data reported 
Sen CAar. Year in which data reported 
4, Day Day the data is reported 

5. Time Morning or afternoon report 
G2 Sky Condition of the sky overhead 
7. Weather 1 Rain or types of rain 

8. Weather 2 Thunder, hail, smoke or snow 
9. Weather 3 Factors affecting visibility 

10. Visible Visibility in miles 

1; tevemp 1 The wet bulb temperature 

12. Temp 2 The dry bulb temperature 

ls Pret se 3 Relative humidity 

14. Windir Wind direction 

15. Vindsped Wind velocity 

16. Gusts Velocity of gusting wind 

te (Gloud of High clouds 

a8. Cloud, 2 Middle clouds 

19. Cloud 3 Low clouds 

20. Cloud 4 Cumulus clouds 

21... Cloud °5 Heavy cumulus clouds 

22% ‘Cloud 6 Cumulonimbus clouds 

VLE cae Condition of the grass 

24. Buildidx The Build-up index 

25. Spreadx The spread index 


26. Remarks Any general remarks 
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cards will cross tabulate the time of day of fire start with the 
cause of the Bios R.S. Miyagawa has used both files in his capacity 
as Sere ee 

The sequence leading up to the creation of the SPSS files 


is shown in Figure 2-4, 


Provincial ftre 
data stored on 
IBM 7070 tape. 


Provincial 
weather data 

stored on IBM 
7070 tape. 


Re-organized 
and examined 
data. 


Re-organized 
and examined 
data. 


Transfered 
data to IBM 
360 tape. 


Transfered 
data to IBM 
360 tape. 


Created file 
FORSTUDY. 


Separated and 
transfered Footner 
weather to IBM 360 
tape. 


Created file 
WEATHER 


FIGURE 2-4 


SPSS FILE CREATION 
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Relevance and Potential of the Approach 


The simulation model described in this thesis can be used 
as a research base for studies on suppression forces, presuppression, 
and fire prevention. For example, a study of forest fire suppression 
would begin with a working knowledge of what sizes of fires the sup- 
pression crews are fighting, where the fires are located, what time 
of day the fires are started, and the rate of spread. The described 
simulation model calculates the growth rate of every fire, the location 
of each fire, the time of day of fire start, and the fuel type feeding 
the fire. Equipped with this knowledge, a designer of suppression 
forces could expand the model so it incorporates the locations of 
alternative suppression headquarters. The elapsed time from discovery 
to suppression crew arrival could be calculated and used as a measure 
of the effectiveness of the alternative suppression force locations. 

A study on fire prevention could also use this model as a 
base. The model calculates and outputs the cause of each fire and, 
as further explained in Chapter 4, the growth rate of a fire is a 
function of the time of day of fire start. A modification of the 
fire cause probability distribution will affect the growth rate 
of the fire which directly affects the total cost of the fire. 
Since the model indirectly describes cost as a function of cause, 
a researcher could measure the effect of reductions in any of the 
10 cause categories recognized by the model. The reduction in cost 
becomes a measure of the funds to be allocated to achieve the cause 


reduction objective. 
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The study of forest fire detection in the Footner Forest 
is by no means complete with the operation of this model. The model 
Opens doors to many other research areas because the results of other 
research can be tested in the model. One such research area could 
be aircraft scheduling. While operating the model as described in 
Chapter 1, it became clear that the number of fire discoveries by 
aircraft in the Footner Forest is very sensitive to the time of 
day the aircraft are flown. For example, the percentage of aircraft 
discoveries is doubled if the aircraft are flown at 1 p.m. instead 
of 10 a.m. What is suggested is the need for an optimum scheduling 
routine which considers the site of fire at discovery. 

Another area of research suggested by operating the model 
lies in the area of the public detection system. As illustrated in 
Chapter 1, the detection system in a forest is very sensitive to the 
informal public detector and research in this area could prove 
extremely rewarding. 

Information required for studying a specific problem is 
unique. The output from this model is cost orientated since the 
research measured detection effectiveness as a function of total 
cost of detection plus suppression. This does not limit the model 
to cost analysis only. The print out can easily be modified to 


include many other parameters of interest to the researcher. 
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CHAPTER 3 


STATE OF THE ART 
AESASbAS i Decne ioe Soatens 

Research on Forest Fire Detection Systems has taken 
four forms over the years. There have been field trials to test 
and compare alternative detection systems; there have been computer 
based simulations and analytical models to analyse cost-effectiveness 
and describe fire growth; there have been national conferences to 
describe fire problems and set out research guidelines; and there 
have been research studies on new technological areas such as infrared 
detection. 

The field trial studies have usually been conducted during 
the fire season to compare aerial detection and fixed tower detection 
alternatives [4] [14]. Field trials have also been conducted to 
test the effectiveness of infrared scanners [1] [6] [7]. A field 
trial is an experiment and as such can be classified as a form of 
simulation. A typical field trial could leave some of the towers 
unmanned and replace them with aircraft detection. The substitution 
of aerial detection for fixed tower detection could then be analysed 
over the course of the season and a decision made on the effectiveness 
of the new system. 

Field trials are very useful tools for testing some aspects 


of forest fire detection such as visibility or blind areas. The use 
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of smoke generators has been quite effective in this regard. Smoke 

generators are reported to have been used to check on the efficiency 
of the detection system [4] [14]. In one such report, the forestry 

official said the use of smoke generators had three objectives: 

a) To check on the detection Bt 

b) To provide data on visibility conditions. 

c) To keep up the interest of the observers on routine 
flights and familiarize new observers with the grid 
system by locating a smoke and pinpointing it on the 
map [14 pp tol. 

When analysing complete detection systems however, the field trial 

is too slow and perhaps very costly. It could take several years 

to accumulate enough information and test enough design alternatives 
to enable the designer to select the best design. Indeed, the data 
may be erroneous if the independent variables associated with fires 
cannot be controlled or accounted for. Factors such as an extremely 
hazardous season or very wet seasons effect the average results if 

the avensee is not calculated over many years. Field trials are 

very difficult to control. With the invention of the modern computer, 
the designer has been able to simulate many fire seasons in a matter 
of seconds. The computer based simulation of forest fires and the 
accompanying detection systems allows that several design alternatives 
be compared while the independent variables such as cause, time of 


day, fuel, etc. all being accounted for. The most impressive simulations 
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of forest fire detection systems had previously been designed by 
Pato Kourta: [6] [10]. The first simulation Kourtz did was for a hypo- 
thetical forest consisting of 2,000 square miles divided into three 
Pecoree The design variables were the number of lookouts, number of 
aircraft flights over areas not protected by lookouts in each sector 
on each danger class day, and the time of each air patrol. The environ- 
mental variables were daily danger index, daily visibility, number of 
fires in each fire occurrence sector and the time of occurrence of each 
fire. The objective of the Kourtz simulation was to find the alternative 
detection system that maximized effectiveness at a given budget level. 
The measure of effectiveness was average area burned per fire up to 
the time of detection. The Kourtz study was theoretical and used 
hypothetical data. Unlike this study, it did not take into consi- 
deration diurnal effects of fires burning slowly at night and rapidly 
during the heat of the day. The Kourtz simulation assumed an average 
fuel and did not distinguish between fire causes. Table 3-1 illus- 
trates the comparative differences between the Kourtz study and this 
study. In spite of the aforementioned weaknesses, Kourtz displayed 
a valid technique in evaluating detection systems. His work 
provided a basis for this thesis. 

A simulation program to evaluate airborne infrared detection 
was presented in a paper by Kourtz in Missoula, Montana in June, 
1971 [10]. This model incorporated a Monte Carlo technique for 


determining environmental data to operationalize the simulation. 
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TABLE 3-1 


COMPARATIVE CHARACTERISTICS OF THIS STUDY AND THE KOURTZ STUDY 


RELEVANT CHARACTERISTICS OF THE 


KOURTZ STUDY 


Criteria chosen as the Average 
number of acres burned up to 
time of detection. 


Average fuel assumed in 
growth model. 


The growth model used was deve- 
loped to study suppression 
forces. The Parks Model was 
driven by the danger index 
only [15]. 


The forest was hypothetical. 


The air routes were designed 
probabilistically. 


Towers were hypothetical and 
did not have blind areas 
associated with them. Discov- 
eries were probabilistic. 


The public could discover every 
fire based on a poisson distri- 
bution. 


The general level of sophi- 
stication was low. 


CHARACTERISTICS OF THIS STUDY 


Criteria chosen as the mini- 
mum total cost where Total 
cost = Detection cost + 
Suppression cost + Damage 
cost>. 


Four fuels were considered 
in this study. 


A growth model was developed 
for the study area and was 
driven by time, fuel and 
burning index. 


The forest was real. 


Route flying in simulation of 
actual practice and ability 
to spot a fire was related to 
the daily visibility. The 
route locations were deter- 
ministic. 


The towers were real and had 
actual blind areas associated 
each of them. Discoveries 
were deterministic. 


The public found only those 
fires not detected by aircraft 
or towers. The elapsed time 
was calculated using an empi- 
rical distribution. 


This model used over 400 proba- 
bility distributions to simu- 
late the real world. 


1. Damage cost was not included because no reliable data exists. 
However, the simulation is designed to handle such costs. At 
this time it is possible to use the model to do sensitivity 
analysis with respect to ascertainment of the range of damage 
cost that would not change the tentative conclusions reached 


in the economic analysis. 


This is an avenue for future research. 
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Probability distributions were determined from historical data on 
thunderstorm occurrence time, spread dies, fire ipeatdout ignition 
time, hold over time and cloud cover. The growth model used in the 
simulation was developed by Kourtz and O'Regan and will be discussed 
in the next section. Also presented by Kourtz at this seminar were 
examples of existing computer programs to calculate visible area 
from contour maps and store fire data such as that accumulated on 
fire report forms. An outline of the later program was given in 
a December, 1971 publication [9]. 

A study held at Woods Hole, Massachussetts, July 17 to 
August 11, 1961 under the guidance of the committee on Fire Research 
of the National Academy of Sciences - National Research Council 
resulted in a publication of their recommendations and conclusions 
on fire problems in the United States [2]. Appendix VIII of the 
report dealt with operations research in forest fire problems. 
Appendix VIII states: "there is merit in undertaking this report 
as an initial effort to improve the understanding of forest fire 
managers on the potentialities of operations research studies and 
to assist in the recognition of problem areas by operations research 
specialists", The study committee obviously recognized the usefulness 
of operations research in studying forest fires. The committee also 
recognized the need for data to "develop and validate a suitable fire 
spread model" [2, p. 163]. Professors Shephard and Jewell felt that: 


"in all cases of data gathering for operations research studies, it 
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is important to have collections designed to serve precisely 
the model structure involved" (2, py 1.63)% 

The Woods Hole conference marked the beginning of a new 
era in forestry problem solving. In 1961, the U.S. Department of 
Agriculture Forest Service entered into a cooperative agreement with 
the University of California to develop a long range operations 
research program in the field of forest fire control systems [3] 

[5] [13]. The studies done under this program were in the form of 
mathematical models designed to provide optimal solutions to 
suppression methods or establish fire growth in a mathematical 
form that could be used in suppression or detection optimizing 
routines. Simulation, as a technique in analysing forest fire 
detection was not used under this program. 

The research in infrared detection systems has been headed 
by S.N. Hirsch, F.H. Madden and R.L. Bjornsen [1] [6] [7]. The 
study falls under the general heading of Project Fire Scan; a project 
to determine the capability of infrared systems as a forest fire 
detector and to determine the application of infrared systems for 
delineating the perimeter of large fires when smoke obscuration 
prevents obtaining the intelligence by conventional means [6]. 

The necessary infrared equipment had been developed by military 
contract and it was Project Fire Scan's primary objective to 
determine the probability of detecting targets in timber cover as 


a function of vertical angle, timber characteristics and fire size. 
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Field trials have been conducted over various timber canopies and 
various fire eee The preliminary results published by Hirch 
indicate the infrared detection system could be an economically 
feasible system of forest fire peter rion, Infrared detection was 
not considered in this study. This study was limited to analysing 
the existing detection system and infrared was not a part of that 
System. Future research could incorporate infrared once its 


effectiveness and limitations in Alberta's forests are known. 


Fire Growth Models 

Several researchers have had a hand in the development 
of fire growth models. A few of the prominent models will be 
discussed and compared in this section. 

A fire growth model developed by G.M. Parks was designed 
to predict fire growth during the time interval from attack on the 
fire until its control [15] [16]. The model took the form of 
a = G + HT where: A is the area burned in acres. 

T is the elapsed time in hours. 

G is the linear growth rate in acres per hour. 

H is the acceleration in acres per hour squared. 
Integrating this equation with respect to time: 

A = GI + 1/2 HT2 
The values of G and H are estimated by fitting curves to actual fire 


histories for any particular area. Presumably this could be done for 
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different fire danger indices and different fuels although no 
literature was found where the model incorporated different fuel 
types. Although the Parks model was not specifically designed to 
model a fire from start until detection, it has been used for 
this irnoden 

A fire growth model was developed by PH. Kourtz and 
W.G. O'Regan using a Monte Carlo simulation technique [12]. The 
model described a smoldering or creeping ground fire that had an 
area of less than 0.05 acres. The model was capable of simulating 
for a given fuel and weather conditions the expected perimeter 
location and burning area of a small fire at any time after ignition. 
A computer program was written to combine data on fuel grid, variable 
spread rates, persistence times and the optimizing procedure. The 
fuel grid data was calculated by a Monte Carlo sampling technique 
and a set of hypothetical fuel distributions. The variable spread 
rates were determined by fuel type and spread index class for the 
day. The persistence times were assigned to each ground fuel type 
for each spread index. The persistence time was defined as the 
time in hours that a fuel type could sustain fire after the initial 
spreading fire had burned midway through the grid square. The opti- 
mizing procedure incorporated the Dijkstra algorithm to find the 
shortest time path between each square in the fuel grid. The model 
was not verified because of lack of data, however, it was incorporated 
in an infrared detection study design to compare the effectiveness 


of an infrared system versus a visual system [10]. 
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A geometric model of fire growth was developed by C.E. Van 

Wagner [17]. This model assumed a burning fire approximated an 
elliptical shape where: A - area of fire. 

a - long semi axis of ellipse 

b - short semi axis of ellipse 

V - linear rate of spread at the head 

u - linear rate of spread at the flanks 

w - linear rate of spread at the rear 


t - time since ignition 


(1) A =Tab 
2): aven +ow)t. / 25 
(3) _bp=put 


Therefore the area 
(4) A= _IL (v + w)ut2 
2 
The growth rate of the fire is then expressed in the form 
of a velocity component and an acceleration component. The velocity 
Aversa. ea oe 
component is: aoe T (Vv-aw)ue 
and the acceleration component is: d7A =T (v + w)u 
dT2 
The model can be used only if values of v,w, and u are known for a 
given fuel type, danger index and specific burning conditions such 
as fuel moisture and wind. As is the case with most geometric models, 
the model depends on a uniform fuel. 
A summary of the characteristics of the growth models 
previously discussed as well as the growth model presented in Chapter 4 


is shown in Table 3-2. 
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CHAPTER 4 
THE FIRE GROWTH MODEL 


It ronieidon 

In 1951, a report on the analysis of 36,000 forest fires 
in the Northern Rocky Mountains was published [1]. Barrows’ work 
was aimed at gaining a better understanding of occurrence, behavior, 
control, and effects of fire in order to try to answer a fundamental 
question: What are the requirements for adequate fire control at 
least cost [l, p. 137]? Since 1951, there has been a tremendous 
amount of work done in the analysis of prevention-detection-sup- 
pression-reforestation systems. However, even with technological 
developments such as the infrared forest fire detector the economics 
of the current state of technology is such that the design of fire 
detection systems revolves around the effective use of three major 
system components; which are manned towers, surveillance aircraft 
manned with trained observers, and the general public. A simulation 
model has been developed to aid in least-cost design of forest fire 
detection systems within the Province of Alberta. The development 
of the design model required the use of a fire growth model which 
provided information so that management of the Alberta Department 
of Lands and Forests could use the results provided by the simulation 
model with confidence. The purpose of this Chapter is to outline 
the approach taken in constructing a computer-based fire-growth 


model and the results obtained using this model. 
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Fire growth models of potential usefulness for cost 
considerations in systems analysis fall into four major categories. 
These categories include theoretical-physical, geometric, deter- 
ministic-analytical, and probabilistic models derived specifically 
for use in computer Sealers. The theoretical-physical models 
require measurement of physical parameters which exceed the limits 
of currently existing field data and, therefore, offer a degree of 
rigor not matched by existing data [3]. Geometric models which 
assume fire spread according to some well-defined growth law, and 
take a standard shape such as an ellipse represent the other analy- 
tical extreme and are usually not sophisticated enough to provide 
the precision required for use in a complex systems analysis [8]. 

The deterministic-analytical model developed by Parks focuses pri- 
marily on fire growth during the time interval between initial 

fire detection and final fire control [7]. What was needed for 

this study was a fire-growth model which provided a realistic 

growth in size up to the time of fire detection. The model presented 
by Kourtz and O'Regan was developed to portray creeping ground fires 
of a size less than 0.05 acre [4]. The simulation outlined in Chapter 
1 required a model applicable to a range of fire sizes extending from 


less than .25 acre to greater than 500 acres. 


Development of the Fire-Growth Model 


Given a fuel type, burning index, and the time of day of 


fire start the growth model computes the size of the fire at the 
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time of discovery. This information is then used in a Monte Carlo 
fashion to compute the cost of fire suppression. The development 
of the growth model is described in detail below. The idea for this 
approach was provided by Miyagawa and Stashko [5]. 

The first task was to develop rate-of-spread classifica- 
‘tions for Footner Forest fuel types. An available map of fuel types 
for Footner was broken into the four categories of coniferous, 
mixed wood, lichen-muskeg, and grass. Table 4-1, provided a guide- 
line for a rate-of-spread classification for Footner [2, p. 40]. 
In developing the rate-of-spread classification shown on Table 4-2, 
the following two basic assumptions were made: 

(1) That the partitioning of fuel types, as shown on 
Table 4-2, is reasonable. 

(2) That the forested lands of Footner are approxi- 
mately equally distributed with respect to valley 
bottom, benchland, slopes, and ridgetops. 

Next, the fire-danger rating referred to as the Burning 

Index, was developed. We assumed that the Burning Index was a 
function of two other fire-weather indices called the Buildup Index 
and the Spread Index. These two indices are calculated according 
to the instructions provided with the set of tables for use at 
Alberta Forest Service Weather Stations [9]. Daily records of 

the Buildup Index and the Spread Index have been kept, in Alberta, 


during fire seasons for the last nine years. SPSS subroutines were 
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TABLE 4-1 


RATE-OF-SPREAD CLASSIFICATIONS FOR ALASKAN FUEL TYPES! 


Valley bottom Benchland Slopes 


Fuel type Reha mEDIa Sem Ley a. SauthenvaaNorhaly. Ridgetops 
White spruce or 

birch-spruce M H H E E H E 
White birch or 

birch-aspen M H H E 3 H E 
Black spruce H E H E-F E-F H H 
Aspen M M M M H M M-H 
Cottonwood L M L M M M 
Willow-alder M M M ie H M M 
Grass E F IP. F F F F 
Muskeg M H H E F H Ec 
Tundra M H H E E-F H E-F 


'Rate of spread: L=low, M=medium, H=high, E=extreme, F= flash. 


Based on BI of 40: 3 m.p.h. wind; 30 percent relative humidity; severity index 8; today’s slat moisture content 6 percent. 


used to analyse this data [6]. Cumulative probability distributions 
for the Buildup Index and the Spread Index were developed for each 
day of the fire season in the Footner Forest. The average fire 
season in Footner spans a time interval of 153 days, therefore, 306 
probability distributions were used, in Monte Carlo fashion, to model 
the behavior of fire-weather conditions in Footner. Given a daily 
Spread Index and a daily Buildup Index, as generated by the simu- 
lation, the daily Burning Index was determined by a table look-up 
procedure. Data provided in the Buildup and Spread Index Tables 

for use at Alberta Forest Service Weather Stations provided a guide- 


line for establishing the relationships shown in Tables 4-3 and 4-4. 
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BURNING INDEX AS A FUNCTION OF SPREAD INDEX AND BUILDUP 


INDEX FOR THE FOOTNER FOREST DURING THE MONTHS OF 


AUGUST AND SEPTEMBER 


Difficulty Of Control Of Fires Should Increase 
As Buildup Index Increases 


OO 


Forward Rate of Spread of Fires Should Increase 


As Spread Index Increases 


Spread Buildup Index 
Index 0-9 10-25 26-50 51-100 100+ 
0-1 1-20 1-20 1-20 1-20 1-20 
Low Low Low Low Low 
2-9 1-20 1-20 21-40 21-40 21-40 
Low Low Moderate Moderate Moderate 
10-19 20 21-40 21-40 41-60 61-80 
Low Moderate Moderate High Very High 
20+ 1-20 21-40 41-60 61-80 81-100 
Low Moderate High Very High Extreme 
* Values in the Table represent the applicable range 
of Burning Index. In terms of fire danger 1-20 low, 


21-40 = moderate, 41-60 


81-100 = 


extreme. 


high, 61-80 


very high, 


“ie ee eran A 
a we } 
hs a m 
Lai 
hi | | 
ny. 
A 
7 
J | qua rave 1H aeons anante bia OT TOMI wiua 
eal Cr eee ey ae 
. 40 aT MOM sur arn c. ‘Pao \aMHDOOR aut 5 Ya cto a 
Pa RR Ci a?) i ak te 
| Savana We Gs TEHOUA ” em ie 
5 ; ea a " hae » 
ra am ee oT eae 
. . ie ya ees See 
a ,, “ Ae > A > 
szec9t>nl Bliod? eerlt 16 forth) © 20 visa) : we 
aan ntont Xabwl qubliva tA i ah 
~ 


~ A ee — 


z ’ 


4 _'s 
yy 3 ) 
gs 

o 

Baie eeahes dali, | i a). 5 
; jwn 2's Ons Obes" | tket a) os a 8 
eputobol 62 <%>bo4 asuacboM wad wo mit » (age 
” y c ? P 7 : 2 
tl gupnit ee eremneiannngp stare et: Slee St eg rr ah a 
erie , - ; Pd nt re 
04-L28 %-Lh O}~-£5 Oh-£5 » OXeL "SO Ts08 ¢ a 
feilt vray aot ‘statobo SieTebolv woul 
é : at bay 2 
nempe ner meemininnn, 10st ovantiandscli a Ce - 
- ; wt, ai at ; 2) 
1-18 . o8-te 8-18 =—s Hees Ose ss SS > 
geors4 Agth vrov r agit ofetaboM vod 1; pe Uta 
¥ - ¢ ae , 
* i hg ‘ ws, 7 + 
‘ ie) ue § ere 
Ne eps 
at a * 1 Ca bide a 
pg tis oidoot iyqa ody! saantgor ot det ody at 20008¥) i ‘ 


wok = O5-L, togmeb” $212 2o emr03 al .xebrl aqparptt to ae 
| 
J 


ligt ytew = 00-16 we = O8-t. 33 prohom «| CS : 
| ; 1 OAO1I KD, 001-28 Ti 


a : : . ahs 7 P } x a 


S38, 


Given the fuel type and Burning Index provided by simulation, 
the average initial rate of fire spread can be ascertained. The data 
shown in Table 4-5 was used as a starting point for estimating the 
fire perimeter increase applicable to the Footner Forest [2, p. 42]. 
Experiments were run with the design model, continually modifying 
the values shown in Table 4-5, until a set of data was obtained which 
produced what seemed to be realistic average initial rates of spread 
for the Footner Forest. The final set of data which was used to 
drive the growth model in subsequent analysis is shown in Table 4-6. 
The values shown in Table 4-6 have been adjusted for diurnal fluctu- 
ation by using the information for Alaska*Stick as presented on 
Figure 4-1 [2;p.48]. We were, therefore, able to develop a fire 
growth model based on fuel type and burning index which produces a 
different growth rate for each hour in the day. In this sense, 
the growth model is dynamic. When one examines Table 4-6, it is 
evident that minimum values of perimeter increase shown on Table 4-5 
provided lower bounds for values adjusted for diurnal fluctuation. 

In order to use the values obtained for fire perimeter increase to 
compute the size of the fire, in acres, at discovery, the following 
three assumptions were made: 

(1) That the fire burned through a single fuel type 

that was either coniferous, mixed wood, lichen- 
muskeg or grass. 

(2) That grass fires occurred only during the month 


of May in the Footner Forest. 
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(3) That the area of burn was circular in shape and 
could be computed according to the simple equation 


derived below. 
TABLE 4-5 


AVERAGE INITIAL RATE OF SPREADL ACCORDING TO 


FUEL TYPE, SLOPE STEEPNESS, AND BURNING INDEX AT SITE OF FIRE“ 


Fuel rote of Slope Burning index 
spreod type steepness} 1-10 11-20 21-30 31-40 41-50 51-60 61-70 71-80 81-90 91-100 
Percent Perimeter increase (chains per hour) 
Low 0-10 O 1 1 1 J Pe 2 2 3 4 
11-25 i¢) ] 1] 1 2 2 3 3 4 6 
26-50 1 ] 2 ¥ 3 3 4 4 6 8 
51-75 1 2 3 3 4 5 6 Zz 9 13 
Over 75 2 3 4 5 6 7 8 10 15 20 
Medium 0-10 0 ] 1 1 Pe 2 2 3 4 5 
11-25 1 ee 1 1 2 2 3 3 4 6 7 
26-50 1 2 2 3 Ss 4 5 6 8 10 
§1-75 2 3 3 4 5 6 7 9 13 16 
Over 75 3 4 5 6 8 10 12 15 20 25 
High 0-10 0 1 2 3 4 5 6 Sis 9 12 
11-25 ] 1 3 4 6 7 9 11 13 17 
26-50 2 2 4 6 8 10 12 16 18 24 
51-75 3 3 6 9 13 16 19 25 28 38 
Over 75 4 33 10 15 20 25 30 40 45 60 
Extreme 0-10 1 S 4 5 6 8 10 13 16 19 
11-25 1 4 6 7 9 11 14 19 22 27 
26-50 2 6 8 10 12 16 20 26 32 38 
51-75 3 9 13 16 19 25 32 4) 50 60 
Over 75 5 45 20 25 30 40 50 65 80 95 
Flash 0-10 1 5 12 15 19 24 30 37 46 S¥é 
11-25 1 vA 17 21 27 34 42 52 65 81 
26-50 2 10 24 30 38 48 60 74 92 114 
51-75 3 16 38 48 60 76 95 117 146 181 
Over 75 5 245; 60 75 95 120 150 185 230 285 


lAverage initiol rate of spread refers to perimeter increase between discovery of fire ond first attack. This rate of spread may 
be onticipated during the first 4 to 5 hours. 

2Table based upon study of 2,955 fires in Notional Forests, R-1, 1936-44. Values for very high and very low burning index 
have been estimoted. 

3General descriptions used in slope descriptions are: Level, 0-10 percent; Gentle, 11-25 percent; Moderate, 26-50 percent; 
Steep, 51-75 percent; Very steep, over 75 percent. 
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FIGURE 4-1 


BURNING INDEX FLUCTUATION VS. TIME 


! | i 
BURNING INDEX 
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C = Perimeter increase in chains per hours 
A = Size of fire at discovery in acres © 

t = Time in hours 

CE= Zhe 
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A = 4356C2 /41 square feet 


A=c* / 407 acres 
If a fire burned for t hours prior to discovery the total 
area burned is computed according to equation 7. 
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For example refer to Table 4-6. If a fire started at 0500 
hours, in a coniferous fuel, on a day with a burning index 
that reached a maximum in the 21-40 range, and burned 2 
hours prior to discovery the size of the fire at discovery 
would be computed according to equation (7), above, as 
follows: 


A = (8.1)* / 407 + (8.2)2 / 407 = 1.06 Acres 


Results 

The fire growth model is a small portion of a computer 
program developed to aid in the design and evaluation of forest 
fire detection systems. As outlined in Chapter 1, the model was 
validated by comparing results obtained from a computer run, designed 
to model the present operating decision and physical processes 
existing in the Footner Forest, with statistics obtained from analysis 
of nine years of historical data. Table 4-7 shows the results obtained 
with the fire-growth model when it was used to compute the size of 
the fire at discovery for 2,667 fires spread over the 7,956 days 
of 52 fire seasons. The distribution of fires within fuel categories 
is shown in Table 4-8. The visual fit is quite close. When the 
fire sizes are grouped into categories as shown in the lower portion 
of Table 4-7, the results obtained from the fire growth model pass 
the Pee eee ni es test for goodness of fit. This standard 


statistical procedure provides a test of the hypothesis that the 


a0 rz 
086: 48 besz59—, ort? su O-8 afd hei aera 


ey | 
pn gehamd a A3ty, eb gabe ak ad 


ay 


oy. Ce, a. ca, “ 
: S pensud | ine sare ONES its at muse os seit : 
: Byki aw ae ea 7 
gisveseth 36 ost ] ‘ost to. osie’ ie cravozetb Od tokIg wi08 et i 


e& ,svods ,(%) noksaups oct, gatbr00 2° basuqaoo a star z 
7 ne v1 
i we ffot 


- oe oa, ; 4 : 
asibA 20.2 © WOd \ 748.8) + WOR '\ Ser 8) =A i 


‘) ie)” ni 


a ' . 
‘ j 
' . 25 » | 
. 1) ‘ 
Ss “ ; oF f ube 
aie = : 
< : 


owe \ th 


sostuqmoo 8 fo noidxoq Clams ist ddwo tg, ont oft 


sn? 


teoro1 To motsaulavs bos agtasb Sa ak bre of - beqote: rob nts 
esw [sbom ad3 1 reaqe mt benkisuo ail -2medeNs notsaaseb ¢ yt 
a 


bengiesb .aut isj3uqmos & dis boateido « ah uaey sate. wd bess sbi v 
: re 

B9B8e990%q fsoteyde, bas potetseb ; galss aI ah ct orld Labor mod 

i = ‘ 7 

ete, Kioos mort Doalsido eotsatiare fitiw .426 ft aria oF satan 
ev 


® Pa 4 


. > 
“bent sado le ak sit ewode t-a oldet s60Rb inabaogebt "20, axtey ome 


“ a 
jo nel oni3 i o3 ‘eee’ 2au at nadw Lebom ‘iii sgrontt ~~ 
o 


—eyAb ace, t ate rove oe oe, (28, £ vot seplbeins: is 8 


rxogess2 Loud uae asst ai jase roath nit sore sath 


} of i hibits oat 2 


ake ‘al : 


#02 Howl: as ie he 
mol3 104. ‘xawol at nah anon ‘aaa baq query 
[Pr 7 rar 
eesq tonne ont Yad 3 ee bie edi ie 4h 


ai Da 
Pa er ee 
biebasde ehit $ ah “ 280b 


= 
= 


a 7 


<5 ' Fi eee 
8, ae 
ai ats oodsoa at We 
i € sir ion - i 4). 
ye et 
Jj : 


Sinoy 71 Ajai wmyxoadde een pazrdde 210A BanteA osay. YOTYR 02 WINq 203 IFWTT eWFa YO puNog saddn oyz tT 


Ne eee a a ES 


0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°G6 0°66 0°66 0°66 0°66 0°66 OOT-T8 
0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 08-19 


0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°66 0°GG 0°GG 0°66 0°66 0°66 0°G6 0°66 0°GG 0°66 0°G6 0°66 0°66 0°66 0°66 09-Tr (yseta) 
1°06 0°26 S*26 0°S6 S°b6 0°96 $°L6 0°66 0°66 0°66 0°66 0°86 0°46 0°96 0°S6 0'¥6 0°S6 S'76 8°16 Z°T6 8°06 ¥'06 0°06 0°06 = OF - 12 ssei9 
er a oOnL Orie taat ra Ooy (07m, OP 8S Se 282 OL CONT (Oot Ort Ot Ot Ok Olt ON orl OFT 02-1 
Se a a IS a 
5°15 0°6S 0°29 0°S9 0°F9 0°S9 0°99 0°L9 0°49 0°L9 0°99 0°S9 O'¥9 O°F9 0°29 O'T9 0°09 S*6S O°6S SBS O°BS S*LS OLS OLS HOT -TS 

S°Sp OLE 0°0S O'TS O°ZS O°SS O°PS O°SS O°SS 0°SS OPS OSS O°2S S*TS 0°05 O'6b O'8h S*Lp OLY S°9H O'9Y S*Sh O'S O'SP 08-19 (eesae 
ee OOS O'ZP O'SP O'bD O'SY O°9P O'Ly O'Lb O'LY O'9D O'S O'by O'S O'Zy Ol OOF S°GE O'GS S*8E OSE S*LE OLE O'ZE 09-1t a 
a Glee t Gee Gg OTR PRM G 8 69) OG Oe) OG 672 878 2°78 98 S°@ #°8 $°8 iz's (2:8 TJs Ts) Ors O'8 “b°8 Ob-T2 sno1ajtuo} 
POL OC Ear Tar Ot Obie 2) SS EE G8 272 Ol OE. Vk OL OL OOK OT O'R OO Ot) Oy 04 02-1 

SUE OAS SSS Rae eI Be Sele SS IRS ASA Sa SO at ec eae skae aneie nn Laan aete tk eet arene tase Dime: Bierter aie Tyee ie st 
OI O'ZE H'ZE Z'EE BEE Z'VS O'SE O'SE O'S Z°PS OVS BES HSE OSE O'ZE H'ZE Z°ZE O'ZE OTE HTS OTE 9°OS OOS O'OE OOT-T8 (y3ty) 
PLU B°LT O'ST Pi8T 8°8T 2°02 9°02 O°TZ O'IZ O'IZ 9°02 Z'O2 BGT v'GT O'6T O°ST Z°SI O'RT BLT O°LT HLT Z°LT O'LT OL 08-19 ee 
STL E 2LSe EL SFT SSE SST S291 SOLES SO StOl SOT 8° St CSE SPT Syl SST SST O Sa ZT pact CUI SLi Seite sett 09-Tr Stee 
Gica SOD. 1, 2° 9mC 0679, SHO MSE GO) PO. Grd ZO 179 O°b 6S 8°S Z°SeaL"S Ops 9°S™ 67S: SS) 57S Or-12 WPT 
TOR ORT CORIO TIO 1) O°k nO tmmOat) OCT 08k 0° 0°, 60°R OOP O° Ot OT (Oct OM OT SOP oct 02-1 

0°61 0°02 9°02 O°T2 b°TZ O°ZZ 9°22 O72 O"EZ O°EZ O°7Z H°ZZ O'ZZ Z°IZ-8'OZ. ¥°OZ O°OZ O°6T H°GT OST O°SE B°ST O'ST ORT  OOT-TE 

pik Peo sl SLY ST. OST SCL UCPIMOGON OPT O°ST OST PET Z°ST OTST GCI 8°21. 9°21 O21 Cee Coe OlZ1 Oct Or zt 08-19 canta 
boo b 96) 0°G) 2-6) t° Ga OG 86 OROTROUOT O°0T 8°6 976 $°6 “226 0°6 8°S 9°8- 9° pS zee 2°e) o's Of O78 09-Tr Pan) 
Gouin e = (Occ) CTO Zee Gn Hoc! Ce MEGS KG? oh C $5'g ZIG ERE OS! GSe SR NZ. 22 ORE 972 S82. SIZ S*2 Ob-12 PpooM PpaxTx 
Dole Get MO ba O2t 0cT, OF ROT AO 0 Ost: «OC OE OTT OT. Of OCR OCT O°T OTT OMT OT. 025 .05T. 0°T 02-1 

DOr? O0SZ 0022 OOTZ 0002 O06T OO8T OOLT 009T 00ST OOPT OOSE OOZT OOTT 000T 0060 0080 0020 0090 0050 00%0 00¢0 0020 O0TO 


ae Salah eq ET) 41noH xopu] addy peaids 


InOH 2Joq SuyeYy): Ul- eseoiduy Ia.eWti9sg Buyuing jo a2ey Tang 


yisddod MANLOOA YOL AVA AO UNOH ANY ‘SXHGNI ONINUNG 


‘AdAL THNd OL ONIGXOOOV AVAUdS AO ALVU TVILINI 3DVAAAV 


9-7 HTAVL 


y 


i 
t 
t 
7 
‘ 
E 
i 
i 


ry 
ya) 


60 


TABLE 4-7 


FIRE SIZES AT DISCOVERY AS COMPUTED BY THE FIRE GROWTH MODEL 


#In Category Percent of Occurrences Percent of Occurrences 
Size Category Produced. by Expected on Basis of 9. Computed. by Fire Growth 
Acres Simulation ‘Year Historical Analysis ‘Model Sample Size 2,667 

Spot 1,047 44.2 39.3 

aap = 2.45 601 18.3 22D 
1.46 - 2.45 152 4.6 5.7 
2.46 - 3.45 126 4.4 4.7 
3.46 - 4.45 64 1.6 2.4 
4.46 - 5.45 40 4.7 LSS 
5.46 - 6.45 48 Lez 1.8 
6.46 - 7.45 48 0.6 1.8 
7.46 - 8.45 9 1.4 0.3 
8.46 - 9.45 24 0.8 0.9 
9.46 -10.45 11 5.0 0.4 
10.46 -25.00 164 5.0 * 652 
25.01-100.00 186 4.6 7.0 
100.01-500.00 136 3.2 3.1 
500+ ahi 0.4 0.4 
Totals - 2,667 100.02 100.0% 

Regrouped Frequencies 

Spot — “1.45 62.5 61.8 
1.46 - 5.45 TS s3 14.3 
5.46 - 9.45 4.0 4.8 
9.46 -100.00 14.6 3.0 
100.01 -500.00 Shey Died 
500+ 0.4 0.4 


Total 100.0% Total 100.0% 
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set of data generated by the model has the same frequency distri- 
bution as the set of historical data. 

Empirical probability distributions relating fire size 
at discovery to the final cost of fire suppression were obtained 
through analysis of the nine years of available data. Given the 
fire size at discovery as computed by the growth model, it was 
possible to compute, using a Monte Carlo procedure, the cost of 
fire suppression. The cumulative probability distributions used 
in computation of fire suppression cost are shown in Table 4-9. 
The cost results obtained, for 2,667 fires spread over 7,956 days 
of 52 fire seasons, following the procedure described above are 
shown in Table 4-10. The results obtained pass both the Chi-Square 


and Kolmogorov-Smirnov tests for goodness of fit. 


Conclusions 

A fire-growth model has been developed which provides 
accurate information, for cost analysis, when used in conjunction 
with a simulation model constructed to aid in the design and 
evaluation of forest fire detection systems. Although the fire- 
gpowth model was based on operating data applicable to the Footner 
Forest in Alberta, the approach taken in development is entirely 
general and can be applied to any forest where both fire-weather 
data, as obtained from field measurements, and fire data of a general 


nature are available. 
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TABLE 4-9 


CUMULATIVE - EMPIRICAL PROBABILITY DISTRIBUTIONS RELATING FIRE SIZE 


AT DISCOVERY TO THE COST OF FIRE SUPPRESSION - FOOTNER FOREST 


Cumulative Probabilities 


Fire Average Costs Used, $ 
Size 
Acres 25 48 268 742 2059 «6956 14650 34094 60666 87857 78529! 
Spot 203 428 698 820 987 1.00 
og eihas Bice: 259 478 630 967 .989 1.00 
Wr 087 130 .260 £303 825 868 955.8 £7400 
ate P27 72 59 aes .908 1.00 
346. 4eS 001 1125 250 251 875 1.00 
4.46- 5.45 042 043 .168 6251 334 .959 1.00 
S462 wonas eS 334 500 501 833 1.00 
enac- 74s 5333 334 isae 667 1.00 
7.46- 8.45 001 902 143 429 is 1.00 
8.46- 9.45 001 .250 500 501 1.00 
9.46- 10.45 040 .980 200 360 760 .880 881 £0200) 921 960 1,06 
10.46- 23.00 040 160 400 401 720 $80 .960 1.00 
25.61-10C.090 001 .087 174 348 739 740 .870 957) 7958 959 1.06 
100.01-500.09 001 002 063 1138 438 565 27351 .939 940 941 1.0¢ 
500+ 001 002 003 004 005 006 007 .008 .009 .500 L0G 
TABLE 4-10 


COMPUTED SUPPRESSION COSTS USING OUTPUT FROM GROWTH MODEL 


AND PROBABILITY DISTRIBUTIONS SHOWN IN TABLE 4-9 


Category Number In Percent Of Occur- Percent of Occur- 
Fire Cost Average Category rence Expected On rence Computed By 
Range Used Produced By Basis Of 9 Year Simulation Sample 
$ $ Simulation Historical Analysis Size 2667 
No Data? 25 402 APS gukes 15.0 
1-100 48 358 i325 13.4 
101-500 268 543, 20.9 20.4 
$01-1000 ‘ 742 296 12.0 ye | 
1001-5000 2059 790 29.3 29.0 
5001-10000 6956 114 4.8 4.3 
.10001- 25000 14650 88 2.4 So 
25001-50000 34094 52 a On: : 2.0 
50001-75000 60666 1 0.0 0.0 
75001-100000 87857 5 0.4 0.2 
1000000+ 285294 18 0.8 0.7 
Totals 2667 100.0% 100.0% 


1 ‘small Fires That Resulted In Negligible Suppression Cost 
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Experiments with the growth model, to date, have indi- 
cated that the two fire-weather indicators, i.e., the Buildup Index 
and the Spread Index, are not sensitive enough to the probable high 
rate of spread of grass fires during that portion of the year when 
this vegetation is cured and easily Yom Further research is 
required in order to improve our understanding of the spread of 


grass fires. 
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APPENDIX A 
A MODIFICATION MANUAL 


The simulation model is programmed in Fortran IV G. 
The program was run on the IBM 360-67 computer at the University of 
Alberta computing center, Edmonton. The general operation of the 


model is shown in Figure 1-1, on page 6. 


Data Organization 


A listing of the required data is shown in Table A-l. 
The data deck is set up with the variables ordered in the same 
sequence they appear in the main program (see the main program 
listing, Appendix B). Figure A-1l shows the data deck set up. A 
description of each of the variables used in the data deck set up 
can be found in the dictionary in Appendix C. Each of the variables 
in the data deck list, Table A-l, has been categorized as deterministic, 
stochastic or design. The deterministic variables are calculated 
and read into the program to make it general in its operation. These 
variables allow the designer to more easily change the model to apply 
it to another forest. The stochastic variables are determined by 
the history of fires and weather in the forest under consideration. 
They are cumulative probability distributions usually discrete, having 
a unique value for each category or range of the variable. The design 


variables are those that represent the physical characteristics of the 
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SET UP NAME 
NREPET 


NAIR 


CAU 
COSTT 

SPRING, SUMMER 

Aud i. ALT 2 oe, ALT 10 
FORES 1,....FORES 5 
PNFIRE 

MATT 

IC 

IB 

NROUTE 

NTOWER 

ISEX 

NDAYS 

AMONIQ 

PCOST 

ROUTE 

PUBLIC 


IDC, TOWER1 - TOWERS, 
ITC, IXCORD, IYCORD 


RING 
PFDIST 


MONTH, MCODE, IDAY, 
PS1 - PS4, PBl - PBS 


WHERE 
PCDIST 
SPOT 

PT IME 
BTAB1 
BTAB2 
VIS 

SPREAD 


PFIRE 


TABLE A-1 


DATA DECK LIST 


VARIABLE TYPE 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Deterministic 
Stochastic 
Design 
Stochastic 


Design 


Design, Stochastic 


Stochastic 


Deterministic, 
Stochastic 


Stochastic 
Stochastic 
Deterministic 
Stochastic 
Deterministic 
Deterministic 
Stochastic 
Deterministic 


Stochastic 


INPUT FORMAT 

Tl 

12 
F5.4 

19 
11F6.0 
F10.0 
212 

10A4 

5A4 
F10.5 

1 

12 

12 
I2 
12 
14 

13 
11F7.0 
11F4.3 
See main program 
16F5.4/ 9F5.4 


12, 5A4, 12, 
13.13 


F2.0, F1.0, F2.0, 
4F5.4 


PR Sen eI 
9F4.3 


See main program 
10F5.4 
6F8.5 
10F5.4/10F5.4/10F5.4 
511 
Spat 
6F3.2 
24F3.1 


See main program 
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NUMBER OF CARDS REQUIRED 
it 


1 


1 
1 
One for each size range 
12 x NROUTE 
8 


NTOWER 


F4.3 20 x NTOWER 


No. of months x No. of districts 


NDAYS 
No. of districts 
No. of months x No. of districts 
No. of subdivisions 
3 x No. of cause categories 
4 
4 
4 
No. of months x No. of fuel categories 
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FORES | eeeee FORESS 


SPRING, SUMMER 


DATA DECK SET-UP 


FIGURE A-1 
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detection system. Examples of the values each of the variables 
takes on is shown in Appendix D. 

An explanation of how the data is accumulated and key 
punched on IBM cards pale Before reading the procedure for 
setting up the data for a variable, the reader should acquaint 
himself with the definition of the variable found in Appendix C, 
examine the input format found in either Table A-1 or the main 
program listing in Appendix B, and examine the data listing in 


Appendix D. 


Deterministic Variables 

The values of the deterministic variables are usually 
derived from two sources. One, they are heuristically evaluated by 
examining the system or two, they can be calculated from other vari- 
ables evident within the system. Most of the deterministic variables 
in this model are of the heuristic nature and an examination of their 
definitions will explain the values they take on. However, some of 
the remaining deterministic variables are complex and require further 
explanation. 

The variables CAU and AMONIQ are identical. The values 
for these vectors were calculated by averaging the suppression costs 
for all the fires in each of the cost range categories. For example, 
the average suppression cost of the fires falling within the $501 to 


$1,000 range was $742 for the Footner Forest. 


i nN iti: BA ita 
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The values for the matrix SPOT are determined by dividing 
the forest map up into ranger districts then dividing the ranger dis- 
tricts up into rectangular areas. This makes it possible to allow for 
irregular shaped ranger Tatras. The data for each rectangle is 
recorded on one Sarat This portion of the program is unique for each 
forest and will have to be modified for forests other than Footner. 
Figure A-2 shows how the Footner Forest was subdivided into rectangles. 
The rectangular areas are ordered according to ranger districts. The 
first four cards in this instance are for the 4 rectangles composing of 
DF1, the next card is for the rectangle describing DF3, and so on until 
all the rectangles are described. The data card records the values for 
six parameters describing the rectangle. The first eight columns display 
the reciprocal of the length of the Y side of the rectangle. Columns 
9 - 16 are reserved for the reciprocal of the X length of the rectangle. 
Columns 17 - 24 and 25 - 32 are reserved for the length of the Y side 
and X side of the rectangle in miles respectively. Columns 33 - 40 
and 41 - 48 are reserved for the Y and X coordinates of the southwest 
corner point of the rectangle. 

The tables BTAB1 and BTAB2 are taken directly from the Alberta 
Forestry Brochure relating spread and build up indices to fire hazard. 
The numbers 1 to 5 replace the low, moderate, high, very high, and 
extreme hazard categories. BTAB1 is a spring and summer table and 
BTAB2 is a fall table. Examples of these two tables are shown in 


Figures 4-2 and 4-3. 
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The SPREAD matrix is described in Chapter 4 since it forms 
the basis of the ipeweliblel. The SPREAD table is shown in Figure 4-6. 
Stochastic Variables 

The stochastic variables are probabilistic in nature and 
determined by examining the frequency distributions found in historical 
fire data. The random variable is divided into several categories 
each taking on a different range of values. The number of times the 
value of the variable falls into each of the categories (the frequency 
distribution) is summed and a probability mass function determined. 
The probability of each range is calculated in the usual manner. 

The first stochastic variable listed is called PNFIRE. 

This variable takes on a single value. This value is calculated by 
dividing the number of fires that have occurred by the number of 
days over which they have occurred. 

The variable PCOST is designed by cross tabulating the 
eleven suppression cost ranges with the fifteen discovery size 
categories. The results are the frequency distributions of fire 
suppression costs for each size class. The cumulative probability 
mass function can then be derived and the Monte Carlo method used 
to determine the suppression cost of each fire once the size at dis- 
covery has been established. It should be noted here that the cost 
ranges are replaced with an average cost for each range. The data 


set CAU describes these average values. 
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The variable PUBLIC describes a time distribution for 
the public detection system. The values used to build the proba- 
bility mass function for this variable are generated by cross tabulating 
the 25 time categories for each day with the 4 fuel categories. 

PFDIST values of the probability mass function are deter- 
mined by a three way cross tabulation. The fuel types are tabulated 
with the ranger districts holding the months constant. For example, 
in Footner, May will have five frequency distributions of fuels. 

One for each ranger district. 

The stochastic portion of the MONTH, MCODE, etc., set up 
deck is PS1 to 4 and PB1 to 5. The values for these probability 
mass functions are determined from cross tabulations of each day in 
the fire season with the 4 ranges of Spread index and the 5 ranges 
of the Buildup index. 

The probability mass function is determined for WHERE by 
cross tabulating the months with ranger districts. The frequencies 
of fires for each district in each month can then be determined. The 
probability mass function for PCDIST is determined by a three way 
cross tabulation. The tabulation of district with each fire cause 
holding the month constant produces the required frequency distri- 
butions. PTIME is determined by cross tabulating the 25 time cate- 
gories with the 10 cause categories. A frequency distribution of 


time is established for each of the 10 cause classifications. 
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VIS is determined by cross tabulating the 4 Spread 
index categories with the 6 visibility categories. The frequency 
distribution for the visibility categories is determined for each 
of the Spread index See corageee 

PFIRE frequencies are determined by a 3 way cross tabu- 
lation of day with cause and holding the month constant. This will 
give a frequency distribution of the number of fires occurring on 


each day of the fire season. 


Design Variables 


The design variables in the model are towers and air routes. 
The tower data is contained on two different sets of cards, the first 
set of cards is labeled IDC, TOWER], etc. . The second data set is 
labeled RING. The air route data set is labeled ROUTE. When modi- 
fying the design alternatives, either the tower data or air route 


data or both can be modified. 


Tower Data 

An example of an IDC, TOWER1, etc., card is shown in 
Figure A-3. The tower data is recorded in the following manner. 
The first two columns are reserved for the ranger district number. 
Columns 3 to 22 inclusive are reserved for the tower name. Columns 
23 and 24 are reserved for the tower number. Columns 25 to 27 inclu- 
sive are reserved for the tower's X coordinate and columns 28 to 30 


are for the Y coordinate. 
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FIGURE A-3 


An Example of the Set Up Deck IDC, TOWERI, etc. 
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Corresponding to this first set of data is a second set 
of data called RING. An example of a RING card is shown in Figure A-4. 
This is the blind area data for each tower in the first set. The 
blind area data must be in the same order as the tower data or else 
the blind area data will be applied to the wrong tower. Each tower 
has 20 cards describing its blind area. A tower has 4 quadrants 
associated with it and each quadrant has 5 concentric quarter 
circles at 5 mile intervals. (See Figure 1-3). One data card con- 
tains information on a sector of each quadrant. Columns 1 and 2 
are reserved for the quadrant number; column 3 for the RING number; 


columns 4 and 5 for the mileage of the RING from the tower; and 
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FIGURE A-4 


An Example of the Set Up Deck RING 
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columns 6 to 9 inclusive for the probability of not being able to 
see the fire. This probability is the blind area in the sector 
divided by total area in the sector. The RING portion of the tower 
design is stochastic. 

It should be noted that a tower called DUMMY had been added 
to the tower code names (See Appendix E). The purpose of the DUMMY 
tower was to allow the system to have a design alternative of zero 
towers. Since there must always be at least 1 tower to allow the 
program to run, a DUMMY tower located well outside the forest and 


having fictional RING data serves the purpose. 
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Air Route Data 
The air route data is located in the set up deck called 


ROUTE. An example of the set up deck is shown in Figure A-5. 


FIGURE A-5 


An Example of the Set Up Deck ROUTE 
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Each air route uses twelve data cards to describe it. The first 
card gives a general description of the air route. Columns 1 and 

2 are reserved for the ROUTE number, columns 3 to 7 inclusive repre- 
sent the 5 months May to September. If the ROUTE is to be flown in 
any of the 5 months, a 1 is punched in the corresponding column. 

If the ROUTE is not to be used during a month, a O is punched in 

the corresponding column. Columns 8 to 10 inclusive are reserved 
for the total mileage of the air route; columns 11 and 12 for the 

24 hour clock time the aircraft begins the flight; i.e. 1 p.m. = 13. 
Columns 13 to 15 are reserved for the aircraft speed in miles per 
hour; columns 16 to 18 for the cost in dollars per hour; and columns 
19 and 20 for the number of legs on the ROUTE. 

The next eight cards describe the air route shape and 
location. (See Figure A-6). The legs are counted clockwise from 
the starting point. Each time the ROUTE changes direction, a new 
leg starts. The program will accept up to 15 legs per ROUTE. Blank 
cards must be used to fill in the data for non-existing legs. For 
example, if there are 7 legs on a ROUTE, 4 blank cards will be 
required to complete the 8 leg cards. There are 38 columns used to 
describe an air route leg. The program is designed to describe the 
path of a ROUTE through a maximum of 9 ranger districts. There are 3 
columns set aside for each ranger district to record the length of 
the leg in that district. If the leg does not pass through a ranger 


district, the columns are left blank. The first 27 columns are used 
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Notice that for leg 2 because it is vertical, the intercept 


EXAMPLE OF A TYPICAL AIR ROUTE WITH 7 LEGS 


Not drawn 
to scale 


LEGEND 


~—— shows the starting and finishing 


coordinated, of each leg 


-——-shows where the Y intercept is 
measured for each leg 


measured is the “X intercept 


If 8 = 90, measure the X intercept, and the slope columns 


ares Lert blank. 


If @ = 0, measure the Y intercept and the slope columns 


are left blank. 


If @ = 0-90, measure the Y intercept and the slope is 


measured as tan @. 
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then, to record the length in miles of the leg in each ranger district. 
Columns 28 to 31 are reserved for the Y intercept (X intercept) of 
the leg. This is the distance in miles from the reference origin for 
the entire forest to where the leg, when extended, crosses the vertical 
(horizontal) axis. Which intercept to ide depends on the value of 0, 
see the note on Figure A-6. The reference origin is located so the 
forest being tested is in the first quadrant (see Figure A-6). Column 
number 28 is reserved for a + or —- sign depending on whether the inter- 
cept is positive or negative. Columns 32 to 36 are reserved for the 
slope of the leg. Again, this number is preceded by a+ or - sign in 
column 32. Columns 37 and 38 are reserved for 8, the smallest angle 
the leg forms with the horizontal. The value of 9 is recorded in 
degrees. Data for the next leg is repeated in this same order between 
columns 39 to 76. In this manner, data for each leg is recorded up 
to teg 15. 

The last three cards in the 12 card ROUTE deck are used 
to store data on the limitations of the air route previously described 
on the first 9 cards. Since the air route legs are described in terms 
of straight lines, of specific intercepts and angles, there must be 
some constraints placed on the boundaries of the routes, otherwise 
the leg will be viewed by the computer as being of infinite length. 

The last three cards, 72 columns on each of the first two 
and 36 on the last, are used to describe the X and Y limits for each 


leg. The first 12 columns are reserved for the starting X coordinate, 
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starting Y coordinate, finishing X coordinate and finishing Y coordi- 
mate of the first leg. The next 12 columns for the second leg and 

so on until data for the 15 legs is recorded. For non-existing 

legs the data card is again left blank. 

The air routes can be modified by changing the time the 
flight originates, the physical description of the legs or by 
adding new routes. The same air routes can be reflown later in the 
day by changing the air route number and the starting time and 
treating it as a completely new route. When the same air route is 
flown more than once, the descriptive statistics calculated showing 
the effectiveness of the air routes in the print out routine will 
be wrong, unless the variable NREPET is changed. The present pro- 
gram limits the designer to flying all the routes at the same 


frequency. The routes must be flown either all once or all twice etc. 
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SSIGNON GAVM *LIST® 
ON AT 23:244.5%4 CN 04-10-72 LAST ON AT 23:12.43 
SLIST *SCURCE* 


5 c FOBTRAN IV G MAIN PROGRAM 
6 
7 Cc FOCTNER FOREST FOREST FIRE DETECTION SYSTEM DESIGN MODEL 
8 
9 Cc THIS IS A PROTOTYPE DESIGN MODEL FOR THE ENTIRE ALBERTA FOREST SYSTEM 
10 
11 
12 ce THE LIMENSION STATEMENTS LISTED BELOW ARE FOR THE ENTIRE PROGRAM 
13 
14 INTEGER TITIME 
15 INTEGER TOWER, FUEL 
16 IRTEGER SPRING, SUMMER 
17 INTEGER SPRDEX,BUPDEX,SCOUNT, DIST, CAUSE,TIME 
18 1,ETRE1(4,5) ,BTABZ (4,5) , BURNDX,HORIZN 
19 DIMENSICN ATOWER (30) -PNFD (30) ,AROUTE (15) , PTFD (30) , PSUFA ( 
20 115) ,PTECA(15) ,AIFUEL (4) ,PGREAT (4) 
21 DIMENSION X(500) 
22 LIMENSION AMAY (15) ,REST (15) 
23 DIMEKSICN GT$({11) 
24 DIMENSICN BROUTE (15), PSNFA2(15) ,PTFDA2(15) 
25 CIMENSION PCOST(15, 11) ,AMONIQ (11) 
26 DIMENSICN CRANGE (11) ,PRANGE(11) ,CT$ (11) ,CAU (11) 
27 LCIMEXSION ROUTE(15,251) , IBOYD (15) , PUBLIC (4,25) 
28 DIMENSION RING (600,84) 
29 DIMENSION IDC (30), TOWER1(30) , TOWER2 (30) , TOWER3 (30) , TOWER4Y (30), 
30 TICKERS (30) , ITC (30) , LYCORD (30) , LXCORD (30) 
31 DIMENSION PEDIST (25,4) 
32 CIMENSION CDRANGE (11) 
33 DIMEKSICN RCS{10) ,PRCS(10) 
34 DIMENSION RSIZE(15) ,PSIZE (15) , PRSIZE (15) 
35 DIMFNSION RV(6),PRV (6) 
36 DINENSION MONTH(153) , IDAY (153) ,MCODE (153) ,PS1(153) ,PS2(153) ,PS3 (15 
a7 13), £S4 (153) ,PB1(153),PB2(153) ,PB3 (153) ,PB4 (153) ,PB5 (153) , WHERE (5,5 
38 1) 
39 DIMENSION SPOT(10,6) ,MCOORD (2) 
40 DIMENSION PCDIST (25,10) ,PTIME(10,25), 
41 1VIS(4,6),PIDAYS(153),PPIRE(153) ,SPREAD(&,5, 24) 
42 1,EDIST (5,10) -XN(24) ,XNORML (24) 
G3 CCMMON IX,IY,YFL 
Ty 
45 ; 
46 Cc THIS SECTION OF THE PROGRAM READS IN ALL THE DATA REQUIRED TO DSIVE 
a7 c THE SIMULATICN MODEL 
48 
Te) Cc NREPEI=1,2,OR 3 IF AIRCRAFT FLY ONCE,TWICE,OR THREB TIMES 
50 REAL (5,4441) NREPET 
51 4441 PCEMAT(I1) 
52 Cc NAIR IS EQUAL TO THE NUMBER OF AIRCRAFT REQUIRED 
53 BFAD(5,6201) NAIR 
Su 6201 FCRMAT(I2) 
55 Cc ZALEHA IS THE PERCENTILE OF THE NORMAL DISTRIBUTION WHICH LEAVES 
56 Cc ALPHA/2 PERCENT FROBABILITY IN EACH TAIL 
57 READ (S,3501) ZALPHA 
58 3501 FCEMAT(F5.4) 
59 c IX IS THE SEED POR THE RANDCM NUMBER GENERATOR 
60 
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FORTRAN IV G MAIN PROGRAM 


READ (5,3500) IX 

FCRMA1(I9) 

TIX = er 

CAU IS A VECTOR OF 11 AVERAGE COSTS 
REAL (5,2034) (CAU(J) ,J=1,11) 

FORMAT (11P6.0) 

CCSTT IS THE SEASONAL COST PER TOWER 
REAL (5, 2017) COSTT 

FCRMAT (F10.0) 


SPRING, SUMMER IS THE NUMBER OF AIR ROUTES OPERATING IN SPRING ASD SUMMER 


FERC (S,2006) SPRING, SUMMER 

FURMAT (212) 

ALT1 TO ALT10 IS THE DESIGN ALTERNATIVE FOR THIS RUN 

REAL (5,2003) ALT1,ALT2,ALT3,ALTU,ALT5,ALT6,ALT7,ALTS,ALT9, ALT10 
FCRMA1 (1044) 

PCRES1 TO FORESS IS THE NAME OF THE FOREST 

READ (S,2001) FORES1,FORES2,FORES3, FORESY,FORESS 

FOFMAT(SA4) 

ENFIEE IS THE PROBABILITY OF HAVING A FIRE FOR EACH DAY IN THE SEASON 
REAC(5,3017) PNFIRE 

PCBPAT(F10.5) 

READ THE FLYING CODE CRITERIA. IF MATT =0,CHECK WEATHER AND CAUSE. IP 
MATTI=1,ALWAYS FLY. IF MATT=2,NEVER FLY. 

READ (5,745) MATT 

REAL TEE CAUSE REQUIRED TO INITIATE AIRCRAFT 

FEAD(5,745) IC 

FCEMAT (12) 

READ THE MINIMUM BURNING INDEX REQUIRED TO INITIATE AIRCRAFT 
FEALC(S,745) IB 

READ TEE NUMBER OF AIRCRAFT ROUTES BEING USED. 

READ (S,748) NRCUTE 

FCRMAT (I2) 

REAL THE NUMBER CP TOWERS OPERATING IN THE SYSTEN 

REAL (S,797) NIOWER 

FORMAT (12) 

KEY= 20 * NTOWER 

READ TEE NUMBER OF FIRE SEASONS BEING CONSIDERED 

READ(S,648) ISEX 

PCBMAT (I4) 

REAC THE NUMBER OF DAYS IN A FIRE SEASON. 

REAL (5,1) NDAYS 

FORMAT (13) 

AMCNIQ IS A VECTCR OF AVERAGE SUPPRESSION COSTS 

FEAL (5, 1405) (AMONIO (J) ,J=1,11) 

FCBMAT (11F7.0) 


PCCST IS A TABLE OF FIRE SIZE AT DISCOVERY VS. AVERAGE SUPPRESSION COSTS 


REAL (5,1400) ((PCOST(I,J) ,J=1,11) ,I=1,15) 

FCRMAT (11F4.3) 

ROUTE IS A TABLE OF DATA DISCRIBING EACH AIR ROUTE 
REAC(S,720) ((ROUTE (I,J) ,J=1,251) ,I=1,NROUTE) 
PCEMAT (F2.0,5F 12 0, F300, F2.0,F 300, F 300,F220/ 
19F3.0,F4.0,F5.2,F220,9F 3.0, Fe 0g FSe2eF200/ 
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122 

2s 

124 

125 Cc FORTRAN IV G MAIN PROGRAM 

126 

127 

128 19F3.0,PF4.0,F5-2,F2-0,9F 3.0, F420, F5e2,F220/ 

129 19F3.0,F%.0,F5.2,F2.0, 9F 3.0, F4.0,F5e2,F2-0/ 

130 IIE CPE st Dec p kh le Oe GES Ope 60, ES. 2p hte Os 

131 19FO 50; F4s0 85627 52-0, 983.0, F420, F5.20P2.07 

132 1SF3.0,F4.0,F5.2,P2.0,9F3.0,F4.0,F5-2,F2-0/ 

133 NSFS.OC F400 F5Sai2,F2.0,9F3.0,F4.0,F5.2,F2-0/7 

134 19F3.0,F4.0,P5-2,P2-0/24F3.0/24F3-0/12F 3.0) 

135 € PUBLIC IS A TAPLE OF CUMULATIVE FROBABILITIES POR BURNING TIME 
136 C IF IFE FIRE IS DISCOVERED BY THE EUBLIC 

137 REAL (5,737) ( (PUBLIC (I,J) ,J=1,25) ,1=1,4) 

138 737 PCFMA1(16F5.4/9P5.4) 

139 c REAL THE LOCATICN NAME AND COORDINATES OF EACH RELEVANT TOWER 
140 CC 798 J = 1,NIOWER 

141 REAL (5,799) ICC(J) ,TOWER1(J) ,TOWER2 (J), TOWER3 (J) ,TOWER4 (J) , TOWERS ( 
142 13) ,I1C (J) -IXCORD(J) ,IYCORD (J) 

143 799 FCEMAT(I2,5A4,12,13,13) 

144 798 CCNTINUP 

145 Cc READ THE BLIND AREA DATA FOR EACH TOWER IN THE SYSTEM 

146 FEAD(5,795) ((RING(I,J) ,J=1,4) ,I=1,KEY) 

147 795 FCRMAT(F2.0,F1.0,F2.0,P4, 3) 

148 c READ THE CUMULATIVE PROBABILITY DATA FOR THE DIFFERENT FUELS 

149 Cc GIVEN THE DISTRICT AND MONTH OF THE FIRE 

150 REAL (5,300) ((PFDIST(I,3) ,J=1,4) gI=1, 25) 

151 300 FCRMAI(4F5S.4) 

152 Cc BEAD FOR EACH LAY IN THE MONTH, THE NAME OF THE MONTH,THE CODE FOR 
183 Cc THE MCNTH,AND THE CUMULATIVE PROBABILITIES OF BUILDUP INDEX 

154 Cc AND SESEAD INDEX FOR EACH DAY 

155 READ (S,2) (MONTH (3) ,MCODE(J),IDAY(J),PS1(J), PS2(J) ,PS3(J) ,PS4(J), 
156 1EE1(J),P22(J),PB3(J),PB4Y(J), PB5(J) ,J=1, 153) 

157 2 FCEMAT(A4,12,12, 9F4. 3) 

158 é REAL THE MONTHLY CODE AND THE CUMULATIVE PROBABILITIES FOR PIRES 
159 Cc CCCUFBING IN AXY GIVEN DISTRICT 

160 REAL (5,4) ((WHERE(I,J),J=1,5) ,I=1,5) 

161 4 POENAT (2X,5F5-4/2X,5F5.4/2X,5F504/2K yo SF504/2K, SF 504) 

162 c RFAC TiE CUMULATIVE PROBABILITIES FOR A FIRE BEING CAUSED IN A 
163 c FARTICULAR DISTRICT BY A GIVEN EVENT 

164 READ(5,5) ({PCDIST(I,J) ,J3=1, 10) ,L=1,25) 

165 5 FCRYAT(10F5.4) 

166 Cc REAL TEE MAP CATA NECESSARY TO SPOT THE FIRES 

167 READ(5,199) ((SPOT (I,J) ,J=1,6),1=1,10) 

168 199 FCRMAT(6F8.5) 

169 c READ LATA FOR CAUSE VS CUMU' ATIVE PROBABILITY OF TIME OF DAY OF FIRE START 
170 READ(5,10) ((PTIME(I,J) -J= 1,25) ,1=1,10) 

171 10 FORMAT (10F5.4/10F5. 4/5F5. 4) 

172 Cc READ VALUES OF BURNING INDEX IN TERMS OF BUILDUP INDEX AND SPREAD 
173 (e INLEX BTAB1 = SPRING AND SUMMER BTAB2 = FALL 

174 READ (5,23). ((BTAB1 (I,J) ,J=1,5) -1=1,4) 

175 REACAS. 23), ((BTAB2 (1,3) 5 J=1-5) ,1=1,4) 

176 23 FCBKMAT(SI1) 

177 c REAL CUMULATIVE PROBABILITY FOR VISIBILITY IN MILES ¥S SPREAD INDEX 
178 REAL(5,24) ( (VIS (I,J) ,J=1,6) ,I=1,4%) 

179 24 FORMAT (6F3. 2) 

180 
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FORTRAN IV G MAIN PROGRAM 


READ VALUES FOR RATE OF FIRE SPREAD WHERE I INDEX REFERS TO PUEL TYPE 

J INDEX REFERS TO THE BURNING INDEX AND THE K INDEX REFERS TO TINE OF DAY 
pC 33 I = 1,4 

DO 38 3 = 1,5 

READ (5,32) (SPREAD(I,J,K) ,K=1,24) 

FCFMAT (24F3. 1) 

CCNTINOE 
CCNTINUE 
pc’ 672 Ir= 
LO 671 J = 
CCNTINUE 
CONTINUE 
REAL THE CUMULATIVE PROBABILITIES FOR A FIRE OCCURRING ON ANY GIVEN DAY 
DURING THE 153 DAY PIRE SEASON 

REAL (5,30) (PPIRF(J),J=1, 153) 
FCBMAT(VOF7.6/10F7.6/10F7.6/10F7.6/10F7.6/10P7.6/10F7.6/10F7.6/ 
110F7.6/10FP7.6/10F7.6/10F 7. 6/10F 7.6/10F7<6/ 10F7« 6/3F7~ 6) 


~~ = 
ne 


THIS SECTION OF THE PROGRAMS ZERO"S ALL THE REQUIRED SEASONAL COUNTERS 


TAPAY=0. 
TTACFKE=0. 
IFEST=0. 
VAF=C. 
XPAB=0. 
SIGMA=0. 
D=C. 
TEIFC2=0. 
TESNF2=0. 
TFAL=0. 
TSNED=C. 
TENED=0. 
TFEIFC = 0. 
TENFA=0. 
TESNFA=0. 
TFEIFCA=0. 
SNFE=0. 
TCIALG=0. 
GREAT=0. 
CAYSA = 0. 
BANCE = 0. 
SUMPR = 0. 
Sas) = Os 
SHSIZE = 0. 
SFSIZE = 0. 
sfcs = 0. 
SFECS = 0. 
SKVY = 0. 
SPBV = 0. 


TACRES = 0. 
DO 4013 J=1,4 
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243 
244 

245 Cc FORTRAN IV G MAIN PROGRAN 
246 

247 

248 AIFUEL(J)=0. 

249 FGFEAT (J) =0. 

250 4013 CCNTINUE 

251 DO 4OC7 J=1,NROUTE 
252 BFCUIE(J)=0. 

253 AFCUTE (J) =0. 

254 PSNFA (J) =0. 

255 PTFLA (J) =0. 

256 FSNFA2(J)=0. 

257 PTFLAZ(J) =0. 

258 4007 CCNIINUE 

259 DO 4C00 J=1,NTOWER 
260 ATICWER(J)=0. 

261 ENED (J) =0. 

262 PIFL (J) =0. 

263 4000 CCNIINUE 

264 DO 2026; 3 = 1,11 
265 GI$(J) = 0. 

266 CT$.(J) = 0. 

267 CRANGE(J) = 0. 

268 PFANGE(J) = 0. 

269 2028 CCNIINUE 

270 Do 2053 I = 1,15 
271 AMAY(I) =0. 

272 REES1(I) = 0. 

273 A SIZE, (10. =F Or 

274 PSIZE(I) = 0. 

275 2053 CCKIINUE 

276 Dc z0£7 J = 1,10 
277 RCS (3): 806 

278 PECS (J) -= 0. 

279 2057 CCNTINUE 

280 Dc 2C€1 J = 1,10 
281 RV(J) = 0. 

282 FRV(J) = O. 

283 2061 CCNTIINUE 

284 

285 

286 Cc 8888 IS THE RANGE FOR THE SEASONAL CALCULATIONS 
287 

288 DO 8€&@ MMMM = 1,ISEX 
289 Cc 7€50 J=1,11 

290 CRANGE (J) =0. 

291 7650 CCNIINUE 

292 BTd=0. 

293 XXEAP=0. 

294 VVAR=0. 

295 AAY=0. 

296 JSEX = ISEX 

297 JSEY="AMM 

298 NCDEX = SSMM 

299 MCOUNT = 0 

300 
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302 
303 

304 

305 Cc FORTRAN IV G SAIN PROGRAM 

306 

307 

308 FIFE = 0. 

309 

310 
311 E THIS SECTION CALCULATES THE NUMBER OF FIRES IN THE SBASON 
EP. 

313 EG 28 J = 1,153 

274 c GREG IS THE NAME OF THE RANDOM NUMBER GENERATOR 

315 CALL GKEG 

316 IP(YFL .LE. PNFIRE) FIRE = PIRE # 1. 

Sid, 28 CCATINUE 

318 c TKIE SECTION ASSIGNS THE FIRES FOR THIS SEASON TO A PARTICULAR DAY 
319 

320 IFIFE = FIRE 

321 BGs 29 Je = eos 

322 29 EXDAYS ()) = 0 

323 DG BtI = LEB E 

324 CALL GFEG 

325 INDEX = 0 

326 rC 66 I = 1,153 

327 INDEX = I 

328 IF(I ~EC.1)GO TO 65 

329 IF(I.FC.153)GO TO 84 

330 IF (YFL.GT.PFIRE(I-1).AND.YFL.LE.PPIRE(I ))GO TO 67 
331 GC 10 66 

332 84 IF(YFL.GT.PPIRE(I-1))GO TO 67 

333 65 IF(YFL.LE.PFIBRE(I))GO TO 67 

334 66 CCKHIIKNUE 

335 67 FICAYS(INDEX) = FIDAYS(INDEX) + 1. 

336 31 CCNTINUE 

337 

338 

339 Cc STATEMENT 9999 IS THE RANGE FOR THE MONTHLY CALULATION BY DAY 
340 ; 

341 XKCUNT=0. 

342 CO 9599 KKKK = 1,NDAYS 

343 I2CB = 0 

344 KKK = KKKK 

345 IFICA = FIDAYS (KKKK) 

346 IF(IFICA .GT.0)GO TO 649 

347 IF (KKK. EQ. 1,.0R. KKK. EQ. 32. 0R. KKK. BQ.62.0R. KKK. EQe 93-0R- KKK. EQe 124) 
348 1G¢ 1¢ 3 

349 GC 10 2018 

350 

351 

352 Cc TRIS SECTION CALCULATES THE SPREAD INDEX FOR THE DAY 
353 

354 649 CALL GREG . 

355 IF (YFL.LE.PS1(KKKK))SPRDEX = 1 

356 IF (YFL.GT.PS1 (KKKK) -<AND.YFPL.LE.PS2 (KKKK) )SPRDEX = 2 
357 IF (YFL.GT.PS2 (KKKK) -AND.¥PL.LE.PS3 (KKKK)) SPRDEX = 3 
358 IF (YPL.GT.PS3 (KKKK) «AND YFL.LE.PS4 (KKKK)) SPRDEX = 4 
359 

360 
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362 

363 

364 

365 c PORTRAN Iv G MAIN PROGRAM 

366 

367 

368 Cc THIS SECTION CALCULATES THE BUILD UP INDEX FOR THE DAY 
369 

370 IF(YFL.LE.PB1(KKKK)) BUPDEX = 1 

371 IF (YFL.GT.PB1(KKKK) -AND.YFL.LE. PB2 (KKKK)) BUPDEX = 2 

Shere IF (YFL.GT.PB2 (KKKK) -AND. YPL.LE.PB3 (KKKK)) BUPDEX = 3 

373 IF (YFL.GT.PB3(KKKK) -AND.YPL.~LE.PBY (KKKK)) BUPDEX = 4 

374 IF(YFL.GT.PBY (KKKK) -AND.YFL.LE. PBS (KKKK)) BUPDEX = 5 

375 IF(IFCB .~.EQ. 1)GO TO 651 

376 

377 

378 (eS THIS IS A CHECK TO SEF IF THERE IS A FIRE ON THIS DAY 

379 

380 IF (KKK. EQ. 1-08. KKK. EQ. 32. OR. KKK. EQ.62.OR. KKK. EQ. 93.O0R. KKK. EQ. 124) 
381 NGC To" 3 

382 GO 1C 6 

383 3 MCCUNT = MCOUNT ¢ 1 

384 BPLIANE = S * (HCOUNT —- 1) 

385 IF(IFILTA .~EQ. 0)GO TO 2018 

386 

387 

388 Cc THIS SECTION LOCATES THE DISTRICT THE FIRE OCCURS IN 

389 

390 6 XKCUNT = XKOUNT *# 1. 

391 CALL GFEG 

392 IF(YFL.LF.WHERE(MCOUNT,1))DIST = 1 

393 IF (YFL.GT.WHERE(MCOUNT,1) -<AND.YFLLE.WHERE(MCOUNT,2))DIST = 3 
394 IF (YFL.GT. WHERE (MCOUNT, 2) -AND. YFL.LE. WHERE (MCOUNT,3)) DIST = 5 
395 IF (YFL.GT.WEERE(MCOUNT, 3) eAND-YFLJLE. WHERE (MCOUNT,4)) DIST = 6 
396 IF (YFL.GT. WHERE (MCOUNT,4) -AND.YFL.LE. WHERE (MCOUNT,5)) DIST = 7 
397 TPE(ETSI SCE Oem) orisor = 1 

398 ME ({CTST SEO so) nu esr 1=i 2 

399 LE (DISH SEOs oy) Le prsT =. 3 

400 IF(CIST .EQ. 6)IDIST = 4 

401 LP (DUST are co eres) => 

402 MAFK = MDIANE ¢# IDIST 

403 

4O4¥ 

405 (e THIS SECTION LOCATFS THE X AND Y COORDINATES FOR THE FIRE 
406 

407 IF (ICIST-1) 168,198,197 

408 198 CALL GFEG 

409 IF (YFL.GT.O.. AND. YFL.LE.. 404) NROW=1 

410 IF (YFL.GT..-4O4.AND.YFL.LE.. 835) NROW=2 

411 IF (YFL.GT..835.AND.YFL.LE..921) NEOW=3 

412 IF(YFL.GT..921.AND.YFL.LE.1.) NROW=4 

413 GCe TOs ngo 

414 197 IF (IDIST-3) 192,196,192 

415 196 CALL GREG 

416 IF (YFL.~GT.0..AND.YFL.LE.j. 242) NROW=6 

417 IF (YFL.GT. -242.AND.YFL.LEee S44) NEOW=7 

418 IF (YFL.GT. -SY4.AND.YFL.LE.1.0) NROW=8 

419 GO 10 155 

420 
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422 
423 
Y24 
425 
426 
427 
928 
429 
430 
431 
832 
433 
439 
435 
436 
437 
438 
439 
Gud 
4ut 
442 
yu3 
4ay 
Gus 
446 
4yu7 
448 
44g 
450 
4S1 
452 
us3 
454 
455 
456 
457 
458 
459 
460 
461 
462 
463 
464 
465 
466 
467 
468 
469 
470 
471 
472 
473 
G74 
475 
476 
477 
478 
479 
480 
481 
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PORTRAN IV G MAIN PROGRAM 


192 IF (ICIST. EQ. 2) NROW=5 
IF (IDIST.EQ.4) NROW=9 
IF (ILIST.EQ.5) NROW=10 
195 DO 130 KK=1,2 
MM=SECT (NROW, KK¢ 2) 
CALL GREG 
DC 131 J=1,Mm 
XJ=J 
INDEX=J 
CUMP 1=XJ*SPCT (NROW, KK) 
CUME2=CUMP1-SECT (NROW, KK) 
TEUIREG. 1) GON TOMI27 
IF (YFL.GT.CUMP2.AND.YFL.LE.CUMP1)GO TO 128 
GO 10 131 
127 IF(YFL.GT.0..AND.YFL.LE.CUMP1)GO TO 128 
131 CCKTINUE 
128 IF (INLEX.EQ.1)GO TO 129 
MCOCRD (KK) =SPOT(NROW,KK¢4) #INDEX 
GC TC 130 
129 MCOORD (KK) =SPOT(NROW, KK¢4) 
130 CCKTINUE 


THIS SECTION CALCULATES THE CAUSE OF THE FIRE 


CALL GREG 
INDEX = 0 
DC 7 J = 1,10 
INDEX = J 


EF (J. O5) 9 t)iGORTORS 

IF (J.EC.10)GO TO 83 

IF(YFL.GT. PCDIST (MARK, J-1) -AND. YFL.LE.PCDIST (MARK,J))GO TO 9 
GOGH 


83 IF(YFL.GT.PCDIST (MARK,J-1))GO TO 9 
8 IF(YFL.LE.PCDIST (MARK,J))GO TO 9 

7 CCKIINUE 

9 CAUSE = INDEX 


RCS (INDEX) = RCS(INDEX) *# 1. 


THIS SECTION CALCULATES THE FUEL THE PIRE IS BURNING IB 


CALL GREG 
INDEX=0 
po 302 J=1,4 
INDEX=J 
IF(J.EC.1)GO TO 304 
IF(J-FQ.4)GO TO 305 
IF (YFL.GT. PFDIST (MARK, J-1) «AND. YFL.LE.PFDIST(NARK,J)) GO TO 306 
Go 1G. 3C2 
305 IF(YPL.GT.PFDIST (MARK,J-1))GO TO 306 
304 IF(Y¥FL.LE.PFDIST (MARK,J))GO TO 306 
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485 Cc FORTRAN IV G NAIN PROGRAM 
486 

487 

488 302 CCKIINUE 

4389 306 FUEL=INDEX 

490 ALFUEL (INDEX) =AIPUEL (INDEX) #1. 

091 ¢ THIS SECTION CALCULATES THE TIME OF DAY THE FIRE STARTS 
492 

493 CALL GREG 

49a INDEX = 0 

495 BGpa ledh=0e25 

496 INDEX = J 

u97 IF(J-EC.1)GO TO 12 

498 IF(J-EC.25)GO TO 82 

499 IF (YFL.GT. PTIME (CAUSE, J-1) -AND. YFL.LE.PTIME(CAUSE,J ))GO TO 13 
500 6Ce tesa 1 

501 82 IF(YFL.GT.PTIME(CAUSE,J-1))GO TO 13 

502 12 IF(YFL.LE.PTIME(CAUSE,J)) GO TO 13 

503 11 CCKTINUE 

504 13 TIME = INDEX 

505 

506 

507 e THIS SECTION REASIGNS TIMES FOR PIRES THAT HAD UNKNOWN STARTING TIMES 
508 

509 IP(TIME .EQ. 25)GO TO 14 

510 GC 1c 15 

511 14 SUM = 0. 

512 LC 16:3) = 41,28 

513 SUM = SUM + PTIME(CAUSE,J) 

514 16 CCNTINUE 

515 fo 7 sents 

516 17 XN(J) = PTIME(CAUSE,J) /SUM 

517 BO 48sd5=— 17:24 

518 IF(J-EC.1)GO TO 19 

519 XNCRAL(J) = XN(J) # XNORAL(J-1) 

520 GC 1C0 18 

521 19 XNORBL(J) = XN(J) 

522 18 CCNIINUE 

523 CALL GREG 

524 INDEX = 0 

525 DOs 2063e=415.24 

526 INDEX = J 

527 IF(J.EQ.1)GO TO 21 

528 IF(J-EC.24)GO TO 81 

529 IF (YFL.GT. XNOBML (J-1) AND. YPL.LE.XNOFML(J ))GO TO 22 
530 GC 10 20 

$31 81 IF(YFL.GT. XNORML (J-1))GO TO 22 

532 21 LF(YEL.LE.XNOFML(J))GO TO 22 

533 20 CCBIINUE 

534 22 TIME = INDEX 

535 

536 

537 Cc THIS SECTION CALCULATES A BURNING INDEX USING THE SPREAD AND 
538 c BOILE UP INDEX TABLES 

539 

540 
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598 
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PORTRAN IV G MAIN PECGRAM 


IF (XKOUNT-—1.) 651,651,1090 

IF (MCODE (KKKK) -EQ.1.0R. MCODE (KKKK) - EQ. 2. O08. MCODE (KKKK) . EQ. 3) BURNDX 
J=E1AE1(SPRDEX,EUPDFX) 

IF (MCODE (KKKK) -EQ.4.OR.MCODE (KKKK) -BQ-5) BURNDX=BTAB2 (SPRDEX, BUPDEX 
1) 


THIS SECTION CALCULATES THE VISIBILITY FOR THE DAY 


CALL GREG 
IKDEX = 0 
DEt2sta°s 
INDEX = J 
IF(J-EC. 1)GO TO 26 

IF(J.EC.6)GO TO 80 

IF (YFL.GT. VIS (SPRDEX,J-1) -AND.YPL.LE.VIS(SPRDEX,J ))GO TO 27 
GEaT6e75 

IF (YFL.GT. VIS (SPRDEX,J-1))GO TO 27 

IP(YFL.LE.VIS(SPRDEX,J))GO TO 27 

CONTINUE 

HCRIZN = (INDEX * 5) - 5 

RV(INCEX) = RV(INDEX) ¢# 1. 

IE (IEBOB .EQ. 1)GO TO 2019 


1,6 


THIS SECTION DETERMINES WHETHEB THE TOWERS ARE CLOSE ENOUGH TO SEE 
wens fees 


INK=1 

CC 793 J = INK,NIOWER 

TNE EX eed) 

LENGTH= ( (IABS (MCOORD (2) -IXCORD (J) ) ) **2+# (LABS (NCOORD (1) ~1YCORD (J) )) 
1997) 99.5 

IPF (LENGTH. LE.HORIZN)GO TO 792 
IF(J.EC.NTOWER)GO TO 791 
CCNTINUE 

ICWFRF = INDEX 

IICWER=TOWER 

GO 10 790 

IPURNT=7777 

GC TO 750 


IF THE TOWERS ARE CLOSE ENOUGH, THIS SECTION CHECKS THE BLIND ABEA 
CF THE TOWER TO VERIFY THAT THE FIRE CAN BE SEEN 


IP (MCCOBD(2)-IXCCRD(TOWER) .GE.0.AND.MCOORD (1) -IYCORD (TOWER) .GE. 0) 
VICUAL=1 

IF (MCCORD (2) -IXCORD (TOWER) .LE.0. AND. MCOORD (1) ~LYCORD (TOWER) .GE.0)I 
icuac = 2 

IF (BCOORD (2) -IXCORD (TOWER) .LE.0. AND. MCOORD (1) ~LYCORD (TOWER) ~LE.0)I 
1CUAL = 3 
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FORTRAN IV G MAIN PROGRAM 


IP (HCCORD(2) -IXCCRD (TOWER) .GE.0. AND. SCOORD (1) ~-ILYCORD (TOWER) «LE. 0) I 
ICUAL = 4 

JCANNE = (20% (TOWER-1) ¢1)+#(S*IQUAD-5) ¢( (HORIZN/S) —1) 
CALL GEEG 

IF (YFL.LE.RING (JOANNE,4))GO TO 730 

ATCWEB (ITOWER) =ATOWER (ITOWER) #1. 

CCMEUTE THE BURNING TIME IF DISCOVERED BY TOWER 

1F (TIME.GT.6-.0R. TIME.LE.19) IBURNT=0 

IF (1IME.LE.6) LEURNT=7-TIME 

IF (TIME.GT.19) IBURNT=24-TIME+*7 

GCRTOM iSO 

LEURNI=7777 

IF (1CWFER-NTOWER) 1001,750,750 

INK = TOWER + 1 

GO 1C 1000 

IECB = 1 

GC 10 649 


GCI TC THIS POINT ONLY IF IPIDA = 0 
CEECK FOR FLYING CODE CRITERIA 


IF (MATT - 1)9998,2099,9999 
IF(BCENDX - I2)9997,2099, 2099 
IF(CACSE -— IC) 9999, 2099,9999 


GCT 10 THIS POINT ONLY IF ONE OR MORE FIRES ON THIS DAY 
THIS SECTION CEECKS THE FLYING CODE TO SEE IF AIRCRAFT FLY TODAY 


IF (MATT — 1)725,2099,727 

IF (ECENDX-IB) 1300,2016, 2016 

IF (CAUSE.NEIC)GC TO 727 

IF(XCOUNT - 1.$2€99,2099, 2096 

IF(CAUSE .EQ.IC) GO TO 2099 

GGiehC 3727 

crO 2020 I = 1,NROUTE 

IF (RCUTE (I, MCOUNT#1)-1.) 2020,2021,2021 
DAY$A = CAYSA ¢ (ROUTE(I,7) /ROUTE (I,9)) *ROUTE(I,10) 
AAS=AA$+# (ROUTE (I,7) /ROUTE (1,9) ) *BOUTE (I, 10) 
CCNTINUE 

IF (IFIDA.EQ.0)GO TO 9999 


THIS IS THE START OF THE AIR ROUTE ROUTIBE 


DC 749 J=1,NBOUTE 

IF (ROUTE (J, MCCUNT#+1)—1.) 747,746, 746 
IeCYD (J) =0 

GO TC 749 

IBCYE (J) =1 

CCRIIBUE 

KIBK=1 
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FORTRAN IV G NAIN PROGRAH 


4844 DO 743 J=KINK,NROUTE 
KDEX=J 
IF (LEOYD(J).E¢.0)GO TO 743 
MAXIEF=RKOUTE (J,11) 
DO 742 LEG=1,MAXIE 
IMAXIE=LEG 


CHECK FOR FLIGHT LEG ANGLE 


ITHETA=23¢ (LEG-1) *12 
IF (RCUTE(J,ITHETA) .EQ.0.)GO TO 741 
IF(ROUITE(J,ITHETA) ~FQ.90.)GO TO 740 


WE KNOW THE LEG IS SLOPED SINCE THE TWO PREVIOUS IF STATEMENTS WERE FALSE 
IF (ROUTE (J, ITHETA)-LE.45.)GO TO 739 
THIS SECTION IS FOR AN AIR ROUTE SLOPED GREATER THAN 45 DEGREES 


BACIAN=ROUTE (J, I THETA) /57. 296 
FIREY=MCOURD (1) 
XLINE= (FIPFY-ROUTE (J, ITHETA-2) ) /ROUTE (J, ITHETA-1) 
FIEEX=MCOOBD (2) 
CDELTAX=ABS (PIREX-XLINE) 
SIGHT=CELTAX*SIN (RADIAN) 
ISIGE1I=SIGHT 
JAB=19244* (IMAXIE-1) 
HAM=193¢4* (IMAXIE-1) 
KAP=194#4* (IMAXIE-1) 
RAP=1S5¢4* (IMAXIE-1) 
IF (HCRIZN-ISIGHT) 742,1451, 1451 
1951 IF (KOUTE(KDEX, JAM) «LE. PIREX.AND. PIREX. LE. ROUTE (KDEX,KAM) .OR. ROUTE 
1(KLEX,JAM) «CE. FIREX. AND. FIREX.GE. ROUTE (KDEX,KAM) «AND. BOUTE (KDEX, 
IMAM) .LF.FIREY.AND.sPIREY.LE. ROUTE (KDEX,NAM) .OR. ROUTE (KDEX, MAM) -GE. 
IPIREY.AND. FIREY.GE.ROUTE (KDEX,NAM))GO TO 724 
GC 10 742 


THIS SECTION IS FCR AN AIR ROUTE SLOPED LESS THAN 45 DEGREES 


739 FALIAN=ROUTE(J,ITHETA) /57-.296 
FIREX=MCOORD (2) 
YLINE=ROUTE (J, ITEETA-2) *ROUTE (J, ITHETA~1) *PIREX 
FIFEY="COORD (1) 
CELIAY=ABS (FIREY-YLINE) 
SIGHI=CELTAY*COS (RADIAN) 
ISIGKI=SIGHT 
JAR=19244* (IMAXIE-1) 
BAM=19344* (IMAXIE-1) 
KAP=19444* (IMAXIE-1) 
KAB=195¢4* (IMAXIE-1) 
IF (HCRIZN-ISIGHT) 742, 1452, 1452 
1452 IF (RCUTE(KDEX, JAM) «LE. PIREX. AND. PIREX.LE. ROUTE (KDEX, KAM) .OR. ROUTE 
1(KCEX,JAM) «GE. FIBEX. AND. PIREX.GE. ROUTE(KDEX, KAN) «AND. ROUTE (KDEX, 


96 


722 

723 

724 

725 =A eer a 

726 c FORTRAN IV G MAIN PROGRAM 

727 

728 ; 

729 1MAM) .LE.FIREY.AND. FIREY. LE. ROUTE (KDEX ,NAM) -OB. ROUTE (KDEX ,MAM) «GE. 
730 1FIBEY.AND. FIREY.GE.ROUTE(KDEX,NAM))GO TO 724 

734 GEnTC 742 

732 

733 c THIS SECTION IS FOR A VERTICAL AIR ROUTE , THETA=90 

734 

735 740 XLINE=ROUTE(J,ITHETA-2) 

736 FIEEX=MCOORD (2) 

737 SIGHI=ABS (XLINE-PIREX) 

738 ISIGHI=SIGHT 

739 MAM= 19344 (IMAXIE-1) 

740 NAS=195¢4* (IMAXTE-1) 

741 IF (HCRIZN-ISIGHT) 742,1453,1453 

742 1453 IF (RCUTE(KDEX,MAM) «LE. PIREY.AND. PIREY. LE. ROUTE (KDEX, NAM) «OR. ROUTE 
743 1(KLEX,MAM) .GE. FIBEY.AND.FIREY.GE.ROUTE(KDEX,NAM))GO TO 724 

744 CE RE WZ 

745 

746 c THIS SFCTICN IS FOR A HORIZONTAL AIR ROUTE , THETA=0 

747 

748 741 YLINE=ROUTE (J, ITHETA-2) 

749 PIFEY=MCOORL (1) 

750 SIGHI=ABS (YLINE-FIREY) 

751 ISIGET=SIGHT 

752 KAM= 19444 (IMAXIE-1) 

7153 JAM=192+4* (IMAXIE-1) 

754 IE (HCFIZN-ISIGHT) 742, 1454,1454 

7155 1454 IF (FCUTE(KDEX,JAM) «LE. PIREX.AND. FIREX. LE. ROUTE (KDEX, KAM) «OR. ROUTE 
756 1(KCEX,JAM) .GE.FIREX.AND. PIREX.GE. ROUTE (KDEX,KAM))GO TO 724 

757 742 CCETINUE 

758 743 CCKTIIRKUE 

759 IF (IEURNT. EQ.7777) GO TO 738 

760 c IF TEE IBURNT VALUE IS 7777, THE FIRE WAS NOT DISCOVERED BY TOWER 
761 IDELI=IBUENT 

762 GCETO"7'31 

763 

764 

765 Cc THIS SECTION DETERMINES WHICH LEG THE FIRE OCCURRED ON, CALCULATES ‘THE 
766 Cc VALUE CF IBURNA, AND ALSO CHECKS TO MAKE SURE THE PIRE STARTED BEFORE 
767 c TEE AIRCRAFT ARRIVED ON THE SCENE 

768 

169 724 IRANK=CIST-1 

770 RECUTE (KDEX) =AROUTE (KDEX) #1. 

771 AIRLDIS=0. 

772 TOTAL = Oe 

773 rO 723 L=1,IMAXIE 

774 JAR=12% (L-1) 

775 EG 722) N=1e9 

776 SARR=(N¢11) JAR 

777 AIFDIS=AIRDIS*ROUTE (KDEX, JARR) 

778 LLL=JARR*IRANK 

779 IF (LEQ. IMAXIEAND.Ne EQ. 1) TOTAL=ROUTE (KDEX, LLL) 

780 722 CONTINUE ; 


781 


723 


G442 


844 3 


729 


726 
727 


738 


V3'5 
736 
733 


734 
732 


8101 
8102 
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CCNTINUE 

AIRDIS=AIRDIS- (TOTAL/2.) 
BURNA=AIRDIS/ROUTE (KDEX, 9) 
CLCCK=FOUTE (KDEX,8) *BURNA 
ICLCCK=CLOCK 

IF (ICLOCK-TIME)44U42, 4443, 4443 
IF (KLEX.EQ.NROUTE) GO TO 727 
KINK=KDEX+1 

GC 10 4444 

IF (ICLCCK.GE. TIME) IBURNA=ICLOCK-TI ME 
IF (IEURNT.EC.0)GC TO 726 

IF (IEURNA.LT.IBURNT)GO TO 729 
ILELI=IBURNT 

Geter yen 

ILELI=IBURNA 

TEAL=1EAD#1. 

BFCUTE (KDEX) =BROUTE (KDEX) #1. 
GG ewe 

ILEL1=0 

GO) JIC) W7/3\1 

IBUENA=9999 

IF (IEURNT.ECQ.7777)GO TO 738 
IDELI=IBUENT 

GG fo wen 


THIS SECTION DETERMINES THE TIME IT TAKES THE PUBLIC TO 


97 


FIND THE FIRE 


CALL CREG 

DO 7326 K = 1,25 

INDEX = K 

TE KOE CRM )IGO -T05735 

IF (K~EC.25)GO TO 733 

IF (YFL.GT. PUBLIC (PUEL, K-1) .AND. YPFL.LE. PUBLIC (FUEL,K))GO TO 734 
GC 1C 736 

IF (YFL.LE. PUBLIC (FUEL,K))GO TO 734 

CCKNTIINUE 


IBURNP = 30 

GC 1¢C 732 
IBURNE = INDEX 
ILELT = IBURNP 
SNFEF=SNPP+1. 


THIS SECTION CALCULATES THE SIZE OF THE FIRE WHEN IT IS DISCOVERED 


ACFES = 0. 

IF(IDELT - 1)8100,8101,8101 
TLELT. =30 

GC 1C 8102 

IDELI = IDELT + 1 

DG 721 IGROW = 1,IDELT 


842 
843 
844 
845 
846 
B47 
848 
849 
850 
851 
852 
853 
854 
855 
856 
857 
858 
es59 
860 
861 
862 
863 
B64 
865 
866 
867 
868 
869 
670 
871 
872 
873 
874 
875 
876 
877 
878 
879 
880 
881 
882 
883 
884 
885 
886 
887 
888 
889 
890 
891 
892 
893 
894 
895 
896 
897 
893 
899 
900 
901 


8504 


8505 
8506 


1403 
1402 
1404 


652 
9999 


7699 
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ISTAR = IGROW - 1 

ITTIME = TIME # ISTAR ; 
IF(TTIIME.GT.24.AND.TTTIME.LT.49)TTITIME = TTTIME - 24 

IF(TITISE .GT. 4Y8)TTTIME = TITIME - 48 

ACRES = ACRES ¢ ((SPREAD(PUEL,BURNDX,TTTIME) ) **2)/(40. * 3.14728) 
CCATINUE 

IF (HCCUNT-1) 8501,86501, 8502 

TACRES=TACRES*#ACRES 

GC TO @503 

TIACRE=TTACRE*ACRES 


THIS SECTION CALCULATES TARE SUPPRESSICN COST OF THE FIRE 


IF (ACEESSLES O25) I=1 

IF (ACKES.G1..25.ANDSACRES.LE. 1.45) I=2 

IF (ACFES.GT. 1.45. AND. ACRES. LE.2.45) 1=3 

IF (ACBES.GT.2.45.AND.ACRES.LE. 3.45) 1=4 

IF (ACFES.GT. 3.45 AND. ACRES. LE.4.45) I=5 
IF(ACFES.GT.4.45.-AND.ACRES.LE.5.45) 1=6 

IF (ACRES.GT.5.4S.AND.ACRES.LE.6.45) I=7 

IF (ACEFS.GT.6.45. AND. ACRES.LE. 7.45) 1=8 

IF (ACFES.GT. 7.45 .AND. ACRES. LE.8.45) I1=9 

IF (ACKRES.GT. 8.45. AND.ACRES.LE.9.45) I=10 
IF (ACKES.GT. 9.45. AND. ACRES. LE. 10.45) I=11 
IF (ACRES.GT.10.45.AND. ACRES. LE. 25.00) I=12 
IF (ACKES.GT. 25.00. AND. ACBES.LE. 100.00) 1=13 
IF (ACFES.GT.1C00.00.AND. ACRES. LE.500.00) 1=14 
IF (ACFES.GT.5C0. 00) I=15 

IF(KCCUNT -—1)€&504,8504,8505 

AMAY({I) = AMAY(I) + 1. 

GC TO 8506 

REST UL) C=. REST (Db) tae 

RSIZE(I) = RSIZE(I) + 1. 

CALL GFEG 

CO 1402 J=1,11 

INDEX=J 

IF (J.FC.1)GO TO 1403 

IF(J.FC.11)GO TO 1404 

IF(YFL.GT. PCOST(I,J-1) .AND. YF) .LE. FCOST (I,J))GO TO 1404 
GC 10 1402 

IF(YFL.LE.~PCOST(I,J))GO TO 1404 

CCNTINUE 

SECCE1T=AMONIQ (INTEX) 

CEANGCE(INDEX) = CRANGE(INDEX) + 1. 

CRANGE (INDEX) =CRANGE (INDEX) +1. 

IP (FILAYS (KKKK)-XKOUNT) 652,652,6 
XKOUK1T=0. 

CCNTINUE 

bc 7699 J=1,11 

GT${J)=CAU (J) *LRANGE (J) 

CCNIINOE 
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902 

903 

904 

905 

906 c FORTRAN IV G SAIN PROGRAM 
907 

908 

909 DC 7698 J=1,11 

910 BIS=E1E+GT$ (J) 

911 7698 CCKIINUE 

912 Cc THIS SECTION CALCULATES THE COST TO OPERATE THE TOWERS FOR & SEASON 
S03 

914 WDLEX = NDEX 

915 WIOWER = NTOWER 

916 Tf = GIOWER * COSTT 

917 

918 Cc THIS SECTION IS RESERVED FOR THE CALCULATION OF DAMAGE 
Si'9 

920 c¢=0. 

921 X (NDEX) =TS*¢AAS+D$+BTS 

922 IP (KCEX.GE.50)GO TO 6499 

923 GC 10 &888 

924 

925 

926 Cc THIS SECTION DEFINES THE STOPPING RULE FOR THE MODEL AWD CALCULATES 
927 Cc THE STATISTICS REQUIRED 

928 

929 6499 CC 6498 J=1,NDEX 

930 XXEAF=XXEAR#X (J) 

931 6498 CCNIINUE 

932 XEAR = XXBAR/WCEX 

933 CO €4S7 J=1,NLTEX 

934 VARR= (X (J) -XEAR) *#2 

935 VVAR=VVAR#¢VARR 

936 6497 CCHIINUE 

937 VAK = VVAR/(WDEX - 1.) 

938 SIGMA=VAR**¥.5 

939 D = ZALPHA *(SIGMA / WDEX**.5) 
940 GDEX=NCEX 

941 BIGN=(ZALPHA**2* VAR) /D¥*2 

942 IF (GLEX.GE.BIGN) GO TO 6496 
943 IF(NDEX.EQ.200)GC TO 5540 

944 GO 10 5539 

945 S540 ZALFHA = (D*WDEX**.5) /SIGMA 
946 GC 1C 6496 

947 

948 

949 Cc THIS SECTION ELIMINATES AUTOCORRELATION IN THE RANDOM NUMBER GENERATOR 
$50 

951 5539 CALL GREG 

952 AX = IIX 

953 AIX = YFL * AX 

954 IIx = AIX 

95/5 IIx = IIX +3 

956 TIX~ =; (TIXK/2) Se2 ye A 

957 Ix = IIx 

958 CALL GREG 

959 8888 CCNIINUE 

960 
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THIS IS THE START OF THE PRINT ROUTINE FOR THE FIRST PAGE OF OUTPUT 


6996 WRITE (6,2000) 
WEITE (6,2002) FCRES1,FORES2, PORES3, FORES4, FORESS 
WRITE (6,2004) ALT1,ALT2,ALT3,ALT4,ALTS,ALT6,ALT7,ALTS,ALT9,ALT10 
WEITE(6,2084) 

2084 FCRMAT(1HO,48X,*CODES USED IN COMPUTER PROGRAM®) 
WPITE (6,2085) 

2085 FCEMA1T(//1H ,6X,*FOREST DISTRICT CODE*,14X,*FUEL CODE*,22X,*SPREA 
1f INLEX CODE",10X,*BUILDUP INDEX CODE‘) 
WRITE (6,2086) 

2086 FCRMAT(1H0,9X¢°1",5X,"DFI%,21X4°1%,5X_¢*CONTFEROUS*,19X,"1",6X, "0-1 
1°,19X,'°1°, 6X, '0-9") 
WEITE (6,2087) 

2C87 FCEMET(1H ,9X,°2",5X,"DF3",21X,°2% g5X,"MIXED WOOD',19X,"2"%,6X, "2-9 
1*,19¥,°2", 5X," 10-25") 
WRITE (6, 2088) 

2088 FCBMAT(1H ,9X,°3",5Xy"DF5%,21X,%3°,5X,"LICHEN AND MUSKEG*®,12X,'3', 
15X,'10-19",18X,"3",5X, 26-508) 
WEITE(6,2089) 

2089 FCRMAT(1H ,9%,°4", 5X, *DP6*,21X, 94" 5Xe*GRASS*,24X,°4*,5X,°20 PLUS! 
116X%,°4",5X,"51-100%) 
WHITE (6, 2090) 

2090 FCEMAT(1H ,9Xe°5%,5X,*DP7*,85X,°5*,5X,°100 PLUS") 
WRITE (6, 2091) 

2091 FCRMAT(//////1H ,6X,*BURNING INDEX CODE*®,16X,*MONTH CODE*,22X,°PIR 
1E CAUSE CODE",19X,*TOWER NAME CODE") 
WEITE (6,2092) 

2092 FCRMAT(1H0,6X%,°1%,6X,'1-20 LOW®,19X,"1%, 3K, °MAY®, 25X59 1", 5X, OTHE 
1R INDUSTRIES" ,12%,°1",4X, YATES") 
WRITE (€,2093) 

2093 FOEMAT(IN ,6X,5°2",5X,"21-4O MODERATE", 15X,°2%, 3X, "JUNE, 24X,°2", 5X 
1, *LIGETENING’,18X,°2",4X,*STEEN RIVER") 
WRI1E(6,2094) 

2094 FCRMAT(1H ,6X,°3%,5X,"41-60 HIGH", 19X,°3%, 3X," SULY®, 24X,°3*,5X,'SE 
TITLEMERT', 10X,°3",UX,*AMBER") 
WRITE (€,2095) 

2095 FORFAT(1H ,6X,°4",5X,°61-80 VERY HIGH*,14X,*4",3X,*AUGUST',22X,'4" 
1,5X,"WCOD INDUSTBIFS*,13X,°4",UX,"ADAIR") 
WRITE (6,5396) 

5396 FCRMAT(1H ,6X%,°5*,5X,"°81-100 EXTREME*,15X,°5", 3X, "SEPTEMBER", 19X,° 
15¢,5X,"RECREATION*,18X,'5*,4X,"ZAMA®) 
WRITE (€,2097) 

2097 FCRKMAT(1H ,73X%,°6"%,5X,"INCENDIARY*,18X,°6",4X,* BASSET’) 
WRITE (6,2098) 

2098 FCRMAT(1H »73X,"7*,5X,"PUBLIC PROJECT", 14X,°7', 4X,*WATT MOUNTAIN®) 
wEITE (€,5399) 

5399 FCRMAT(1H ,73X,°8*,5X, "MISCELLANEOUS KNOWN',9X,°8*,4X,* RAINBOW") 
WRITE (6, 3000) 

3000 FORMAT(1H »73X,°9%,5X,"MISCELLANECUS UNKNOWN*,7X,°9*,4X, "WHITESAND 
1S*) 
BRITE (6,3001) 
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1022 

1023 

1024 

1025 

1026 Cc FORTRAN IV G MAIN PROGRAM 

1027 

1028 

1029 3001 FOFMAT(1AH ,72X,°10°,5X,"RAILROAD*,19X,°10°, 4X," PONTON’) 

1030 WRITE (6,3002) 

1031 3002 FORMAT(1H ,6X,°AIR ROUTE CODE*,21X,*SECTOR CODE*,54X,°11°,4X, "FOGG 
1032 1Y PCUNTAIN') 

1033 WhITE (6, 3003) 

1034 3003 FCFRAT(1H ,106X,°12*,4X,* LAMBERT CREEK") 

1035 WRITE (6,3004) 

1036 3004 FCRMAT(1H ,6X,°1%,5X,*ROUTE 64 SPRING®,15X,°1°,5X,°FIRST QUADRANT 
1037 1*,43xX,°13",4X,*BUFFALO’) 

1038 WRITE (6,3005) 

1039 3005 FCKNAIT(1H ,6X,°2*,5X,*ROUTE 21 SUMMER*®,15X,°2°,5X,"SECOND QUADRANT 
1040 19,42x¥,°14", 4X, "WADLIN®) 

1041 WRITE (6,3006) 

1042 3006 FCEMAI1(1H ,6X,°3*,5X,*ROUTE 22 SUMMER',15X,°3*,5X,*°THIRD QUADRANT 
1043 49, 42xX,°15",4UX,*PANNY®) 

1044 WRITE (6,3007) 

1045 3007 FCRKMAT(1H ,6X%,°4"*,5X,*ROUTE 23 SUMMER*®,15X,°%°,5X,*POURTH QUADRANT 
1046 1*,43X,°16",4X,*TALBOT LAKE") 

1047 WRITE (6,3008) 

1048 3008 FCERMAT(VH ,6X,°5%,5X,"ROUTE 24 SUSMER®,79X,°17°,4X,*BISON LAKE*) 
1049 WRITF(6,3009) 

1050 3009 FCKMAT(1H ,106X%,°18%,4X,"HAWK HILLS®) 

1051 WEITE(6,3010) 

1052 3010 FCKMAT(1H ,106X,°19°,4X,*JEAN LAKE") 

1053 WRITE (6, 3011) 

1054 3011 FCFMAT(1H ,41X,*RING CODE',53X,°20°,4X,* WOLVERINE’) 

1055 WRITF (6,8300) 

1056 8300 FCEMAT(1H ,106X,°21°,4X,* EETITOT’) 

1057 WRITE (6,3012) 

1058 3012 FCRMAT(1H ,42X,°1°,6X,°S MILE RADIUS! ,4UX,*22%, 4X, *DUMBY ‘) 

1059 WRITE (6,3013) 

1060 3013 FCERAT(1H ,42X,°2",5X,°10 MILE RADIUS") 

1061 WRITE (6,3014). 

1062 3014 FCREAT(1H ,42X,*3*,5X,°15 MILE RADIUS") 

1063 WRITE (6,3015) 

1064 3015 FOEMAT(1H ,42X,°4*,5X,*20 MILE RADIUS") 

1065 WRITE (€,3016) 

1066 3016 FOBMAT(1H ,42X,°5*,5X,'25 MILE RADIOS") 

1067 

1068 

1069 Cc THIS IS THE START OF THE PRINT RCUTINE FOR THE SECOND PAGE OF OUTPUT 
1070 

WO WRITE (€,2000) 

1072 2000 FCRMAI(1H1,30X,*PROTOTYPE FOREST FIRE DETECTION DESIGN MSODEL® 
1073 1) 

1074 WRITE (6,2002) FORES1,FPORES2,FORES3,FORES4,FORESS 

1075 2002 FCBBAT(1H ,41X,5A4,2X,*FOREST') 

1076 WRITE (6,2004) ALT1,ALT2,ALT3,ALTY,ALTS,ALT6,ALT7,ALTS,ALTI,ALT10 
1077 2004 FORMAT(1H ,13X,°DESIGN ALTERNATIVE - - - ~ - = ------. *,10A 
1078 14) 

1079 WEITE(6,2005) 

1080 2005 FCBMAT(1H ,41X,*DESIGN CHARACTERISTICS®) 


1081 
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2010 WRITE (6,2011) NTCWER, SPRING, SUMMER, NAIR 
2011 FCRMAT(1H ,3X,"NUMBER OF TOWERS °,14,3X,*NUMBER OF AIR ROUTES SP 
TFING ",14,2X,"SUMMER *,1I4,2X,*MAXIMUM NUMBER OF AIRCRAFT REQUIRE 
10°, 1€) 
WHITE (6, 2012) JSEX 

2012 FCRKMAT(1H ,29X,"NUMBER OF PIRE SEASONS RUN BEFORE CUTOPF®,110) 
WRITE (€, 2013) 

2013 FCFMAT(1HO,39X,° FOREST cost SUMMARY‘) 

IF (KICWER -0) 2015,2014,2015 

2014 AYG$1T = 0. 

Gorie (2022 

2015 ZICWER = NTOWEB 
AVG$T = ZTOWER * CCSTT 

2022 ALEX = JSEX 
AVG$A = CAYSA/ADEX 
LESIGN = AVG$T + AVGS$SA 
WRITE (6,2023) DESIGN 

2023 FCRHEAT(1H ,44X,*AVERAGE TOTAL COST OF THIS DETECTION DESIGN* 

1,14X,"$",F12.0) 
WRITE (6,2024) 

2024 FCEMAIT(1H ,43X,*SUPPRESSICN COST SUMMARY‘) 
WFITE(6,2025) 

2025 FORMAT(1H ,7X,'FIRE*,6X,°$") 

WEITE(6,2026) 

2026 FCEMAT({1H ,4X,*COST*,10X, *CATEGORY',5X,*NUMBER® ,22X, *CATEGORBY') 
WHITE (6, 2027) 

2027 FORMAT(1H ,4X,*RANGE*,5X,*AVERAGE OSED*,6X,'OP FPIRES*,3X,* PROBABIL 
V11Y',6X, "TOTAL $*) 
po 2c2o/,di= AuMA 
BANGE = RANGE ¢ CRANGE(J) 

2029 CCATINUE 
cer 2030, 3 
FRANGE (J) 

2030 CCHTINUE 
BEtzesa's = 4,71 
SUMEF = SUMPB ¢# PRANGE(J) 

2031 CCKIINUE 
PC 2032 3 =a1ga4 
C1£(J) = CAU(J) * CRANGE(J) 

2032 CCKIINUE 
DC 20323 J = 1,11 
ST# = ST$ + CT$(J) 

2033 CCAHTINUE 
WRITE (6,2035) CAU(1),CRANGE(1),PRANGE(1) ,CT$(1) 

2035 FCRMAT(1H ,2X,°NO CATA’, 7X,F7.0,7X,F7.0,6X, 75.4, 9X, P90) 
WRITE (6, 2036) CAU (2) ,CRANGE (2) , PRANGE (2) ,CT$ (2) 

2036 FCRMAT(1H ,5X,_°1-100",6X,P 720, 7X pF 70 0¢6X, F544 9X, F920) 
WRITE (6,2037) CAU(3) ,CRANGE(3) ,PRANGE (3) ,CT$(3) 

2037 FCRMAT{1H 43X,°101-500",6X,P7.0, 7X F720, 6X, F5. 4,91 4F9e0) 
WHITE (6,2038) CAU(4%) ,CRANGE (4), PRANGE (4) ,CTS (4) 

2038 FCEMAT (1H 6 3X%4°501-10009 5X FTO 4 IX oP 1009 6X 5 F504 ,9X, FI.0) 
WRITE (6,2039) CAU (S) ,CBANGE (5) , PRANGE (5) ,CTS$ (5) 
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an 

2039 PCRKPAT(1H ,2X,°1001-5000',5X,P7.0,7X,F720,6X,F5.4,9X,F9. 0) 
WRITE (6,2040) CAU (6) ,CRANGE (6) , PRANGE (6) ,CT$ (6) 

2040 FOCFMAT(1H ,2X,°5001-10000°,4X,F7.0,7X,F7.0, 6X, F5.4,9%,F9.0) 
WRITE(6,2041) CAU(7),CRANGE(7) ,PRANGE (7) ,CT$ (7) 

2041 FORMAT(1H ,1X,*10001-25000',4X,F7.0,7X,F7.0,6X,F5.4,9X »P9.0) 

WRITE (6,2042) CAU (8) ,CRANGE (8) -PRANGE (8) ,CT$(8) 

2042 PCEMAT(1H ,1X,°25001-50000', 4X, F7.0,7X,F7.0,6X,F5.4,9X%,F9.0) 
WRITE (6,2043) CAU (9) ,CRANGE (9) , PRANGE (9) ,CT$ (9) 

2043 FOCEMAT(1H ,1X,°50001-75000' ,4X,F7.0,7X,F7.0,6X,F5.4,9X,F9.0) 

WRITE (6,2049) CAU(10) ,CRANGE(10) ,PRANGE(10) ,CT$ (10) 
2044 FOFMAT(1H ,1X,°75001-100000",3X,F7.0,7X,F7-0,6X,F5.4,9X,F9.0) 
WPITE(6,2045) CAU(11) ,CRANGE(11) ,PRANGE(11) ,CT$(11) 
2045 FCKMAT(1H ,1X,*1C00000 + *,6X,F7.0,7X,F7.0, 6X¢F5.4,9X,F9.0) 
WRITE(6,2046) RANGE,SUMPR,STS$ 

2046 FORMAT (1HO,19X,*TOTALS® ,UX,F8.0,5X,P6.4,7X,°$*,F 10.0) 
BETS = STSE/SALEX 
WEITE(6,2047) ACTS 

2087 PCEMAT(1H ,44UX," AVERAGE TOTAL COST OF SUPPRESSION FOR THIS DESIG 
IN*,12X,°'S",F12.0) 

AC1IC = 0. 
WRITE (6,2048) ACTD 

2048 FCEMAT(1H ,44X,* AVERAGE TOTAL COST OF DAMAGE FOR THIS DESIGN®,17 
1xX,°$*,F12.0) 

BIGICT = DESIGN # ACTS # ACTD 
WREITE(E,5395) BIGTOT 

§395 FCRMAT(1HO,35X,"AVERAGE SEASONAL TOTAL COST OF DETECTION + SUPPRES 
ISICK ¢ DAMAGE* ,9X,°$*,F12.0) 

WRITE (6,2049) 

2C49 FCEMAL(//1H J4X,*FOREST*,3X,*FIRE® ,4X,*SIZE* ,4X,*SUMMARY®,6X,* FORE 
IST*, 7X, "CAUSE® ,9X,*SUMMARY® ,10X,*FOREST® ,5X,*VISIBILITY*, 4X, *SUMNMA 
VEY") 

WRITE (6,2050) 

2050 FCEMAT(1H0,3X,*SIZE*,9X,* NUNBER® , 3X, *PROBABILITY*,8X, "CAUSE", 15X,° 
INUMBEP IN*,6X, *PROBABILITY*,9X,"VISIBILITY*,2%,*NOMBER*,2X,°PROB 
VAEILITY") 

WRITE (6, 2051) 

2051 FCRPAT(1H ,* CATEGORY IN FOR EACH CATEGORY 
1 CATEGORY FOR EACH IN OF 
1 FCR EACH’) 

WRITE (6,2052) 

2052 FCRMAT(1H ,3X,*ACRES*,8X,*CATEGORY',3X, *CATEGORY* , 46X, "CATEGORY *, 1 
11X,"MILES*,7X, "DAYS", 6X, * CATEGORY") 
re 2054 J = 1,15 
SRSIZE = SRSIZE ¢ RSIZE (J) 

2054 CCNTINUE 
re 2055 J = 1,15 
ESIZE{J) = RSIZE(J) / SRSIZE 

2055 CCBIINOE 
co 2086 J = 1,15 
SPSIZE = SPSIZE # PSIZE(J) 

2056 CCETINUE 
bo 205@€ J = 1,10 
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SRCE = SRCS ¢* RCS(J) 
2058 CCKIINUE 
ECs 205900. =e deplO 
FRCS(J) = RCS(J) / SRCS 
2659 CCKIIKUE 
DC 2060 J = 1,10 
SERCS = SPRCS + PRCS(J) 
2060 CCNIINUE 
DG 20€2 J 
SRV = SRV 
2062 CCETINUE 
Berek ayd = 1.6 
EBV (J) =RV(J) /SRV 
2063 CCKIINUE 
CC 2Cé4 J = 1,6 
SERV = SPRV # PRV(J) 
2064 CCKTINUE 
WEITE(6,2065) RSIZE(1) ,PSIZE(1) ,RCS (1) »PRCS (1) ~RV(1) -PRV(1) 
2065 FCRRAT(1H0,UX,"SPOT*,8X,F7.0,4X,F 6.4, 10X,*OTHER INDUSTRIES",5X, F9. 
10,9X%_, F604, 12X," 0°,8X,F7.0,6X,F6.4) 
WRITE (6, 2066) RSIZE (2) , PSIZE (2) , RCS (2) »PRCS (2) -RV (2) »PRV (2) 
2066 FCRBMAT(1H 4 3X,% 25-1.45"',5X,F7.0,4X,F 6.4, 10K," LIGHTENING’, 11X,F9.0 
1,9X,F6.4,12X," 5°,8X,F7.0,6X, F624) 
WRITE (€,2067) RSIZE(3) , PSIZE (3) , RCS (3) »PRCS (3) ,RV(3) ,PRV (3) 
2067 FCRMAT(VH ,2X_1-46-2.45",5X,F7.0,4X,F6.4, 10K, "SETTLEMENT", 11X,F9. 
10,9X,F6.4U,12X,°10",8&X,F 7.0, 6X, F624) 
WHITE (6, 2068) REIZE (4) ,PSIZE (4) , RCS (4), PRCS (4) ,RV(4) ,PRV (4) 
2068 FCRMAT(1H ,2X,°2.46-3.45",5X,P7-0,4X,F6.4,10X,"WOOD INDUSTRIES", 6X 
1,F.0,S5X,F6.4,12X,'15%,8X,F 7.0, 6K, F6.4) 
WFITE (6,2969) REIZE (5) , PSIZE (5) » RCS (5) , PRCS (5) , RV(5) , PRV (5) 
2069 FCRMAT(1H ,2%,°3.46-4LU5",5X,P7.0, 4X, F6.4, 10X,* RECREATION, VIX, FI. | 
10,9X,F6.4,12X,*20',8X,F7.0,6X,F6.4) 
WHITE (6,2070) RSIZE (6) ,PSIZE (6) , RCS (6) , PRCS (6) , RV (6) , PRY (6) 
2070 FCKRAT(1H ,2X%-°4.46-5.45%,5X,F7.094X,P6.4, 10X,* INCENDIARY", 11X,F9. 
10,9X,_F6.4,12X,'25" 8X, F720, 6X,F 6.4) 
WRITE (6,2071) RSIZE (7) ,PSIZE (7) » RCS (7) , PBCS (7) 
2071 FCFMAT(1H 2% 9'5246-6.45",5X,F7.0,4X,F6.4,10X,"PUBLIC PROJECT', 7X 
1,FS.C,9X,F6.4) 
WEITE (6, 2072) RSIZE (8) »PSIZE (8) , RCS (8) , PRCS (8) 
2072 FCEMAT(1H ,2X,'6.45-7.45" ,5X,F7.0,4X,P6.4, 10X," MISCELLANEOUS KNOWN 
1°,2X,F9.0,9X,FEs 4) 
WRITE (6,2073) RSIZE(9) , PSIZE (9) , RCS (9) , PRCS (9) 
2073 FCRMAT(1H ,2X%,°7.46-8.45",5X,F7.0,4X,F6. 4, 10X,"MISCELLANBOUS UNKNO 
WN", FS.0,9X,F6.4) 
WRITE (6,2074) RSIZE (10) g PSIZE(10) ,RCS(10) , PRCS (10) 
2074 FCEMAT(1H 42X,°8.46-9.45",5X oF 70 0,4X,F6.4, 10%," RAILROAD® ,13X,F9.0, 
19X, F€.4) 
WRITE (6,2075) RSIZE(11) ,PSIZE(11) 
2075 PCRMAT(1H ,2X%,°9.46-10.45", UX, P7.0,4X,F6.4) 
WHITE (6, 2076) RSIZE (12) , PSIZE (12) 
2076 FCBMAIT(1H ,1X,5°10.46-25.00*, 4X, F7.0,4X,F6. 4) 
WHITE (6, 2077) BSIZE (13) , PSIZE (13) 
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2077 FCRMAT(1H ,1X,*25.01-100.00', 3X,F7~.0,4X,_ F6.4) 
WRITE (6,2078) RSIZE(14),PSIZE (14) 
2078 FCRMAT(1H ,°100.01-500.00°,3X,F7.0,4X,F6~4) 
WRITE (6, 2079) RSIZE (15) , PSIZE(15) 
2079 FCKMAT(1H ,*5C0.CO #*,8X,F7.0,4X,F6.4) 
WRITE (6, 2080) SRSIZE,SPSIZE, SRCS, SPRCS,SBV,SPRV 
2080 FCEMAT(1H0,6X "TOTALS, 3X ¢F8.0,4X, F604, 22K, "TOTALS", 2X,F10.0,9K,P6 
154,11X, *TOTALS*, 4X ,F8.0,6X,F6.4) 
DC e507 J=1,15 
TAMAY=TAMAY+AMAY (J) 
TREES T=TREST*REST (J) 
8507 CCRIIRUE 
APSAD=TACRES/TAMAY 
EISS=1TACRE/TREST 
WEITFE (6,2081) AFSAD 
2081 FCERMAT(1H ,2X,'MAY AVERAGE PIRE SIZE AT DISCOVERY® ,P8.2,2X, °ACRES 
1") 
WEITE(6,8508) PISS 
8508 FCEMAT(1H ,2X,"AVERAGE FIRE SIZE AT DISCOVERY EXCLUDING MAY FIRES* 
1,F8.72,2X,' ACRES") 
IF (NICWER.EQ.1)GO TO 9297 
LC 4CC1 J=1,NTCWER 
TONFL=TSNFD+ATCWER (J) 
4001 CCNIJKUE 
LC 4CC2 J=1,NTCWER 
ENED (J) =ATOWER (J) /TSNED 
4002 CCNTINUE 
CC 4003 J=1,NTOWER 
TENFL=TENFD+ENED (J) 
4003 CCKTINUF 
CC GOC4Y J=1,NTCWER 
FIED (J) =ATOWER (J) /SRSIZE 
4004 CCKTINUE 
CC 4OCS J=1,NTOWER 
TEIEL=1PTPD¢PTED (J) 
4OOS CONTINUE 
9297 IF(MAT1~EQ.2)GC TO 9299 
DC 4006 J=1,NRCUTE,NREPET 
ISNFA=1ISNFA*ARCUTE (3) 
4006 CCNIINUE 
DO 4CC& J=1,NRCUTE,NREPET 
FENFA (J) =AROUTE(J)/TSNFA 
4008 CCNTINUE 
CC 4009 J=1,NROUTE,NREPET 
TESNFA=TPSNFA+*ESNEA (J) 
4009 CCHTINUE 
[CG 4010 J=1,NROUTE,NREPET 
EFIFDA (J) =AROUTE (J) /SRSIZE 
4010 CCKTINUE 
CO 4011 J=1,NROUTE,NREPET 
TEIFLCA=IPTFLA*PTELA (J) 
4011 CCHTINOE 
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9299 TETFCP=SNPP/SRSIZE 
CC 4012 J=1,4 
ITCTALG=TOTALG AI FUEL (J) 
4012 CCNTINUE 
DC 4014 J=1,4 
PGREAT (J) =AIFUEL (J) /TOTALG 
4014 CONTINUE 
DO 4015 J=1,4 
GREAT=GREAT+PCBEAT (J) 
8015 CCKTINUE 


THIS SECTION STARTS THE PRINT ROUTINE FOR THE THIRD PAGE OF OUTPUT 


BRITE (6,2000) 
WRITE (6, 2002) PORES1, PORES2,FORES3,FORESY, PORESS 
WEITE(6,2004) ALT1,ALT2,ALT3,ALT4Y,ALTS,ALT6,ALT7,ALTS,ALT9,ALT10 
WRITE (€,4016) 

4016 FOBMAT(1HO,48X,*TOWER DISCOVERY SUMMARY‘) 
WEITE (6,4017) 

4017 PCEKAT(1HO,9X,*TOWER®, 10X,*TOWER*®,18X,"NUMBER CF ‘*,*FIRES DISCO 
WEREL*',7X,*PERCENT OF TOWER*,5X,*PERCENT OF TOTAL‘) 
WRITE (6,4018) 

4018 FCFMAT(1H ,10X,"CODE',11X,"NAME® ,57X,*DISCOVERBIES' ,9X,*FIRE DISCOV 
1EFIES',//) 
TC 4019 J=1,NTOWER 
WRITE (€,4020) ITC (J) , TOWER1(J) ,TOWER2 (J) ,TOWER3 (J) ,TOWERY (J) , TOWERS 
143), ATCWER(J) ,ENEL (J) ,PTPD (J) 

4020 FCRMAT(1H ,11X,12,5X,5A4, 21X,P7.0,23X,F7.2, 15X,F 7-3) 

4019 CCANITIRNCE 
WRITE (6,4021) ISNFD,TPNED,TPTFD 

4027 PORMAT(1HO,Y7X,*TOTAL® ,6X,F8.0,9X, ‘TOTAL', 9X, F7.2, 7X, "TOTAL", 3X,F7 
1.3) ‘ 
WRITE (6,4022) IPTED 

8022 FCRKMAT(1HO,72X,'PERCENT OF TOTAL FIRES DISCOVERED BY TOWERS',9X,F7 
1.3) 


THIE SECTION STARTS THE PRINT ROUTINE FOR THE FOURTH PAGE OF OUTPUT 


WRITE (6,2000) 
WRITE (6,2002) FCRES1, PORES2, PORES3, FORESS, PORESS 
WRITE (6,2004) ALT1,ALT2,ALT3,ALTY,ALTS,ALT6,ALT7,ALTS,ALT9,ALT10 
WEITE (6,4023) 

4023 FORMAT(1HO,47X,*AIRCRAPT DISCOVERY SUMMARY‘) 
WRITE (€,4024) 

8024 FOFMAT(1HO,26X,"AIR ROUTE*,6X,*NUMBER OF PIRES*,6X,*PERCENT OF AIR 
1 RCUTE',5X,*PERCENT OP TOTAL‘) 
WRITE (6,4025) 

4025 FORMAT(1H ,28X,*NUMBER* ,6X,*OCCURRING ON ROUTE*,10X, °OCCURRENCES', 
18%, "FIRE OCCUBBENCES',//) 
CC 4026 J=1,NROUTE,NREPET 
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WRITE (6,4027) ROUTE (J,1) ,AROUTE(J) ,PSNPA(J) , PTFDA (J) 
4027 FCEMAT(VH ,29X,F4.0,13X,F7.0,17X pF5.2517X,F5~ 2) 
4026 CCHTINUE 
WRITE (6,4028) TISNFA,TPSNFA,TPTFDA 
GO28 PCRMAT(1H ,38X,*TOTAL* ,2X,F8.0,8X,"TOTAL®, 2X,F7.2,9X,*TOTAL', FS. 2) 
WRITE (€,4029) TETFLA : 
4029 FCRMAT(1H ,66X,*PERCENT OF TOTAL FIRES OCCURRING ON AIR ROUTES", 3X 
VA AAA) 
IF (MATI.LEQ.2)GC TO 9298 
LC 5£98 J=1,NRCUTE 
ESNEAZ (J) =BROCTE (J) /TEAD 
5598 CCKTINUE 
LO S597 J=1,NROUTE 
PIFDCAZ (J) =BROUTE (J) /SRSIZE 
5597 CCKNTINUE 
9298 FEPFCI=1TEAD/SRSIZE 
WRITE (6,4024) 
WRITF(6,4025) 
Cc 5599 J=1,NROUTE 
WRITF(€,4027) KOUTE(J,1) ,BROUTE (J) ,PSNFA2(J) ,PTFDA2 (J) 
5599 CCNTINUE 
LC $596 J=1,NRCUIE 
TESNF2=TPSNF2¢ESNFA2 (J) 
5596 CCNIINUE 
CO 5595 J=1, NRCUZE 
TPIFDZ=IPTFD2¢+ETELA2 (J) 
5595 CCKTINUE 
WEITE(6,4028) TEAL, IPSNF2,TPTFD2 
WEITE (€,5594) PEFFCT 
5594 FCRMA1(1HO,€2X,*EFFFCTIVE PERCENT OF FIRES DISCOVERED BY AIR ROUTE 
1S",11X,F7.3,///) 
WRITE (6,4030) 
6030 FCERAT(1HO,46X%,*PUBLIC DISCOVERY SUMMARY") 
WRITE (6,4031) SNFE,IPTFDP 
4031 FCRMAT(1H0,9X,"TOTAL NUMBER OF FIRES DISCOVERED BY PUBLIC®,9X,F9.0 
1,3X,"‘PERCENT CF ICTAL FIRES DISCOVERED BY PUBLIC’,9X,F7. 3) 
SUMMY=TPTEFD+PEFFCT+TPTPDP 
WRITF (6.4932) SUMMY 
4032 FCRRAT(1HO,112%, "TOTAL", 7X,F7.~ 3) 


THIS SECTION STARTS THE PRINT ROUTINE FOR THE FIFTH PAGE OF OUTPUT 


WRITE (6,2000) 
WFITE (6,2002) PCRES1,FORES2, FORES3, FORESY, FORESS 
WRITE (6,2004) ALT1,ALT2,ALT3,ALT4,ALTS,ALT6,ALT7,ALTS,ALT9,ALT10 
WEITE(6,4033) 
4033 FCRMAT(1HO,52X,*FIRE IN FUELS  SUMMABY*) 
WRITE (6,4034) 
8034 FOBMAT(1HO,34X,"FUEFL TYPE*,7X,*NUMBER OF FIRES*,8X,*PROBABILITY® 
1) 
WRITE (6,9035) AIPUEL (1) , PGREAT (1) 
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1442 

1443 

WHY 

1445 

1446 (e FORTRAN IV G MAIN PBOGRAM 

1447 

1448 

1449 4035 FCEMAT(1HO0,33X,*CCNIFEROUS',13X,F8.0, 14X,P6. 2) 
1450 WRITE (€,4036) AIFUEL (2) , PGREAT (2) 
1451 4036 FCRMAT(1HO,33X,"MIXED WOOD',13X,F8.0, 14X,F6. 2) 
14£2 WRITE (6,4037) AIFUEL (3) , PGREAT (3) 

1463 4037 FORMAT (1HO,33X,'LICHEN - MUSKEG*,8X,F8.0, 14X,F6. 2) 
1454 WRITE (€,4038) AIFUEL (4) , PGREAT (4) 

1455 4038 FCKMAT(1HO,33X%,"GRASS",18X,F8.0,14X,F6.2) 

1456 WHITE (6,4039) ICTALG,GREAT 

1457 4039 FCRMAT(1HO,47X," TOTAL", 3X,F9.0,6X, ‘TOTAL’, 2X,F7. 2) 
1458 BEITE(6,5550) 

1459 5550 FCERAT(////1H0,43X,*STATISTICS FOR THIS RUN") 
1469 WRITE (€,5549) XPBAR 

1461 5549 FCRMAT(1HO,13X,"XBAR*,10X,*=",3X,F13.0,8X,'= *,*AVEBAGE TOTAL SEAS 
1462 VONA LCOS 1") 

1463 WRITE (6,5548) VAR 

1464 5548 FCRKAT(1HO,13X,"VARIANCE® ,6X,°=",3X,F 13.0) 

1465 WRITE (6,5547) SIGMA 

1466 5547 FCEMAT(1H0,13X,"STD DEVIATION’, 1X, *=",3X,F 13.0) 
1467 WRITE (€,5546) CD,ZALPHA 

1468 5546 FCRMAT(1HO,13X,"LENGTH OF CONFIDENCE INTERVAL =*,3X,F12.0,5X, ‘AT 
1469 1Z =", F724) 

1470 STCP 

1471 EXC 

1472 

1473 

1474 c THIS SECTION IS THB RANDOS NUMBER GENERATOR 

1475 

1476 SUBFCUTINE GREG 

1477 CCHMCN IX,IY,YFL 

1478 LY=1X*€5539 

1479 IF (1Y)5,6,6 

1480 5 LY=1Y+#2147483647+1 

1481 6 YFL=I1Y 

1482 YFL=YFL*.46566135E-9 

1483 Ix = -IY 

1484 FETUES 

1485 ENC 


END OF FILE 
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The N value in the stopping rule equation, 


Average seasonal total cost of suppression + detection + 


damage, 


Same as AROUTE except for an individual season only, 


The spread index versus buildup index table for the spring 


and summer. 
columns are the 5 buildup index categories. 
in the matrix are the 5 burning index categories. 


The rows are the 4 spread categories and 


the 


The. values 


The spread index versus buildup index table for the fall. 
The rows are the 4 spread categories and the columns are 


the 5 buildup index categories. The 5 burning index 


categories are in the matrix. 
The seasonal cost of suppression. 
The buildup index for the day. 


The time it takes 
location. 


to fly from the airport to the fire 
The burning index for the day as a function of spread 
and buildup. 


A vector of average suppression costs for each of the 
cost categories used. 


The cause of the fire starting. 


The time of day the aircraft arrives at the spot the 
fire occurs. 


The average yearly cost to operate a lookout tower. 
The number of fires falling in each cost category. 
The total suppression dollars for each cost category, 


An intermediate variable calculated to identify the 
upper range of SPOT data. 


An intermediate variable calculated to identify the 
lower range of SPOT data. 
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ELE 


D A variable in the stopping rule called the confidence 
interval. 

DS The cost of damage for a season, 

DAYSA The totalizer of the cost of aircraft detection for the 


total-run. 


DELTAX The horizontal distance from the fire to the air route 
on the map. 


DELTAY The vertical distance from the fire to the air route 
on the map. 


DESIGN The average total detection cost. 
DEST. The ranger district being considered. 
DRANGE The number of fires in each cost category for a season: 


FIDAYS The daily fire count for all the days in a season: 


FIRE The fire count for the number of fires in a season. 
FIREX See MCOORD (2), 

FIREY See MCCORD (1) - 

FORES1 - 


FORES2 The name of the forest. This is printed at the top of 
each page of output. 


FUEL The fuel the fire is burning in. 

GREAT Total percentage of fires in each fuel: i.e. 100Z%- 

GREG The subroutine that calculates a random number between 0 and 1: 
GTS The seasonal cost of suppression for each category + 

HORZN The visibility in miles. 

I An index used to categorize the fire sizes, 

IB This is the minimum value of the burning index that 


initiates aircraft flying , 


IBOB A check to ensure the aircraft fly even on days with no fires. 
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IBOYD 


IBURNA 
IBURNP 
IBURNT 
Ic 
ICLOCK 
IDC 


IDELT 


IDIST 


IFIDA 


IFIRE 


IGROW 


TIX 
IMAXIE 


INK 


IQUAD 


IRANK 


ISEX 
ISIGHT 


ISTAR 


A vector of 
route used, 


The elapsed 
The elapsed 
The elapsed 
This is the 
See CLOCK, 

This is the 


The elapsed 
used by the 


The recoded 


na 


Chespemonths filled with lyor 0. If 1 
if O not used, 


time from fire start until an aircraft discovery, 
time from fire start until a public discovery, 
time from fire start until a tower discovery. 


value of the cause that initiates aircraft, 


district the fire is located in. 


time between fire start and discovery as 
growth model. 


value of DIST so the numbers run from 1 to 5 


for indexing. 


The variable that takes on the value of the number of 
fires for each day. 


Integer value of FIRE. 


An intermediate variable calculated to determine the 
value of TTTIME. 


See AIX. 


The current 


air route leg being flown. 


A counter used to ensure all the towers are used in 
searching for the fire. 


The tower quadrant the fire is located in. 


An intermediate variable used in determining the value 


of TOTAL. 


The maximum 


See SIGHT. 


number of fire seasons the simulation will run. 


An intermediate variable calculated to determine the 
value of TTTIME . 
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The tower code number, 


The minimum angle between an air route leg and the X 
axis (horizontal line). 


See TOWER. 

The value used to initiate the random variable generator GREG, 
The X coordinate of the Hower! 

An intermediate variable used in the random number generator. 
The Y coordinate of the tower, 


The index calculated to check to lower X limit of the 
air route line. 


An intermediate variable created to locate'the correct 
AIRDIS for each leg. 


An intermediate variable created to locate the correct 
AIRDIS for each route. 


An intermediate variable calculated to locate the correct 
row in the RING MATRIX. 


See ISEX, 


The index calculated to check the upper X limit of the 
air route line, 


The current air route being tested. 
A counter used to aid in reading in RING data, 


An index counter used to check that all the air routes 
are flown during the day, 


An index that takes on the value of 1 or 2 to aid in 
determining SPOT, 


The instantaneous counter of the day of the season. 
The total counter for the number of days in a season. 


An index counter used to determine the angle of the 
air route leg. 
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LLL 


MARK 


MATT 


MAXIE 
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MCOORD (1) 
MCOORD (2) 
MCOUNT 


MDIANE 


NAM 


NDAYS 
NDEX 


NREPET 
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The distance from a tower to the fire. 


An intermediate variable calculated to aid in determining 
TOTAL. 


The index calculated to check to lower Y limit of the 
air route line. 


An intermediate variable used to calculate the cause 
and fuel values. 


The flying code criteria used to determine when aircraft 
fly. If MATT = 0, check the values of IB and IC; if MATT 
= 1, fly every day; if MATT = 2, never fly. 

The number of legs on an air route. 

The code for the month from 1 to 5. 

The X coordinate of the fire location. 

The Y coordinate of the fire location. 


The month counter. 


An intermediate variable used in determining the correct 
row in the PCDIST matrix. 


An intermediate variable used to limit in the determination 
of fire location. 


A seasonal counter to keep track of the number of seasons run. 
The month. 


The number of aircraft required to fly all the air routes 
at the designed times. Used in print out only. 


The index calculated to check the upper Y limit on the 
air route line. 


The number of days in the fire season. 
The number of the current fire season run. 
Correction variable to offset the double counting of fires 


on air routes. When NROUTE = 5, NREPET = 1, When NROUTE = 10, 
NREPET = 2. When NROUTE = 15, NREPET = 3. 
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NROUTE 


NROW 


NTOWER 


PB1-PB5 


PCDIST 


PCOST 


PEFFCT 


PFDIST 


PFIRE 


PGREAT 


PISS 


PNFD 


PNFIRE 


PRANGE 


PRCS 


PRV 


PS1-Ps4 


PSIZE 


PSNFA 


dh 


The number of aircraft routes in the design. 


A variable calculated to determine which row to use in 
the SPOT data. 


The number of towers in the design. 

The cumulative probabilities of each of the buildup indices. 
A matrix of the probability data for fire cause given 

the district and month. The rows are the districts for 
each month and the columns are the cumulative probability 


data for the 10 cause categories. 


A matrix of the cumulative suppression probabilities for 
each of the 15 size categories. 


Effective percent of fires discovered by air routes. 
This is the fuel probability data. The matrix has the 
fuel type for each district in each month as the rows and 


the cumulative probability for the fuels as the columns. 


A vector of the cumulative probabilities for a fire for 
each day of the fire season. 


Percentage of fires occurring in each fuel. 
The average fire size at discovery excluding May fires. 
Percentage of tower-discovery fires discovered by each tower. 


This is the cumulative probability for the number of fires 
occurring during the fire season. 


The probability of a fire falling in each cost category. 


The probability of a fire being caused by each of the 
cause categories. 


The probability of a day having each of the visibility 
categories. 


The cumulative probabilities of the spread indices. 
The probability of a fire occurring in each size category. 


Percent of fires falling on each air route. 
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PSNFA2 Percent of aircraft-discovery fires falling on each 
air route, 


PTFD Percent of total fires discovered by each tower. 

PTFDA Percent of total fires falling on each air route. 

PTFDA2 Effective percent of total fires falling on each air route. 
PT IME A matrix of the probability data for the starting time 


of each fire given the cause of fire. The rows are the 
10 causes and the columns are the 25 time categories. 


PUBLIC A matrix of the 4 different fuel categories versus the 
cumulative time probabilities for each of 25 time cate- 
gories the public requires to detect a fire. 


RADIAN The value of theta in radians. 


RANGE The total number of fires over the entire run. 

RES The number of fires falling in each cause category. 

REST An accumulator for fires in each size category not in May. 
RING This is the blind area data for each tower. It is in 


the form of a matrix. 


ROUTE A matrix of the air route data for each air route. The 
matrix will accomodate 15 air routes with a maximum of 
15 legs on each route. 


RSIZE The number of fires in each size category. 

RV The number of days that fall into each visibility category. 
SIGHT The perpendicular distance from a fire to an air route. 
SIGMA The standard deviation about the total mean cost. 

SNFP The total number of fires discovered by the public. 

SPCOST The average cost of suppression. 

SPOE A matrix of coordinates describing the dimensions of 


the 5 districts of the Footner Forest. 
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SPRCS 
SPRDEX 


SPREAD 


SPRING 


SPRV 


SPSIZE 
SRCS 
SRSIZE 
SRV 
STS 
SUM 
SUMMER 
SUMMY 


SUMPR 


TACRES 
TANAY 
T$ 
TEAD 
TIME 


TOTAL 


TOTALG 


PLE 


The total probability of each cause category; ise: 1.00 


The daily spread index.- 


The cubic matrix for the C values of the fire growth model. 
The I values are the 4 fuels. The J values are the 5 
burning indices and the K values are the 25 time categories. 


The number of air routes flown in the spring. 


The total probability of all of the visibility 
categories; i.e. 1.00. 


The 


The 


The 


The 


The 


The 


The 


The 


total 


total 


total 


total 


total 


total 


probability of all the size categories; i.e. 1.00. 
number of fires during the entire run. 

number of fires in all size categories. 

number of days in the entire run. 

suppression dollars spent during the entire run. 


of the frequency distribution when TIME was unknown. 


number of summer air routes being used. 


total 


percent of all fires discovered; i.e. 100%. 


The total probability of fires falling in each cost 
category; 


ie.e 1.00. 


The total acres burnt during the entire run. 


An accumulator that sums the AMAY values. 


The cost of operating the tower for a season. 


The total number of fires discovered by aircraft. 


The time of day the fire starts. 


The length of the discovery leg in the district in which 
the fire occurred. 


Total number of fires in the 4 fuel categories. 
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TOWER The tower that finds the fire. 


TOWERI - 
TOWER5 This is the name of the tower; i.e. Yates. 


TPNFD The. total percent of tower discoveries made by towers; 
176 OOK: 


TPSNFA The total percent of fires on air routes that were 
located on routes; i.e. 1002. 


TPSNF2 The total percent of aircraft discovery fires falling 
On air routes: ise. LOOZ. 


TPTFD The total percent of fires discovered by towers. 
TPTFD2 The percent of all fires that were discovered by aircraft. 
TPTFDA The percent of all fires that fell on air routes. 


TPTFDP Percent of total fires discovered by the public. 


TRESE An accumulator that sums the REST values. 

TSNFA The total number of fires starting on air routes. 
TSNFD The total number of fires discovered by towers. 
TTACRE The total acres burnt during each season. 


TTTIME The hours through which the growth model operates. 


VAR The variance of the total cost distribution over the 
total run, 


VARR An intermediate variable used to determine the value of VAR. 


VIS The visibility data given the spread index for each day of 
the fire season. The rows are the spread index categories 
and the columns are the cumulative probabilities for the 
6 visibility categories. 


VVAR An intermediate variable used to determine the value of VAR. 
WDEX See NDEX. 
WHERE A matrix of probability data for the. district the fire is 


located in. The rows are the months, the columns the 
cumulative probability for each district. 
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See NTOWER. 

The total costs for each season. 

The mean of the total cost distribution. 
A daily counter internal to each season. 


The X value of the route line when the Y value is 
coincident with the Y coordinate of the fire. 


The percentage of frequencies for each of the time 
categories. 


The summation of the percentages of the time categories 
equals 100%. 


The summation of all the total seasonal costs. 
The value of the random variable being generated. 


The Y value of the route line when the X value is coinci- 
dent with the X coordinate of the fire. 


The percentile of the normal distribution which leaves 
alpha/2 % probability in each tail. 


See NTOWER. 
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APPENDIX D 


A LISTING OF THE DATA 


12% 


SSIGNON GAVM *LIST! 


OW AT 


23:12.43 OW 04-10-72 LAST ON AT 09:318.33 


$LIST *SOURCE* 


= 
owoananw 


DATA DECK SET UP FORTRAN IV G 


THIS IS A LISTING OP THE DATA THAT WAS USED TO DRIVE THE MODEL. 
UNDER EACH VARIABLE NAME IS THE NUMERICAL VALUES THAT WERE USED FOR 
THE FOCTNER FPOREST.THEY ARE SHOWN HERE EXACTLY THE WAY THEY APPEAR 
ON THE DATA CARDS. 


NREPET 
1 
NAIR 
2 
ZALPHA 
09500 
IX 
76231491 
CAU 
25 88 268 742 2059 6956 14650 34094 60666 87857285294 
COSTT 
6000 
SPRING, SUMMER 
14 
ALT 1------- were nn nnn nen - === ALT10 
19 TOWERS - FLY EVERY DAY 
FORES 1------------ 2--------PORESS 
POOTNER 
PNPIRE 
33115 
MATT 
1 
Ic 
2 
IB 
2 
WROUTE 
5 
HTOWER 
19 
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DATA DECK SET UP PYORTRAW IV G 


ISEX 
0200 


NDAYS 
183 


AMONIQ 
Zo 48 268 742 2059 6956 14650 34094 60666 87857 285294 


PCOST 
203 428 698 820 987100020003000400050006000 
152 239 478 630 967 98910002000300040005000 
087 130 269 303 825 868 9551000200030004000 
227 272 317 453 90810002000 3000400050006090 
001 125 250 251 875100020003000400050006000 
042 043 166 251 834 95910002000300040005000 
333 334 SCO 501 833100020003000400050006000 
333 334 335 6671000200030004000500060007000 
001 002 143 4829 715100020003000400050006000 
001 250 500 50110002000 30004000500060007000 
040 080 200 360 76C 880 881 920 921 9601000 
040 160 400 401 720 880 9601000200030004000 
007 087 174 348 739 740 870 957 958 9591000 
001 002 063 188 438 563 751 939 940 9411000 
001 C02 003 004 005 006 007 008 009 5001000 


ROUTE 
6410000249413120 48 7 
054000022000000000000000000009500007040080600000000000000000000000 1890000090 
26000000000000000000000000 30 3419 40000000000000000000000000- 67 9343 
24 10 +082-0016 9 37 ¢005¢ 6533 
24 4165-0064 33 


10610318110718110718 1099 1810991580901580901 30063 1300630960680 96068126089 
126089166 1030000000000000000000000000000000000000600000000000000000000000000 
00000000000000000000000000000000000 ; 
210111141714120 4810 

0000000430000000000000000000261- 17360027000000000000000000000000 151- 7336 


65000000000000000000000000 55 00 31000000000000000000000000 219 90 
45 +086 00 15 +173 90 
31 +071 00 30 #195 90 
60 #014¢001207 38 38 +306-021465 


102110129069129069149040 1490402 180402140402 18086214086 166086166086 166070 
1660702000702000 7020004 1200041142034 142034 1021100000000000000000000000000 
000000000000000000000000000000000000 
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122 
123 
124 
W258 
126 
127 
128 
129 
130 
131 
132 
133 
134 
135 
136 
a) iz/ 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
asa 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 
181 
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DATA DECK SET UP FORTRAN IV G 


220111127411120 48 7 

0000000140000000000000000000078 19110000000260000000080000000000168- 7336 

00000 0030000000005000000000006200000 9000000004 800000000000000000000900000000 
30 5 +015 90 58 #121+003821 
18 2) aS #216-011549 


102110090 1050901050601 260601260600 900600900 14090014090014127014127068145 

068145102110000000000000000000000000000000000000000000000000000000000000 

000000000000000000000000000000000000 

230111136411220 48 8 

0060000 13000037000000000000-058 1485603400000000002900000000000001520000000 

010000000000000000000000000020600 0009005 1000000000000000000000000 01390000000 
14 +150 90 63 #127 00 
15 #218 90 82 24 #101+000503 


102110129150129 150193150 193 15019314119319114114 114914114112 7141127207127 
207127207114207114102110000000000000000000000000000000000000000000000000 
00000000000000000000000000000000000 

240111135014120 48 9 

00000 000000000000000 35000000 16000019110000000000000000470100000000 149000 3620 
00000600000000000 1602800000000 14000009000000000000000000005600000001990000000 


UB +239-006232 74 #192 00 
19 #168 90 55 #171 00 


9 + 309-013353 


10117906617 106613710 121540 12154012197012197067197067197077190077190152190 
152190152172 152172106172 106172101179000000000000000000000000000000000000 
000000000000000000000000000000000000 


PUBLIC 


3790 4600 £410 676C 7300 7301 7302 7570 7840 8110 8111 8380 8381 8382 8650 8651 


8920 8921 8922 8923 9190 9460 9461 973010000 


8197 6133 7744 8388 9354 9355 9356 9357 9358 9359 9360 9361 9362 9363 9364 9365 


9366 9367 9368 9677 9678 9679 9680 968110000 


2382 2859 3814 4769 5724 5725 5726 5727 5728 5729 5730 6200 6201 6202 6203 6204 


6205 6206 6207 6676 6677 7630 7631 953210000 


5292 7644 7645 7646 7647 8233 8234 8235 8236 8237 8238 8239 8240 8241 8242 8243 


8244 8245 8246 8247 8248 8249 8822 941110000 
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182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
M29 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
210 
211 
242 
203 
214 
215 
216 
217 
218 
219 
220 
7472.) 
222 
223 
224 
2205 
226 
227 
228 
229 
230 
231 
743322 
233 
234 
235 
236 
Pes | 
238 
239 
240 
241 


DATA DECK SET UP 


TYATES 

TSTEEN FIVER 
6AMBER 

6ADAIR 

3ZAMA 

3BASSET 

3WATT MOUNTAIN 
JRAINBOW 
WHITE SANDS 
1PONTON 

1FOGGY MOUKTAT 
JLAMBERT CREEK 
1BUFFALC 
TWADLIN 

OPANNY 

S9TALBOT LAKE 
SBISCN LAKE 
QHAWK HILLS 
QSJEAN LAKE 


11 50000 
11100000 
11150000 
11200000 
11250000 
12 50000 
12100000 
12150000 
12200000 
12250000 
13 50000 
13100000 
1315000C 
13200073 
13250238 
14 50060 
14100000 
1415000C 
14200000 
14250000 
21 50204 
21100646 
21150388 
21200131 
21250164 
22 50000 
22100239 
221506864 


IDC, TOWER1,ETC 


RING 


1131197 
77180 
19148 
55147 
30108 
64 80 
91111 
11 90 
9184165 
10137131 
11181794 
12206 70 
13137 64 
14159 52 
15212 10 
16150 22 
17118 «#5 
18 92 43 
19217 32 


DBNHNWLSWN 


FORTRAN IV 


124 


. » © | ua ae : 


ee 
Tet 
7 mm : 
Dee 
j 7 


z" en 
| 


ot aos 


c 


242 


251 


256 
257 
258 
259 


262 
263 


284 
285 


288 
289 


301 


22200815 
22250755 
23 5000¢ 
23100408 
23150612 
23200378 
23250272 
24 50593 
24100556 
2415000C€ 
24200000 
24250000 
31 50286 
31100102 
31150102 
31200036 
31250196 
3250255 
32100409 
32150429 
32200679 
3225100C 
33 50000 
33100000 
33150000 
33206000 
3325C00C 
34 50000 
34100000 
3415000C 
34200060 
342500CC 
41 500CC 
41100000 
41150368 
41200254 
41250294 
§2 500CC 
42100000 
4215000C¢ 
42200204 
42250385 
43 50000 
43100000 
43150327 
43200319 
43250181 
44 50000 
44100000 
44150265 


DATA DECK SET UP 


PORTRAN IV 


G 


125 


302 
303 
304 
305 
306 
307 
308 
309 
310 
311 
312 
33 
314 
315 
316 
317 


Shite 


319 
320 
321 
3} 7472 
323 
324 
325 
326 
327 
328 
329 
330 
331 
332 
333 
334 
335 
336 
337 
338 
339 
340 
341 
342 


343 - 
344 


345 
346 
347 
348 
349 
350 
351 
352 
353 
354 
355 
356 
357 
358 
359 
360 
361 


44200291 
44250294 
51 SOOCC 
51100000 
51150000 
51200102 
51250113 
52 5C0C0 
5210C068 
52150082 
5220CC58@ 
S225 0051 
53 50000 
$3100136 
$3150122 
§3200123 
53250113 
54 50000 
54100000 
5415C00C 
5420C000 
5425000C 
61 50120 
61100273 
61150000 
61200000 
6125000C 
62 50000 
62100000 
62150048 
6220CC95 
62250068 
63 5000C 
631C€0000 
63150500 
63200760 
63250500 
64 SOCOC 
64100000 
€415C000C 
64200000 
6425000C 
71 50613 
71100727 
71150674 
71200393 
71250600 
72 50306 
72100375 
72150000 


DATA 


DECK SET UP 


FORTRAN IV G 


126 


362 
363 
364 
365 
366 
367 
368 
369 
370 
371 
372 
373 
374 
375 
376 
377 
378 
379 
380 
381 
382 
383 
384 
385 
386 
387 
388 
389 
390 
391 
392 
393 
394 
395 
396 
397 
398 
399 
400 
401 
402 
403 
404 
405 
406 
407 
408 
409 
410 
411 
412 
413 
414 
415 
416 
417 
418 
419 
420 
421 


72200000 
7225C0CC 
73 50407 
73100153 
73150000 
73200000 
7325000C 
74 50816 
74100272 
7415C000 
7420C000 
7425C00C 
81 50407 
€1100000 
811501€4 
81200181 
81250045 
82 50204 
82100204 
82150250 
82200500 
82250500 
83 5C00C 
83100545 
83150250 
8320C000 
€325000C6 
€4 50203 
84100340 
8415C00C 
84200000 
8425000C 
91 £5000C 
91100000 
91150500 
91200333 
91251000 
92 50000 
9210000 
9215000C 
9220C000 
92250000 
93 50000 
931000CO 
931500060 
93200145 
93250475 
94 scooc 
94100000 
$4150C8&2 


DATA DECK SET UP 


PORTRAW IV G 


127 


422 
423 
424 
425 
426 
427 
428 
429 
430 
431 
432 
433 
434 
435 
436 
437 
438 
439 
440 
441 
442 
443 
YRY 
4uS 
446 
447 
4u8 
4ug 
450 
451 
452 
453 
454 
455 
456 
457 
458 
459 
460 
464 
462 
463 
464 
465 
466 
467 
468 
469 
470 
471 
472 
473 
474 
475 
476 
477 
478 
Q79 
480 
ast 


94200043 
94250306 
101 50000 
1071100000 
101150000 
101200218 
101250362 
102 50510 
102100693 
102150775 
102200635 
102250351 
103 50900 
103100409 
103150000 
103200000 
10325c0C0 
104050407 
104100625 
104150204 
104200276 
104250113 
111850795 
1111008€4 
111150886 
11120Ce55 
111250835 
112 5000C 
112100102 
112150530 
11220C680 
112251000 
113 50000 
11310C000 
113150123 
113200CE7 
113250068 
114 50204 
114100000 
114150000 
11420000C 
121425C00C 
121 50000 
1211000C0 
121150000 
1212000060 
121250000 
122 50000 
122100000 
12215C000 


DATA DECK SET OP 


FORTRAN IV G 


128 


482 
483 
484 
485 
486 
487 
488 
489 
490 
491 
492 
493 
494 
495 
496 
497 
498 
499 
500 
501 
$02 
503 
504 
505 
506 
507 
508 
509 
510 
511 
Siz 
513 
514 
515 
516 
517 
518 
519 
520 
sZt 


522, 


523 
524 
525 
526 
yr) 
528 
529 
530 
531 
58) 72 
533 
534 
535 
536 
537 
538 
539 
540 
541 


122200C87 
122250136 
123 50000 
123100000 
123150307 
1232C07S€ 
123251000 
124 5000C 
12410000C 
12415C00C 
1242000C0 
124250000 
131 50000 
131100000 
131150000 
131200000 
1312560C0 
132 50000 
132100000 
132150000 
132200000 
13225C0C0 
133 50000 
133100728 
133150755 
133200766 
133250750 
134 59000 
134100443 
134150673 
134200826 
134250886 
141 5090C 
141700340 
1411500C0 
1412000C0 
1412500C0 
142 5038€€ 
142100795 
142150450 
1422002C3 
142250153 
143 500C0 
143100136 
143150367 
143200752 
143250637 
144 50500 
144100136 
146415000C 


DATA 


DECK SET UP 


PORTRAW 


Iv G 


129 


542 
543 
544 
545 
546 
547 
548 
549 
550 
551 
5512 
553 
554 
S55 
556 
557 
558 
559 
560 
561 
562 
563 
564 
565 
566 
567 
568 
569 
570 
S71 
572 
573 
574 
575 
576 
577 
578 
579 
580 
581 
582 
583 
584 
585 
586 
587 
588 
589 
590 
591 
DIZ 
593 
594 
$95 
596 
597 
$98 
599 
690 
601 


144200000 
14425COCC 
151 510006 
151101000 
16115C00C 
15120C00C 
15125010C 
152 51000 
152101000 
15215CCCO 
15220C400 
152250700 
153 5100¢ 
153101000 
153151000 
153201000 
153251000 
154 5100C 
154101000 
15415100C 
154201000 
15425100C€ 
161 5000C 
161100CO0C 
16115000G 
1612000C0 
161250020 
162 50400 
162100500 
162150800 
162200950 
16225C8C0 
163 51¢0C0 
163101000 
163151000 
163201000 
16325100C 
164 5100C 
164101006C 
1641510CC 
16420100C 
1€42510CC 
171 51000 
171101000 
171151000 
17120C000 
17125C0C0 
172 51000 
1721010060 
172151000 


DATA DECK SET UP 


FORTRAN IV G 


130 


172200000 
17225C000 
WATSIOCO 
17310 1000 
17315100C 
173201000 
173251000 
174 §10C0 
174101000 
17415100C 
174201000 
1742510CC 
181 51000 
1€1101000 
186115C0CC 
181200000 
181250000 
182 5100C 
182101000 
1€2151C000 
182201000 
182251000 
183 510C0 
183101000 
18315100C 
1€32010CC 
183251000 
184 §1000 
18410100C 
18€415C000 
18420000C 
1@425C0CC 
191 51000 
161101000 
1911510CC 
191201000 
1912510€C 
192 51000 
19210C0CO 
1921500090 
192200100 
192250330 
193 51000 
19310C0C% 
193150059 
19320C0CO 
193250C00 
194 51000 
194101000 
194151000 
19420100¢ 
1942510CC 


DATA DECK SET UP 


FORTRAN IV G 


134s 


662 
663 
664 
665 
666 
667 
668 
669 
670 
671 
672 
673 
674 
675 
676 
677 
678 
679 
680 
681 
682 
683 
684 
685 
686 
687 
688 
689 
690 
691 
692 
693 
694 
695 
696 
697 
698 
699 
700 
701 
702 
703 
704 
705 
706 
707 
708 
709 
710 
711 
712 
713 
714 
715 
716 
717 
718 
719 
720 
721 


DATA DECK SET UP FORTRAN IV G 


PFDIST 
4400 6300 630110000 
3070 6140 845010000 
5000 5001 500210000 
2500 2501 500010000 
8000 80011000020000 
5000 6750 925010000 
3100 7940 966010000 
417C 5C00 917010000 
4122 €472 941210000 
3620 4540 970010000 
§360 6074 921410000 
7850 8&5651000020000 
875010C002000930000 
3570 50C5 $28510000 
6000 6C011000020000 
3€40 9C90 909110000 
00000 5560 890010000 
3751 7502 906210000 
3751 7502 904210000 
6250 75001000020000 
00000 e€5711000020000 
00000 &£711000020000 
00000 &5711000020000 
00000 85711000020000 
00000 85711000020000 


MONTH, NCODE,ETC 
MAY 01 71714 &5710002000 9291000200030004000 
MAY 01 2 667 73410002000 9331000200030004000 
MAY 01 3 £56 557 8341000 9441000200030004000 
MAY 01 4 526 527 6311000 7891000200030004000 
BAY 01 5 565 652 9561000 7831000200030004000 
BAY 01 6 619 714 8571000 762 905100020003000 
BAY 01 7 667 668 8751000 792 917100020003000 
MAY 01 8 679 680 8581000 750 964100020003000 
BAY 01 9 €55 724 9311000 690 966100020003000 


MAY 0110 515 545 7871000 667 970100020003000 
MAY 01117 444 583 8331000 694 944100020003000 
MAY 0112 417 473 7791000 667 917100020003000 
MAY 0113 703 757 8651000 676 865100020003000 
MAY 0114 622 838 9461000 649 865100020003000 
MAY 0115 526 684 9211000 632 869100020003000 
MAY 0116 342 395 6841000 605 789100020003000 
BAY 0117 306 487 8721000 513 795100020003000 
BAY 0118 564 641 7441000 590 821100020003000 
BAY 0119 4OS 513 6481000 486 810100020003000 
MAY 0120 350 575 7251000 400 825100020C03000 
MAY 0121 282 538 7941000 333 769100020003000 
MAY 0122 225 425 9001000 200 675 97510002000 


1372 


722 
723 
724 
725 
726 
727 
728 
729 
730 
731 
732 
733 
734 
735 
736 
msi7/ 
738 
739 
740 
741 
742 
743 
7T4y 
745 
746 
747 
748 
749 
750 
751 
752 
753 
754 
755 
756 
757 
758 
759 
760 
761 
762 
763 
764 
765 
766 
767 
768 
769 
770 
771 
va2 
773 
774 
775 
776 
777 
778 
779 
780 
781 


BAY 0123 
BAY 0124 
BAY 0125 
BAY 0126 
NAY 0127 
BAY O12€ 
MAY 0129 
MAY 0130 
MAY 0131 
JUNEO2 
JUNEC2 
JUNEO2 
JUNEOD2 
JUNEO2 
JUNEO2 
JUNEO2 
JUKFO2 
JUNEO2 
JUNEOZ10 
JUNEONZ11 
JUNEO212 
JUNEOZ13 
JUNEOZ14 
JUNEO215 
JUNEO216 
JUNEO217 
JUNEO218 
JUNEQ219 
JUNEN220 
JUNEO 221 
JUNEO222 
JUNEN223 
JUNEO224 
JUNEO22§ 
JUNEOZ26 
JUNEO227 
JUNEO228 
JUNEO229 
JUNEO230 
JULYO3 1 
JULYO3 
JULYO3 
JULYO3 
JULYO3 
JULYO3 
JULYO3 
JULYO3 
JULYO3 
JULYO310 
JOLY0311 


WDMDNIAMNE WA = 


WO WIAME WA 


175 
167 
Sieg 
571 
£24 
429 
429 
1990 
071 
091 
114 
227 
295 
295 
341 
114 
136 
273 
140 
250 
182 
227 
114 
W222 
220 
262 
415 
190 
190 
PAKS, 
381 
395 
186 
214 
190 
41S 
279 
186 
146 
140 
233 
364 
318 
341 
116 
250 
386 
295 
318 
369 


CATA DECK SET UP 


325 8501000 
405 7861000 
571 7851000 
714 9521000 
714 8811000 
691 9051000 
691 9051000 
595 9051000 
357 7861000 
296 8311000 
387 8421000 
727 9321000 
590 9311000 
772 9311000 
705 9101000 
523 8871000 
681 9311000 
88710002000 
698 9071000 
72710002000 
€37 $051000 
79510002000 
263 8641000 
390 8781000 
708 9031000 
810 9051000 
781 9761000 
690 9521000 
690 88C10C0 
£14 9771000 
€3310002000 
79010062000 
674 9761000 
571 9281000 
666 $761000 
768 $771000 
9€710002000 
651 9301000 
702 $511000 
745 9781000 
74510002000 
£1910002000 
79510002000 
659 9771000 
628 8841000 
864 10002000 
886 9771000 
79510002000 
75010002000 
728 9781000 


vas) 
190 
143 
214 
238 
214 
190 
190 
167 
205 
182 
159 
182 
159 
159 
159 
136 
159 
140 
136 
114 
114 
114 
098 
146 
019 
122 
119 
119 
093 
095 
116 
116 
095 
095 
116 
163 
163 
122 
093 
093 
114 
136 
136 
116 
136 
159 
159 
205 
182 


650 
€42 
548 
STA 
643 
595 
595 
619 
596 
637 
614 
591 
614 
591 
523 
500 
454 
432 
419 
431 
387 
364 
319 
318 
268 
286 
366 
357 
333 
302 
333 
S72 
342 
357 
381 
369 
349 
349 
342 
326 
302 
387 
363 
318 
349 
363 
454 
500 
546 
546 


FORTRAN IV 


97510002000 
952100C€2000 
90510002000 
90410002000 
90510002000 
92810002000 
92810002000 
88110002000 
88210002000 
91010002000 
90910002000 
90910002000 
88710002000 
86410002000 
81810002000 
79510002000 
77210002000 
79610002000 
79110002000 
79510002000 
77310002000 
75010002000 
70510002000 
68410002000 
68310002000 
71510002000 
75610002000 
76210002000 
71410002000 
74410002000 
76210002000 
79110062000 
76710002000 
738 9521000 
714 9521000 
625 9511000 
675 9541000 
744 9531000 
757 9521000 
698 9541000 
628 9541000 
728 9551000 
772°9C02000 
795 9771000 
768 9771000 
81810002000 
84010002000 
81810002000 
88710002000 
€8710002000 


G 


Wee 


aed ds tk Cid 


ny 
os 
ae 


: ae : 


See Cheese aa ecg - 


782 
783 
784 
785 
786 
787 
788 
789 
790 
791 
1912 
793 
794 
795 
796 
797 
798 
799 
800 
801 
802 
803 
804 
805 
806 
807 
808 
809 
810 
811 
812 
813 
814 
815 
816 
817 
818 
819 
820 
821 
822 
823 
824 
825 
826 
827 
828 
829 
830 
831 
832 
833 
834 
835 
836 
837 
838 
839 
840 
841 


JULYO312 
JULYO313 
JULYO314 
JULY0315 
JULYO0316 
JULYO317 
JULYO318 
JULY0319 
JULY0320 
JULYO321 
JULY0N322 
JULY0N323 
JULYO324 
JULYO325 
JULY0326 
JULY0327 
JUDLYO32R 
JULY0329 
JULYO330 
JULYO331 
AUG O04 1 
AUG 04 
AUG 04 3 
AUG 04 4 
AUG 04 5 
2UG O04 6 
AUG O04 7 
AUG 04 @ 
AUG 04 9 
AUG 0410 
AUG 04171 
AO0G 0412 
AUG 0413 
AUG 0414 
AUG 0415 
AUG 0416 
AUG 0417 
AUG 0418 
AUG 0419 
AUG 0420 
AUG 0421 
AUG 0422 
AUG 0423 
AUG 0424 
AUG 0425 
AUG 0426 
AUG 0427 
AUG 0428 
AUG 0429 
AUG 0430 


341 
636 
727 
Sri 
302 
4€5 
341 
209 
2e2h7 
500 
256 
235 
442 
419 
250 
4€8 
500 
341 
341 
386 
477 


BSN 


227 
4S2 
273 
405 
€s9 
256 
190 
227 
5s 
220) 
136 
233 
372 
409 
119 
273 
JE4 
386 
256 
209 
167 
310 
419 
674 
659 
707 
$37 
537 


CATA DECK SET UP 


659 9541000 
931100C2000 
$7710C€C2000 
90410002000 
790 9761000 
907 9771000 
81810002000 
744 9771000 
704 9311000 
€4110002000 
698 $761000 
818 9771000 
761 947100C 
861100C2000 
72710002000 
90710062000 
955 9781000 
773 $781000 
€1810C002000 
79510002000 
€6310002000 
881 9521000 
9321CCC2000 
$C4 10002000 
9321000200C 
929100C 2000 
$77 9781000 
791100C2000 
€42 9521000 
568 9541000 
€82 9551000 
863 9771000 
636 9541000 
75410002200 
907100C€2900 
93210002000 
690100C2000 
773 9781000 
659 9541000 
545 9091000 
768 $771000 
72110002000 
738 9761000 
786 10002000 
861 9771000 
9€7100C02000 
92710002000 
95110002000 
927 9761000 
85410002000 


PAS) 
318 
386 
381 
302 
372 
341 
279 
PLP 
250 
256 
250 
349 
419 
432 
372 
364 
409 
341 
455 
432 
381 
227 
286 
250 
190 
159 
163 
143 
114 
114 
114 
114 
116 
116 
136 


119 


114 
114 
114 
116 
116 
119 
143 
209 
326 
268 
293 
268 
244 


590 
591 
681 
667 
697 
628 
727 
721 
659 
682 
675 
591 
651 
675 
682 


FORTRAN IV 


93210002000 
93210C02000 
93110002000 
93310002000 
93010002000 
93010C02000 
93210C02000 
68110002000 
86410002000 
€87100C2000 
88410002000 
86410002000 
88410002000 
95410002000 
95510002000 


721100020C03000 
705100020C03000 
708100020C03000 


705 
682 
682 
691 
704 
667 
636 
€42 


95510002000 
955100C2000 
93210CC2000 
95310002000 
95410002000 
95310002000 
93110002000 
90410CC2000 
97710002000 
95410C02000 
95210002000 
93210C€C2000 
91010002000 
88610002000 
90910002000 
90810002000 
88610002000 
88610002000 
88110002000 
88710C02000 
91010002000 
84210C02000 
88410002000 
688310CC2000 
881 0002000 
88110C02000 
90710002000 
9€8100C2000 
90210002000 
91710002000 
92710002000 
90310002000 


G 


134 


842 
BuY 
844 
845 
846 
847 
848 
849 
850 
6851 
852 
853 
854 
855 
856 
857 
858 
859 
860 
861 
862 
863 
864 
865 
866 
867 
868 
869 
870 
871 
872 
873 
874 
875 
876 
877 
878 
879 
880 
881 
882 
883 
884 
885 
886 
887 
888 
889 
690 
891 
892 
893 
894 
895 
896 
897 
898 
899 
900 
901 


AUG 0431 436 
SEFTOS 466 
SEPTO5 486 
SFFTOS 595 
SEPTOS 622 
SEPTOS 694 
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APPENDIX E 
COMPUTER OUTPUT GENERATED BY THE MODEL 


The computer print out from a single simulation run 
consists of five pages of output. The first page displays the 
codes used throughout the program. The second page displays the 
relevant cost, size, cause and visibility information. The third 
page illustrates the information generated on the tower system and 
the fourth page shows the information generated on the aircraft and 
air route system. The fifth and last page displays information on 
fuels and the confidence interval for the average total seasonal cost 
as it is calculated by the sequential stopping rule. Each page of 
output was designed for a different application of the illustrated 
information. Table E-1l summarizes the uses of the information on 
each page. 

An example of the output from a single run is shown at the 
end of this appendix. This is an example of a design alternative 
consisting of 19 look-out towers and one aircraft flying the spring 
route and four aircraft flying the summer routes. The aircraft 
operated on days having a burning index of 2 or greater or on days 
having lightening caused fires. The simulation stopped after 50 
seasons as shown at the top of the second page of output. The 
confidence interval shown on the fifth page of output reveals that 


the probability is 95% that the average total seasonal cost for this 
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Statistics For Run 
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TABLE, E-i 


USAGE OF THE OUTPUT 


OUTPUT USAGE 


Codes for Variables used throughout main program. 
Information required for main program deciphering: 


Test the frequency distribution for validity. 
Compare the total costs in each category. 


Test the frequency distribution for validity. 

Presents a probability distribution to operationalize other 
research. 

Compare May fires with fires occurring in other months. 


Test the frequency distribution for validity. 
Compare the number of fires in each category. 


Test the frequency distribution for validity. 

Compare the average total cost of detection with other designs. 
Compare the average total cost of suppression with other designs- 
Compare the average seasonal total cost with other designs. 
Evaluate the number of fires discovered by each tower. 

Compare the number of fires discovered by towers with other 
detection systems. 

Evaluate the number of fires occurring on each air route- 
Evaluate number of fires discovered on each air route. 

Compare the effective percent of fires discovered by aircraft 


with the percent of total fires occurring on air routes. 


Compare the percent of total fires discovered by public with 
the other detection agencies. 


Test the frequency distribution for validity. 


Evaluate the average total seasonal cost and the length 
of the confidence interval. 
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run is $601,826 + $40,360. The confidence interval is partially 
determined by ZALPHA. For ZALPHA equal to .95, the researcher 
requires a 95% confidence level. If ZALPHA was set equal to .99, 
the researcher would require a 99% confidence level and the confi- 


dence interval would be larger. 
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